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DYNAMIC SIMILARITY OF INSTABILITIES OF DISPLACEMENT 
FRONTS IN POROUS MEDIA! 


A. E. SCHEIDEGGER? 


ABSTRACT 


It has been observed that instabilities may occur when a fluid penetrates into 
a porous medium which contains another, more viscous fluid. The problem of 
describing these instabilities analytically is very difficult so that recourse has 
been taken to model experiments. In the present paper, the scaling relationships 
(dynamic similarity conditions) for the interpretation of such model experiments 
are derived, based upon a variety of theories of the instability phenomena. 


1. INTRODUCTION 


It is well known that instabilities may occur when a fluid contained in a 
porous medium (such as rock strata) is displaced by another, less viscous one. 
It has been observed that, under the conditions just mentioned, protuberances 
may occur which may travel through the porous medium at comparatively 
great speed. This phenomenon is called “‘fingering’’. 

The mechanics of finger formation is extremely complicated. Theoretical 
descriptions have been attempted but not one is entirely satisfactory. Generally, 
the displacement equations are so complicated that only approximate solutions 
have been obtained. 

Whenever it is difficult to treat a problem theoretically, one generally 
attempts to investigate it by means of scale-model experiments. In order to 
interpret the observations on scale models properly, the appropriate scaling 
relationships must be known. The latter can be deduced if the basic defining 
equations of the system are known. A knowledge of the specific solutions of 
the equations of motion is not necessary. 

Since there is a variety of descriptions of the mechanics of fingering, one 
will also find a variety of scaling relationships. We shall deduce in this paper 
the scaling relationships that correspond to such theories of the fingering 
phenomenon as have been proposed in the literature. 

‘Manuscript received May 23, 1961. 
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2. THE THEORY OF SCALING 

Prior to an attempt at deriving the scaling relationships appropriate to the 
fingering problem, it might be well to review some of the basic principles of 
the theory of dynamic similarity. 

We shall be concerned with the case where the basic equations of motion 
of the systems under consideration are known, or essentially known. These 
equations of motion will contain a series of variables (some of which can be 
regarded as independent variables, the others as dependent variables) and a 
series of parameters. 

The method for arriving at scaling relationships is to make such a trans- 
formation of variables that the new variables are dimensionless. In this way, 
the equations of motion also become dimensionless. All systems whose 
equations of motion can be transformed into identical dimensionless equa- 
tions are dynamically similar. The scaling conditions are obtained simply 
by setting the expressions for the dimensionless variables as equal to each 
other in the systems under consideration and by expressing the condition that 
the remaining dimensionless parameters in the dimensionless equations of 
motion be equal for the equivalent systems. The latter condition will ensure 
that the dimensionless equations of the equivalent systems are identical. 

An interesting fact in connection with the above procedure can be derived 
from the Pi-theorem. It is well known (‘‘Pi-theorem’’, see, for example, Focken 
1953) that any complete system of equations containing ” physical quantities 
(parameters and variables) with m independent dimensions can be expressed 
in terms of »—m products of these physical quantities. With regard to our 
scaling problem this means that a number of parameters equal to the number 
of independent physical dimensions can be eliminated, i.e. can be absorbed 
into the dimensionless variables. , 

In our calculations, we shall denote the dimensional variables by their 
usual symbols. The corresponding dimensionless variables will be denoted 
by a bar over the symbol in question. The ratios of corresponding variables 
(say, x) in prototype and model will be denoted by R(x) so that 

x (prototype) 


9 2 - oa oe 
(2.1) R(x) x (model) ~* 


It is for these ratios that one wishes to establish relationships. These are the 


scaling conditions. 


3. THE NONLINEAR DISPLACEMENT EQUATIONS OF THE MUSKAT MODEL 

The first model of the explanation of instabilities which we shall consider 
is the Muskat model. This model assumes that there is a sharp front between 
the displacing and the displaced fluid, and that on each side the flow proceeds 
according to Darcy’s law, based on the corresponding viscosity. The Muskat 
model does not specify whether the fluids are miscible or immiscible, since 
only saturations (concentrations) of zero or unity are allowed. 

The Muskat model represents an oversimplification of the actual facts 
inasmuch as some mixing will always occur at the interface between two 
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fluids. Nevertheless, it is known that it yields acceptable results in many 
cases. 

The equations of motion for a two-dimensional displacement process (along 
the x axis) have been investigated by the writer in an earlier paper (Scheidegger 
1960a). They can be deduced as follows. 

Let us assume that at time ¢, the front, if the displacement were stable 
(piston-like), would have reached the position x’ = 0. Let us assume that the 
over-all injection takes place with a constant pore velocity v,. Then we have 
(3.1) x’ = x—dpl. 

Owing to the instability of the front, however, the interface is not given by 
x’ = 0, but by an equation 

(3.2) x = filly) 

where y is a co-ordinate orthogonal to x. Because of the continuity conditions 
and the initial conditions one has 


(3.3a) | fily)dy = 0, 


0 


(3.30) fy, to) = fo(y), 


where w is the width of the system. The pressure pq in the displacing and the 
pressure p; in the displaced fluid must satisfy Laplace equations: 


dpa 4 Opa 


3. = 0, 

oe) Ox ~ oy" 

} a 9° 

(3.5) es [oe 
Ox oy 


Furthermore, continuity at the front requires 
(3.6) pa(interface) = p;(interface) 


and from Darcy's law: 


_ : : k ‘ : 
(3.7) v, interface) = P grad, Pa(interface) 
Ma 


k ; 
2a grad, p;(interface) 
Mt 


where grad, denotes the component normal to the interface of the gradient. 
Furthermore, & denotes the permeability and P the porosity of the medium, 
u is the viscosity of the fluid, and v the velocity of the interface. Finally, in 
order to maintain a constant over-all injection velocity v,, one must have 


; k , k 
(3.8) lim —_ uaP grad pg = lim _ uP grad py = vp. 


z z 
The above set of equations (3.1 to 3.8) constitute the displacement equations. 
The variables are x, x’, y, Pa, P1, and t, and the parameters are vp, k/(uaP), 
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k/(u:P), w, and the function fo(y). There are three independent dimensions. 
According to the predictions of the Pi-theorem, it will be evident that three 
parameters can be eliminated. 

Following the general procedure outlined in Section 2, we make a trans- 
formation of variables so that the new (barred) variables will be dimension- 
less. This may be accomplished as follows by setting: 


(3.9) t= 2/0, x = WE; 

(3.10) Yy = y/w, y= wy; 

(3:11) z’ = x'/w, x’ = wi’: 

(3.12) Pa = —*—b,, a = PatyuawP /k 

P VplawP ’ pM¢ / 

7 : k - ; 

(3.13) pp. = = 5 Pt Pr = PwUpu wl /k: 
Vpu we 

(3.14) D = 9 Up» , = VUp; 

(3.15) i = tv,/w, t = lw/d. 


Then, the system representing the equations of motion becomes 


(3.1’) gf =z-i 
(3.2’) t’ = fi(9), 
a1 
3.3a’) f.(g)dg = 0, 
( J. 
(3.30') FG, bo) = fo/w, 
Opa Oda 
: : ace a 0, 
(3.4°) aR” + ag 
Fe O De Ope 
3.5° es maa eG) 
(3.5") Ox" Og” 
(3.6’) pa(interface) = er Pr(interface), 
Ma 
(3.7’) V(interface) = — erad, pa(interface) 


where the bar over the gradient indicates differentiation with respect to 


barred co-ordinates; 
(3.8) lim — grad pag = lim — grad p; = 1. 


r'’+4—a Z’->c0 


As required by the Pi-theorem, three parameters have disappeared from 
the dimensionless equations. Thus, we have the result that, provided the 
Muskat model be valid, all two-dimensional displacement processes with 
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the same viscosity ratio ua/u, and the same initial conditions are dynamically 
similar to each other. Denoting the ratios of a variable or parameter x in 
the prototype (x,) and model (x,) by R(x) = xp/xm, the scaling conditions 
can be stated as follows: 


(3.9’) R(x) = R(w), 

(3.10’) R(y) = R(w), 

(3.11’) R(x’) = R(w), 

(3.12) R(pa) = R(vp)R(ua)R(w)R(P)/R(R), 
(3.13) R( Pr) = Rv) Ru) R(w)R(P)/R(R), 
(3.14’) R(v) = R(wp), 

(3.15’) R(t) = R(w)/R(vp). 


In addition, as noted above, the viscosity ratios u;/uq must be identical in 
model and prototype; this can be expressed as follows 


(3.16’) R(ut/ua) = 1. 
Similarly, the initial condition must be the same: 
(3.17’) R(fo/w) = 1. 


Conditions (3.9’)-(3.15’), (3.16’), and (3.17’) define the scaling require- 
ments completely. 


4, SPECTRAL ANALYSIS OF INCIPIENT INSTABILITIES 

We shall now consider a different theoretical model for the description of 
fingers in a displacement process. If one assumes the fingers as small (incipient 
fingers), then it is possible to linearize the displacement equations of the last 
section and to describe the front by its spectral components. The equations 
of motion can then be formulated in terms of these spectral components. 
The general theory of this procedure has been outlined by the writer in an 
earlier paper (Scheidegger 1960d). 

Accordingly, we consider again a two-dimensional displacement process 
(parallel to the x axis). At the time ¢, if the displacement were stable, we 
assume that the front would have reached the position x’ = 0. If we surmise 
that the over-all injection takes place with the constant pore velocity v), we 
have again 
(4.1) x’ = x—9,f. 


Again, owing to the instability of the front, the interface is not given by 
x’ = 0 but by an equation 


(4.2a) x’ = fily), 


(4.2) x’(to) = foly), 








1258 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


where y is a co-ordinate orthogonal to x. We now assume that f; is small and 
represent it by its Fourier analysis: 


(4.3) fat J a(K, ¢) sin (Ky)aK + = J b(K, t) cos (Ky)dK 
0 0 


with 


Den 
(4.4a) a(K,t) = J fi(y) sin (Ky)dy, 
ie 
(4.40) b(K,t) = i fily) cos (Ky)dy. 
The equations of motion can then be written as follows: 
(4.5a) g a(K, t) = aga, 
dt 
2 d ‘ 
(4.50) dt b(K, ti) = axb. 


The quantity ax is a function of K which may contain a series of parameters. 

The above equations (4.1-4.5) represent the defining equations of the 
system. The variables are x, y, ¢t, x’, a, b, K, the parameters are vp, do, bo 
(initial values of a, b) and such parameters as may be needed to define ax. 
There are two independent dimensions. 

According to the Pi-theorem, one knows that it is possible to write the 
system of basic equations in dimensionless form and to eliminate two of the 
parameters. The first parameter to be eliminated is v,; as the second we take 
a length parameter w inasmuch as we assume that the function ax depends 
on such a length parameter. It should’ be noted that the width of the system 
does not now enter into the equations as it did in the last section. 

We thus make the following transformation of variables to obtain dimension- 
less variables (the dimensions of the original variables are shown in square 


brackets): 


(4.6a) Cx! fy, oe = 0t! {L]; 
(4.60) ~E = x/w, Gwe: {L]; 
(4.6c) y = y/w, y = wi ies 
(4.7) d = 0/U, Vv = dup, ier 
(4.8) t = tv,/w, t = tw/v, [T]; 
(4.9) f, =f,/w, fh. = uf, [L]; 
(4.10) K = Kw, K = K/w, ie}: 
(4.11a) a =a/w’, a = dw , {L?]; 
(4.110) b = b/w’, b = bw, [L2]; 


(4.12) a = aw/%, a = W,/w, FE). 





4 
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In terms of the above variables, the equations of motion will look as follows: 


(4.1’) £’ = z-i, 
(4.2a’) #’ = f,(9), 
(4.20’) z'(i = 0) = fo/w, 


(4.3’) fr= tf a(R, t) sin (Kg)dK + tf -6(K, D) cos (Ky)dK 


with 
— ea ‘ike » - - 
(4.4a’) a(K, i) = J f(¥) sin (KG)d9, 
- - a Foo = 
(4.40’) 6(K, t) = | f(¥) cos (Kg) dz, 
(4.5a’) = a(K, ft) = a(K)a, 
dt 
ey G2 See 
(4.50’) — b(kK, t) = a(kK)b. 
dt 


The system of equations (4.1’)—(4.5’) represents a dimensionless form of 
the equations of motion. It is free from parameters, but the functions a(K) 
and #(g) must be the same in model and prototype (for all K and g). The 
scaling relationships are accordingly: 


(4.6a’) R(x’) = R(w), 
(4.65’) R(x) = R(w), 
(4.6c’) R(y) = R(w), 

(4.7’) R(v) = R(v,), 
(4.8’) R(t) = R(w)/R(v,), 
(4.9a’) Rf.) = R(w), 
(4.9b’) R(fo) = R(w), 
(4.10’) R(K) = 1/R(w), 
(4.11a’) R(a) = [R(w)), 
(4.110’) R(b) = [R(w))?, 
(4.12’) R(a) = R(vp)/R(w). 


With regard to the last condition, we should like to emphasize that a@ is a 
function of K, and that the condition must be satisfied for all values of K. 
The length ratio R(w) is arbitrary in the above equations. 
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The above scaling relationships are incomplete inasmuch as the form of 
ax is not specified. The best-known expression for ax is one proposed by 
Chuoke et al. (1959); it may be written as follows (neglecting gravity): 

r 73 
(ur—Ma)¥pK —o K 


ke for K < Ky 
- Mtt+ua . 


(4.13) ) ax = 0, forsKk > K~ 


: ee 
K, = “ee (ur—Ha)Up - 


We now see that further parameters have been introduced. In the above 
formulas uw, and yg are the viscosities of the displaced and displacing fluid, 
respectively, and o* is an effective interfacial tension. Instead of the latter, 
it may be preferable to introduce the critical wavelength \, = 27/K, as an 
independent parameter. We then have 





(4.14) CEES ealia A «| ac x] ; for K < 2m/X., 
| ax = 0, for K > 27/d., 

| / 2 
(4.14’) ag = —] a 2s ], for K < 2rw/de, 
ax = 0, for K > 2rw/)¢. 


Since A, is a given quantity in any system, it is obvious that w is no longer 
arbitrary, in fact one ought to set A, = w. With this in mind, the scaling 
conditions (4.6’) to (4.12’) will be completed if one adds thereto 


(4.15) R(w) = R(.), 
(4.16) Rips) pia) = A. 
The conditions (4.6’) to (4.12’) with (4.15) and (4.16) represent the complete 
scaling conditions of the system. 
5. HETEROGENEITIES IN POROUS MEDIA 
The spectral theory of incipient instabilities can be generalized by making 
allowance for the heterogeneities which are present in a porous medium (see 


Scheidegger 1961). The basic equations (4.1)—(4.4) remain the same, but the 
equations of motion then become 


(5.1a) “ a(K,t) = agat+ ¢(t, K), 


q b(K, t) = axb+¢(t, K), 


(5.10) 2 


where ¢(t, A) [L?T-'] is a random function describing the heterogeneities. 
It may also be given by its Fourier transform: 
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(82) o(K,) = + J ” xn(M) sin (Mt)dM+ + wx(M) cos (Mt)dM. 
0 0 


Tv 
In a macroscopically homogeneous medium one has x = w [L?]. Little can 
be done with the general form of the equations of motion as given in (5.1). 
It is useful to make a specific assumption by writing x as follows: 


a) x = AM, Ke Pe PR, 


This introduces two further parameters, viz. A [L?T?] and y [T], which are 
characteristic of the microheterogeneities in the porous medium. Then, the 
equations of motion can be integrated and one obtains (with a general ax): 


2 


: 1 1 2s IS pe eae 
(5.4) a(K) = — Av,Ke™"e" * e**! 
Tv 


x4] ~o%*5i(—aby) [4 —F axe” **[1—erf (ax) 


and a similar expression for b(AK). Thus, eq. (5.4) replaces (4.5). The variables 
are now x, y, t, x’, a, b, K and the parameters vp, Ac, ur/ua (using the expression 
4.14 for ax), as well as A and y. One has, therefore, five parameters; since 
there are two independent dimensions, two of these five parameters can be 
eliminated. We introduce the same dimensionless variables as in Section 4 
(but with A, = w, cf. equations (4.6) to (4.12)). Then the equation of motion 
for d@ becomes (see 5.4): 


Avy ~/7 TVA _ vv, _2 1 Yo. \ 
+ yx ~~ oa yi &q exp (=e ak E — erf (2 ax) || , 


The parameters which occur in this expression are 


ne . UpY 
(5.6) S,= >, 
Xe 
ae Avp 
(o.4) Se = = . 
AcY 


Thus, the scaling relations are exactly those that were found in Section 4, 
but, in addition, the two dimensionless expressions (5.6) and (5.7) also have 
to be identical in model and prototype. 

The occurrence of the parameter y is not very satisfactory, inasmuch as 
it represents a characteristic time which cannot be directly related with any 
properties of the porous medium. It has been shown, however (Scheidegger 
1960), that y can be related to a characteristic wavelength Ay of the hetero- 
geneities in a porous medium by the relation: 


(5.8) lun = V8 @ dy. 
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Introducing (5.8) into (5.6) and (5.7), the latter equations are changed to read 


-~ar | An 

5. + a 

( 9.0 ) Si T /8 Ne ’ 
(5 = So = VJ/8 Avy 
of) 2 .. 


The significance of these relations lies in the fact that they state that the 
form of the fingers depends on the relation between A, and dg, i.e., on the 
relation between the critical wavelength and the wavelength of the hetero- 
geneities. In Section 4, the form of the fingers depended on the relation 
between A, and fo, i.e., between the critical wavelength and the initial form 
of the front, and in Section 3 it depended on the relation between the width 
of the system and the initial form of the front. 


6. STATISTICAL BEHAVIOR OF INSTABILITIES 

It is possible to treat the problem of instabilities in displacement processes 
in porous media statistically (see Scheidegger and Johnson 1961). Accordingly, 
one may again consider a two-dimensional displacement process (parallel to 
the x direction) but now only the average cross-sectional area occupied by 
the fingers is taken into account, the shape and size of the individual fingers 
being disregarded. The displacement of one fluid by the other is regarded as 
complete with no mixing occurring at the finger tips. The relative cross- 
sectional area occupied by the fingers is then analogous to a fictitious 
“saturation” s of the displacing fluid in the porous medium. It turns out that 
the equation of motion for this “‘saturation”’ is similar to the Buckley—Leverett 
equation for the saturation change in immiscible, stable displacement with 
a special expression for the relative permeability which becomes, in the 
present theory, a ‘‘fictitious’’ relative permeability. The system of the basic 


equations is then 


A » OS rz. OS a 
(6.1) I at r (S)qr a 0 
with 
(6.2a) r(s) = 

m—1-+(1/s) 
and hence 

9 Lite shmatncch ee aman 
(6.20) r’(s) [s@m—1) 41] 


where s is the ‘‘saturation”’ of the displacing fluid, q, is the filtration velocity 
of injection, P is the porosity, ¢ is time, and m is the viscosity ratio (of dis- 
placed versus displacing fluid). 

The problem is determined if the initial conditions are given, e.g. 


“s=0 fort <0, 
= hat =@ fort S10: 


(6.3) 4 
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Thus, the variables are s, x, and ¢; the parameters are P, m, and q;. There 
are two independent dimensions. 
We introduce the following dimensionless variables: 


(6.4) = = x/(wm), x = wmg; 
(6.5) i = tq,/(wP), t = wPi/q:; 
(6.6) § = s(m—-1), s = 5/(m-1); 


where w is an arbitrary length. Then, the equation of motion becomes 


Os 1 Os 
ait (§+1)° dt * 


(6.7) 
It is seen, thus, that one parameter (the viscosity ratio) remains. This is 
in keeping with the Pi-theorem inasmuch as of the three parameters (q,, m, P), 
two can be eliminated. 
The scaling relationships are 


(6.9) R(x) = R(w) R(m), 
(6.10) R(t) = R(w) R(P)/R(qz), 
(6.11) R(s) = 1/R(m-—1), 

and 

(6.12) R(m) = 1. 


The ratio of the lengths R(w) is arbitrary. 


REFERENCES 
Focken, C. M. 1953. Dimensional methods and their application (E. Arnold & Co., 
London). 
SCHEIDEGGER, A. E. 1960a. Phys. Fluids, 3, 94. 
19606. Geofis. pura e appl. 47, 41. 
SCHEIDEGGER, A. E. and JoHnson, E. F. 1961. Can. J. Phys. 39, 326. 








THE FRICTION OF HEATED BISMUTH! 


C. D. NIVEN 


ABSTRACT 

An apparatus is described to measure the friction of heated bismuth at loadings 
up to 640 p.s.i. A coefficient for pyrex glass on bismuth near its melting point of 
0.13 was indicated as a fair average value although figures down to 0.1 were 
also obtained. On surfaces where the metal was molten a value of 0.06 was 
measured. While bismuth has ice-like friction qualities up to a point, the author 
failed to prove that bismuth disobeys Amontons’ law at high loading in the 
manner that ice does. To account for this, the crystal structures of bismuth and 
of ice are compared. 

INTRODUCTION 

The idea of carrying out this investigation was prompted by studies on 
the friction of ice from which it was concluded (Niven 1959) that the peculiar 
friction qualities of frozen water originate from the unique structure of that 
material. This of course accounts also for the great expansion of water on 
freezing. Electrostatic forces compel the oxygen atoms of the water molecules 
to arrange themselves at the four corners of a tetrahedron with an oxygen 
atom in the center of these four atoms (Bernal and Fowler 1933; Owston 1958). 
In the liquid state this structure is essentially retained. The structure of 
bismuth (James 1921) consists of pairs of layers of atoms. The spacing between 
a pair of layers is much smaller than the spacing between the bottom layer 
of one pair and the top layer of the next lower pair. This layer structure of 
bismuth has been very clearly described diagrammatically by Krebs (1959). 
Each atom in the top layer of a pair is bonded to three atoms in the lower 
layer of that pair, while each atom:in the lower layer is bonded to three 
atoms in the upper layer. Bonding to the next pair of layers is much weaker. 
After bismuth has melted, however, a different situation arises, because as 
Krebs et al. (1959) have pointed out, the bismuth structure in the liquid state 
resembles to some extent that of rock salt and the distances between all 
layers of atoms tend to become nearly equal. The structures of bismuth and 
of ice differ quite markedly but both of these substances expand on freezing. 
Information, therefore, on the friction of bismuth near its melting point may 
conceivably contribute to a better understanding of ice friction. 

No work, to the author’s knowledge, has been done on the friction of heated 
bismuth at fairly high loading. The work which was done by Bowden and 
Freitag (1960) on bismuth was undertaken with a view to observing the effect 
of high speed rather than of high loading, but since ice does not obey Amontons’ 
law at high loading the friction phenomenon cannot be investigated thoroughly 
at temperatures near the melting point of one of the surfaces involved without 
providing for high loads per unit area of apparent contact. 


'\Manuscript received May 15, 1961. 
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METHODS 


(a) Apparatus 

In designing the apparatus the following two main requirements had to be 
taken into account: (1) Turntable temperatures had to be allowed to rise to 
271.3° C (m.p. bismuth) if the work were to be of any value whatever in 
comparing bismuth with ice because the interesting friction effects for ice 
are in the range of temperatures near and below its melting point. (2) In 
consequence of this requirement measurements had to be made in an inert 
gas atmosphere; this implied that all the experimental arrangements for 
loading the slider and for measuring the drag had to be done inside a working 
space which could be evacuated. 

Since cold bismuth does not oxidize quickly, the bismuth surface could be 
prepared in air by means of a removable lathe-head which could be clamped 
to the frame of the apparatus when required. Instead of the miniature sleigh 
used for the ice experiments a ring of the slider material was used. Arrange- 
ments were made to load this in vacuum which will be described below, and 
the drag was measured by determining the minimum force necessary to prevent 
the ring from being carried round with the surface of the turntable. 

How these operations were carried out from outside the gas-filled space 
will now be described with the help of Fig. 1. 

A motor M drove the main shaft S on which the turntable rested by means 
of the worm and gear g. A was the surface of the turntable; wires e led up 
through the shaft, which was hollow in its upper portion, to the heater h; 
these wires were connected to two rings r on which carbon brushes slid, and 
the power leads from a variac were connected to these brushes. The turntable 
was mounted on a cylinder C,;, which in turn was securely mounted on the 
end of the shaft S. This cylinder had to be reinforced later with an outer 
cylinder C,,; in order to increase the strength of the revolving unit. 

The shaft S had an end bearing at Bo and an ordinary ball bearing at By, 
mounted in a brass block P;, which could be water-cooled if desired. A Klosure 
Kx, provided a vacuun seal where the shaft entered the evacuated space. The 
slider F was made of the material whose friction on bismuth was to be deter- 
mined. In the tests this consisted of either pyrex glass or teflon, because since 
ice is non-metallic and bismuth is metallic, it seemed desirable to use a non- 
metallic slider if the results were to be compared with those obtained with 
stainless steel on ice. The slider was attached to the end of the cylinder C, 
by means of a screw with a very large flat head and could be removed and 
dipped in nitric acid to clean off the bismuth film which faintly silvered the 
polished glass after it had been running on the bismuth. The cylinder C, was 
screwed onto the loading rod L and clamped by means of a lock nut n. To 
permit the loading rod to turn through a small angle with negligible friction 
a teflon bushing t; was provided in the brass block Ps, which like P; could be 
water-cooled if desired. The upper end of the loading rod was drilled to permit 
the insertion of a small teflon bushing t2; a pin, which was secured firmly to 
the lid 1 of the weight box, dipped down into the bushing te. 
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Gas-tight turntable for measuring friction at high loading on a heated surface. 
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The small angle through which the sample F was permitted to turn was 
determined by the travel of a lever arm a in a notch cut in the blocking plate Z. 
This arm was firmly attached to the loading rod with a taper fitting and Allan 
set-screw biting into a flat on the loading rod. The extreme end of the arm a 
entered a notch in a right-angled plate Y which could swing around a pivot 
at q in such a manner that the horizontal movement of the arm a was con- 
verted into a vertical movement of an attachment x to a calibrated spring X. 
This spring was suspended from a rod passing through the Klosure seal K» 
out to a knob which could be pulled up by hand and clamped. This put a 
tension on the spring and pulled the arm over to one side of the notch in the 
blocking plate Z if the tension was greater than the friction force. If the 
tension were less than the friction force, the arm would be pulled over to the 
other side of the notch. When the arm was in the middle the tension of the 
spring, when multiplied by the ratio of the length of the arm a to the radius 
of F, gave a measure of the friction force. 

The weights were let down on the loading rod by means of the screw J 
which could be raised or lowered by the gears G. The lower end of J was not 
screwed where it passed through a Klosure seal K; in the lid of the weight box. 
The weights consisted of rings of lead Wi, We, Ws, W, resting on brass supports 
U, shaped like circular top hats. The total weights were designed to give loads 
on the loading rod of 10, 30, 60, and 100 pounds. The figures actually used 
in the calculations were 2 pounds heavier owing to the weight of L, Cz, and F. 
Three rods with suitably spaced stops i; on them could lower the weights one 
at a time onto stops iz on the loading rod by merely turning the gears G. The 
lifting rods R,; passed through the tops of the supports U through holes con- 
siderably larger than the rod diameter so that the weights which were down 
and loading F could turn through the small angle permitted by the blocking 
plate Z, but in order to keep the weights from working round the loading rod 
until they touched the rods R;, tapered pins were provided to guide the weight 
down so that it always sat in the correct position when down. The vacuum 
pump was connected at V. The small v’s at various places on Fig. 1 indicate 
evacuated space. The turntable was surrounded by a 4-in. diameter brass 
cylinder C; with large windows (not shown in Fig. 1) behind and in front of 
the turntable. This cylinder was clamped between the two O rings O in grooves 
in the blocks P; and P. which were drawn together by means of the nuts 
N on the ends of the rods Ro. 


(b) Measurements 

Before making a set of measurements the following routine had to be 
followed: (1) the weight box had to be lifted high by the screw to permit 
clamping on the lathe-head and machining of the bismuth; (2) the slider 
had to be cleaned in acid; (3) the cylinder C; had to be placed on the O ring, 
the weight box lowered, the vacuum pump started, the argon filled in, and 
the heat turned on under the turntable; (4) after a definite time the motor 
was started with negligible load—it was very important never to start the 
turntable when the load was on the slider and, after it had once been started, 
never to stop it until the load was taken off. 
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The method of applying the weights permitted very rapid changes in loading 
because a weight would come down by simply turning the gear on the screw 
about 23 revolutions. One could tell when a weight had come down by merely 
listening for a change in the sound of the running. For each load the tension 
of the calibrated spring was noted when the connecting rod was in the middle 
of the hole in the lid. Tension values in grams were transformed by means 
of a graph to friction force in pounds between the slider and the bismuth 
surface. The transformed figures or drag values could then be plotted against 
the load values, which were in pounds (see Figs. 2, 3, 4). 
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Fic. 2. Drag values at various loads for }-in. pyrex glass ring on solid bismuth heated 


above 240° C. 


Figure 2 refers to results using a ring of pyrex glass 23 in. in diameter and 
§ in. wide. Some of the points plotted in Fig. 2 were repeated, but no attempt 
has been made to “weight” these points. The two lines which form an angle 
containing most of the points act roughly as limits for the range of coefficients. 
A line bisecting this angle intersects the 100-pound load ordinate at 12.5; 
the results therefore suggest 0.13 as a fair value to take for the coefficient of 
friction. 

These two lines have been drawn also on Figs. 3 and 4 to allow the reader 
to see at a glance how the results plotted there compare with the results in 


Fig. 2. Figure 3 refers to a narrow castellated glass ring, that is to say, a glass 
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bic. 8. Drag values at various loads for castellated pyrex glass on solid bismuth heated 
above 240° C 
hic. 4. Drag values at various loads for teflon ring on solid bismuth heated above 240° C. 
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ring which was ,, in. instead of § in. wide and which had been ground out to 
leave three little runners each about ? in. long. In this way a pressure of 640 
p.s.i. could be applied at the interface when the 100-Ilb load was used. Figure 
4 refers to a ring of teflon about 22 in. in diameter and originally 3 in. wide. 
This width increased considerably as it wore down. In Figs. 3 and 4 it has 
been possible to show the results of different runs by different symbols. 

Figures 2, 3, and 4 all refer to the same approximate temperature on the 
bismuth surface. This temperature had to be estimated from the length of 
time current passed through the heater because the possibility of getting a 
thermocouple poked in on the moving bismuth surface was out of the question. 
At the start of the work an attempt was made to ascertain the temperature 
by pressing a thermocouple against the side of the turntable when it was not 
turning but the temperatures measured in this manner were obviously too 
low and were not really of any use in finding out the temperature of the 
bismuth surface. This was proved when the bismuth melted. Accordingly the 
thermocouple was placed under the slider and actually on the bismuth with 
the table not revolving. Temperature was then read as the surface heated 
up at a particular setting on the variac delivering the current. By plotting 
thermocouple reading against time when the table was still, a rough estimate 
could be made of the temperature when the table was turning or prepared for 
turning. The method is of course only approximate, firstly owing to changes 
in line voltage, secondly owing to the difference in thermal contact between 
the slider and the bismuth when the couple was in and when it was out, and 
thirdly because the method takes no account of the frictional heating due to 
prolonged running. The speed of the slider over the bismuth was 1.4 f.p.s., 
so when the friction force was 6 lb, the work done was 8.4 ft-lb per sec or 
12.6 watts was being supplied. To minimize this last error the readings were 
taken as quickly as possible and then the motor was stopped if another set of 
readings had to be taken. In the later stages of the work the table was allowed 
to heat up for about 70 to 75 minutes at a particular variac setting. One or 
two sets of readings were taken as quickly as possible and then the apparatus 
was opened immediately after to examine the condition of the bismuth track. 
This was of special interest when the castellated pyrex sample had been running 
and exerting a pressure of about 640 p.s.i. on the bismuth surface at the 
interface. Bismuth ordinarily is regarded as quite a brittle metal, but the 
warming of the track by means of the heater coupled with the frictional heat 
and the pressure resulted in a deep grooving of the track and pushing up of 
the metallic bismuth on each side of the track to a height of 4 in. or jt in. 
above the original level. The original level incidentally could easily be seen 
farther out from the track because the tool marks were still quite visible on it. 
The bottom of the groove was perfectly smooth and the running under heavy 
load on these three small glass runners about $ in. long and }. in. wide was very 
quiet and steady. The plastic flow of the bismuth to the side of the track agrees 
with an observation of Bowden and Freitag (1960) that plastic deformation 
could occur in bismuth near its melting point. The pyrex ring with the } in. 
bearing surface polished the track but caused negligible grooving. 
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In Fig. 3, three drag values for the 102-lb load are very close together, 
around 9 or 10 lb. These were obtained on three entirely different consecutive 
runs. Such consistently low values for the drag at high loading suggested at 
first that there might be a curving round of the drag vs. load curve towards 
the load axis at high loads for bismuth much as had been found for ice even 
at as low a temperature as —17°C (Niven 1954). However, after plotting 
the results the writer decided that there was really insufficient evidence to 
substantiate this, partly because some quite low drag values were recorded at 
low values of load and partly because on two occasions only a few minutes 
earlier higher values had been obtained for the 102-lb load. These are plotted 
with the same symbol. 

The drag values for teflon shown in Fig. 4 compare fairly well with those 
for pyrex glass shown in Fig. 2 suggesting that the non-adhesive quality of 
teflon is not an important factor near the melting point of bismuth. At much 
lower temperatures the friction of the pyrex glass was much higher than that 
of the teflon. Indeed at 185° C or near this temperature, the pyrex glass would 
not run smoothly and at high loads the friction was so great that the apparatus 
was endangered. Teflon on the other hand appeared to have a coefficient of 
friction almost independent of temperature in that temperature range. This 
suggests that at the higher temperatures the strength of the asperity is the 
important item, while at lower temperatures when the asperity must be 
strong and its structure very rigid the adhesive bond at the interface must be 
the important item. 

At the start of the work some results were obtained on bismuth surfaces 
which had obviously melted by frictional heating as a result of prolonged 
running. On one such occasion friction coefficients as low as 0.06 to 0.07 were 
obtained; such working conditions were, however, strictly avoided in the 
later stages of the work because they do not represent the correct condition 
for dry friction between pyrex glass and solid bismuth. The value given by 
Bowden and Freitag (1960) for the friction of bismuth when surface melting 
occurs is 0.04; the agreement is not bad considering that Bowden and Freitag 
were using speeds as high as 200 m/s which would favor surface melting. The 
value 0.04 can be taken as a lower limit for the friction of bismuth as the surface 
gets nearer and nearer to the melting point. 

DISCUSSION OF RESULTS 

There is quite a large spread in the drag values. This is believed to be due 
to the nature of dry friction measurements. The indefiniteness of the tempera- 
ture does not account for this spread because drag values would sometimes 
change suddenly during a test. The indefiniteness of the temperature discussed 
above refers more to whether the measurements were made near 240° C or 
270° C. They were probably made in a small range of temperature quite 
near 270° C because there was always a little anxiety lest the bismuth would 
melt and thereby cause trouble in preparing for the next run, but because an 
accurate method of temperature measurement was not worked out the results 


in Figs. 2, 3, and 4 may only be specified as applying to bismuth above 240° C. 








1272 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


The heat should not have been put inside the turntable because it heated 
the bottom of the bismuth first. To apply the heat from above and to control 
it to within a few degrees of 271°C obviously involves major changes in 
design, and when the apparatus was built the main problem was to vary the 
load, to apply very considerable loads, and to measure the force in a gas-tight 
space from outside. 

The coefficient of 0.13 is lower than the coefficient for stainless steel on ice 
at —17° C using a large apparent area of contact; this was around 0.2 (Niven 
1954). The very low coefficient for hot bismuth, as well as the smooth running 
of the table at heavy loading and the highly polished track with the plastic 
flow of metal at the edge, point to a great similarity between the friction of 
ice well below 0° C and of hot bismuth near its melting point. However, the 
coefficient of 0.04 on bismuth with surface melting is higher than the co- 
efficient of ice around 0° C; at —2° C the writer actually measured a coefficient 
of 0.003 on ice using high loading (Niven 1954). The castellated pyrex on 
bismuth did not show any coefficient nearly as low as this and furthermore 
there is, as pointed out above, insufficient justification to claim that Amontons’ 
law was violated at high loading as it is in the case of ice. Accordingly bismuth 
has ice-like qualities only up to a point. 

Such reservation fits in satisfactorily with the descriptions of the structures 
of ice and of bismuth outlined in the Introduction. Since each atom in the top 
layer of a pair of layers which are close together is bonded to three atoms in 
the lower layer of that pair, there exists a large plane substructure which cer- 
tainly cannot be termed open. Therefore, although the bismuth structure 
taken as a whole is open, because the pairs of layers are relatively far apart, 
yet one cannot expect these large substructures to swivel easily in the manner 
in which the open tetrahedral ice structure might be expected to swivel. 
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THE ABSORPTIVITY SPECTRUM OF A UNIFORM, 
ANISOTROPIC PLASMA SLAB! 


I. P. FrENcH, G. G. CLOUTIER, AND M. P. BACHYNSKI 


ABSTRACT 


Using a generalized formulation of Kirchhoff’s law it is possible to relate the 
equilibrium electromagnetic radiation spectrum of a body to its absorptivity 
spectrum. The microwave absorptivity of a uniform anisotropic plasma slab to a 
normally incident electromagnetic wave is obtained by matching fields at the 
boundaries and is solved completely for the model chosen. 

The absorptivity spectrum of the plasma slab is computed for waves propa- 
gating parallel to and normal to an applied static magnetic field for various 
electron densities, electron collision frequencies, and slab thicknesses. Peaks in 
the absorptivity spectrum occur around the edges of the ‘‘stop-bands’’. In 
general, the absorptivity increases with slab thickness and collision frequency. 
The effect of internal reflections is included and gives rise to undulations in the 
absorptivity spectrum. 


I. INTRODUCTION 

The passive electromagnetic radiation from various types of plasmas is of 
great current interest in astrophysics, thermonuclear fusion, communications, 
and laboratory discharges. Of particular interest here are fairly low tempera- 
ture (less than 10,000° K) plasmas such as the plasma sheath of a re-entry 
vehicle. The passive radiation in the radio-frequency range of the spectrum 
will alter considerably the effective noise temperature of an antenna located 
in such a vehicle (Bachynski, French, and Cloutier 1961) and may also lead 
to a missile detection system. 

There are two main possible approaches to the problem of radiation from 
plasmas. The first considers microscopic phenomena such as bremsstrahlung 
and cyclotron radiation from charged particles in a magnetic field and is 
particularly convenient when the plasma is so tenuous that the effects of 
self-absorption can be neglected (Rosner 1958; Hirshfield 1959). The second, 
more general approach, which is used in this paper, derives the radiation from 
the macroscopic electromagnetic properties of the plasma by means of a 
general formulation of Kirchhoff’s law. 

The absorptivity of a uniform anisotropic plasma slab to an electromagnetic 
wave incident normally upon it has been computed for various electron 
collision frequencies and slab thicknesses. Propagation parallel to and per- 
pendicular to a high-intensity static magnetic field is considered and the 
various possible electromagnetic ‘‘pass’’ and ‘‘stop-bands”’ for the different 
wave types are discussed. One of the main characteristics of the spectrum is 
the enhanced absorptivity which occurs around the edges of the ‘‘stop-bands”’. 
The internal reflections of the electromagnetic energy in the slab give rise to 
undulations in the absorptivity and radiation spectra. 

1Manuscript received April 13, 1961. 
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II. THE ABSORPTIVITY OF A PLASMA SLAB 
Rytov (1953) and Levin (1955) have shown that the emissive power of 
a body at a frequency w may be expressed by 


(1) Pi = =—s A, 
T 


where A, is the power absorptivity of the body (an absorptivity of unity cor- 
responds to a black body), 

K is Boltzmann’s constant, 

T is the temperature of the body, 

d is the free-space wavelength, 

P., is the radiated power per unit frequency interval, per unit solid 
angle, per unit surface area for a given polarization. 
This general expression is valid for any degree of absorptivity so that the 
equilibrium radiation spectrum of a body can then be predicted from the 
spectrum of the absorptivity. The absorptivity, 4, is defined as the fraction 
of energy absorbed by a body when a plane-polarized electromagnetic wave 
is incident upon it, and is given by 
(2) 4, = ArbReG-E*)do_ 

| .4Re(E; XH ;*)ds 
where J and E are the current density and the electric field in the body due 
to an incident electromagnetic field (E; and H,) on the body. The numerator, 
which represents the energy absorbed by the body and dissipated as joule 
losses, can, in principle, be deduced from the electromagnetic properties of the 
body. The denominator is the energy incident on the body. The evaluation 
of A, from eq. 2 is a boundary-valué problem which has only been solved 
for very simple geometries (Rytov 1953; Levin 1955). Such solutions are even 
more complicated for anisotropic and inhomogeneous bodies. 
In terms of electromagnetic parameters the absorptivity may be expressed 


(3) A, = 1-R,-T, 


where RX, and 7, are the power reflectivity and power transmissivity respectively 
of a plane wave incident normally on the surface of a body. R, and 7\, are 
evaluated by matching at the boundaries the incident and transmitted fields 
with fields within the body. 

Propagation in an unbounded anisotropic plasma in any direction can, using 
Maxwell’s equations, be described in terms of two waves whose properties 
are a function of the orientation and magnitude of the static magnetic field. 
For mathematical simplicity orientations of magnetic field parallel to and 
perpendicular to the direction of propagation are usually chosen. The so-called 
electron and ion cyclotron waves propagate parallel to the static magnetic 
field and the ordinary and extraordinary waves correspond to waves propa- 
gating perpendicular to the static magnetic field (see, for example, Bachynski, 
Shkarofsky, and Johnston 1961). The ordinary wave is identical with that in 


an isotropic plasma. 
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Using a simple model, the absorptivity spectrum for a plasma in equilibrium 
in a constant d-c. magnetic field is derived. The plasma geometry of a slab 
whose dimensions are infinite along the x—y axis (Fig. 1) and of thickness d 
along the z-axis is chosen. This configuration of plasma permits diffraction 
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Fic. 1. The co-ordinate system. 


effects, which are important when the plasma dimensions are not much greater 
than the wavelength of the radiation, to be neglected. It is assumed that the 
plasma is in thermal equilibrium and the electron density is uniform throughout 
the plasma, implying infinitely sharp boundaries. The effect of internal re- 
flections is included. 

Let a plane electromagnetic wave (£,ye/*') be normally incident on the 
plasma slab with its E-vector along the x-axis. The fields in region 1 are 


Fi, = Ez ele FLT e”*], Ey, = Ei, = 0, 
(4) 


Eo iw 02 . 02 
Hy = Fee" —Toe"}, Hi = Hy = 0. 


The fields in region 2 are 


Bee = Eye" le ve r,e"'), Eo, = Es, = 0, 
(5) : 

H2, = z ele? — Te], Ho, = Hy, = 0. 
The fields in region 3 are 

E3 = T. Bie é 7 Ex, = E3. = 0, 

) che Eo iw 2 
(9) Hy = To Zte™e™, Hig 0 By wit 
si) 


where yo and y are the propagation coefficients in free-space and the plasma, 


respectively, with 
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(7) y = jkK'? = at+jp 


where k = 27/)X, 

K = K,—jK, is the dielectric coefficient of the plasma, 

a and @ are the attenuation and phase coefficients of the electromagnetic 
wave in the plasina, 

I, and [, are the amplitude reflection coefficients for the waves incident 
at boundaries 1—2 and 2-3, respectively, 

To is the amplitude transmission coefficient for the waves incident on 
boundary 1 and transmitted through boundary 3, 

Z, and Z are the characteristic impedance of free-space and the plasma, 
respectively, with 

1 Zo 


(8) Ze ee 
V (equcK ) VK 


€9 is the permittivity of free-space. 
The permeability of the plasma is taken equal to the permeability of free- 
space, mo. 

The real component (K,) and the imaginary or loss component (K,) for 
the various waves are listed in Table I. lon effects which give rise to very 
low frequency resonances are not included. 

Matching the fields at boundaries 1-2 and 2-3 gives 


Zee o 
( T ~~vod a 20 . 
(9) To ¢ | cosh yd+s5 eae Z ) sinh vd | : 
es > sinh yd 
(10) Pik coe 2\Z Z ; — 


cosh yd+5 2 2), sinh yd 
0 


By making the substitutions 


P= B (ke ‘ta‘+6°) fa = (Ri =a" 6 ) 
(11) ee ee he ee 
Ta 8 (k'—a'—8') ss (Re tal +6" +) 
2 ate Ra te 
then 
Z Lo 
ee Same oem . 
ZZ 2(L+jM), 
(12) 72 
Fe EEN ey 7) 
ZZ 2(P+jQ). 


The power transmissivity 7, and reflectivity R, are given by the products 
of their amplitude coefficients and their complex conjugates. This leads to 
(from eq. 3): 
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(13) A, =1- 
1+(L?+M")(sinh’ad+sin" 6d) 





r+ Q*) (sinh’ad+sin*Bd) +2P cosh ad sinh ad+sinhad+cos*8d - 
— 2Q sin Bd cos Bd 


This relation gives the absorptivity of a uniform plasma slab for an incident 
plane electromagnetic wave in terms of the electromagnetic properties and 
thickness of the plasma. 

Simple approximations (Table II) are possible for both high and low radio 
frequencies relative to the plasma frequency wp, = ne?/meo and the electron 
cyclotron frequency w, = eBo/m, where n is the electron density, e and m are 
the electronic charge and mass, respectively, and By is the static magnetic 
field in m.k.s. units. 

At frequencies much larger than both the plasma frequency (w,) and the 
electron cyclotron frequency (w)) the plasma is essentially isotropic and the 
absorptivity is independent of the orientation of the static magnetic field. In 
this portion of the spectrum the plasma appears tenuous to the electromagnetic 
wave. The absorptivity is then directly proportional to the collision frequency 
and the slab thickness and varies inversely as the radio frequency squared. 
This implies (from eq. 1) that the radiation spectrum is independent of fre- 
quency in the regions where self-absorption is negligible, as predicted by 
Dellis (1958). Bekefi, Hirshfield, and Brown (1959) also arrive at the same 
conclusion. In this part of the spectrum the free-space— plasma interface has 
a very small reflection coefficient, implying a good match between free-space 
and the plasma, and a reflectivity close to zero. 

On the other hand, at frequencies close to and below the plasma and 
electron cyclotron frequencies ‘‘stop-bands” are possible for the various waves 
(Table III). In a stop-band the reflection coefficient of the free-space — plasma 
boundary is large, resulting in a bad mismatch and low absorptivity since 
little energy can penetrate the boundary. 

At very low radio frequencies relative to the plasma and electron cyclotron 
frequencies the plasma remains anisotropic and the absorptivity is strongly 
dependent on the magnitude and orientation of the magnetic field. The 
approximations are listed in Table II. 

An effect which is likely to modify the absorptivity spectrum is that due 
to multiple internal reflections of the energy between the boundaries of the 
slab. For a continuous change in slab thickness, d, or plasma phase coefficient, 
8, the functions sin Bd and cos 6d will go through maximum and minimum 
values every time 8d passes through an integral number of 2/2 in eq. 13. 

If the approximation (8/k)? > (a/k)? is made with m an integer, 


_ 2tanh ad(P+tanh ad ) 


A, = (1-+-P tanh ad)? when Bd = nz, 


(14) 


2 tanh ad 


Ay P+tanad 


when 6d = (n+34)z. 
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TABLE III 


Stop-bands in a plasma 











Type of wave Lower frequency limit Upper frequency limit 
Electron cyclotron Wb V (Jon)? +03 +4or 
Ion cyclotron 0 V (Jon)? +08 — Jun 
Ordinary 0 Wp 
Extraordinary V ot+w3 V (Jon)? +02 +4or 

(two stop-regions occur) 0 v( Jun)? +w2 — wp 





From eq. 11, P > 1 and therefore Bd = nz gives a local maximum in absorp- 
tivity and Bd = (n+4)z a local minimum. Thus it is to be expected that 
undulations in the absorptivity spectrum will occur since the phase coefficient, 
8, is a function of the radio frequency. 


Ill. THE ABSORPTIVITY SPECTRUM 

The absorptivity spectrum has been computed (using eq. 13) as a function 
of normalized frequency (w/w,) for a uniform static magnetic field both along 
the direction of propagation and transverse to the direction of propagation. 
Four spectra result, corresponding to the four possible waves discussed earlier; 
the dielectric coefficient for each wave is listed in Table I. The effect of electron 
density and collision frequency on the absorptivity is investigated by com- 
puting a set of curves for a fixed plasma thickness (d/A, = 2.5) and values 
of v/w, ranging from 10-* to 1 (Figs. 2-5). The variation of absorptivity 
with plasma thickness is for a value of v/w, = 10-' and values of d/\, ranging 
from 0.05 to © (Figs. 6-9). A value of w/w) = 2 was chosen. (A, is defined 
as the wavelength corresponding to the plasma frequency and is hence a 
measure of electron density in the plasma.) 

In any practical case the total absorptivity of an anisotropic plasma slab 
in a direction normal to its boundaries is a combination of the absorptivities 
of the two possible waves in this direction. For example, if the externally 
applied magnetic field By is along the direction of propagation the absorptivities 
of the electron (-1_) and ion (.1,) cyclotron waves must be combined. 

The linearly polarized uniform plane wave carrying unit power considered 
in the first section can be regarded as two equal-amplitude circularly polarized 
waves rotating in opposite directions. One of these circularly polarized waves 
is the electron cyclotron wave and the other the ion cyclotron wave. Since 
each wave carries half the power the absorptivity is the mean of A_ and 
Ace: 


(15) Ay = 3(A_+A4). 
Similarly, when the applied magnetic field is perpendicular to the direction 


of propagation the total absorptivity is the mean of the absorptivity of the 
ordinary wave and extraordinary waves if a randomly polarized incident wave 


is used. 
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Fic. 2. Absorptivity (A) of an anisotropic plasma slab as a function of radio frequency 
for the electron cyclotron wave (d/Ap = 2.5, wp/wp = 2) for various electron collision frequencies. 

Fic. 3. Absorptivity (Aw) of an anisotropic "plasma slab as a function of radio frequency 
for the ion cyclotron wave (d/Ap = 2.5, w»/wp = 2) for various electron collision frequencies. 
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Fic. 4. Absorptivity (Aj) of a plasma slab as a function of radio frequency for the ordinary 
wave (d/A, = 2.5) for various electron collision frequencies. 
_ Fic. 5. Absorptivity (Ay) of an anisotropic plasma slab as a function of radio frequency 
for the extraordinary wave (d/Ap = 2.5, w,/wp = 2) for various electron collision frequencies. 
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Fic. 6. Absorptivity (A,) of an anisotropic plasma slab as a function of radio frequency 
for the electron cyclotron wave (v/wp = 107', w)/wp = 2) for various slab thicknesses. 

Fic. 7. Absorptivity (Aq) of an anisotropic plasma slab as a function of radio frequency 
for the ion cyclotron wave (v/w,p = 107, w,/wp = 2) for various slab thicknesses. 
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Fic. 8. Absorptivity (A) of a plasma slab as a function of radio frequency for the ordinary 
wave (v/w, = 107) for various slab thicknesses. 
_ Fic. 9. Absorptivity (A,) of an anisotropic plasma slab as a function of radio frequency 
for the extraordinary wave (v/wp = 107, w,/wp = 2) for various slab thicknesses. 
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(a) The Effect of Stop-bands 

For a lossless plasma (zero collision frequency, v) the stop-bands are regions 
where the normalized phase coefficient (8/%) becomes zero and/or the attenua- 
tion coefficient @ is infinite implying that no propagation can occur (Table 
III). The term “stop-band” for a plasma which has a small but finite number 
of collisions (v Kw,) implies a region where the normalized phase constant 
(8/k) is very much smaller than unity. It has little meaning when the collision 
frequency is high (v ~ w)). Generally there is a maximum in absorptivity 
just outside a stop-band, since (8/k) is not very different from unity and 
most of the energy penetrates the first boundary and is absorbed due to the 
high attenuation coefficient, a. 

In the stop-bands the plasma is badly matched to free-space (8/k not close 
to unity) resulting in a large reflection from the first boundary and low 
absorptivity since little energy enters the plasma. 

Outside the stop-band the absorptivity is in general greater since the match 
to free-space is much better and most of the energy can enter the plasma, 
although the attenuation coefficient, a, is much smaller. 


(b) The Effect of Electron Collision Frequency 

In general, a high collision frequency implies a very lossy plasma and high 
absorptivity over the whole spectrum. A low collision frequency implies a 
bad mismatch to free-space in the stop-bands and therefore very low absorp- 
tivity in this region. Over the rest of the spectrum the absorption path length, 
ad, is small and decreases as the electron collision frequency decreases, implying 
low absorption. 

For example, at the center of the stop-band, 6/k = 0.06 for the electron 
cyclotron wave when v/w, = 10 and d/d, = 2.5. Decreasing the collision 
frequency by a factor of 10 to v/w, = 10-* reduces 8/k to 0.006 implying a 
much higher reflection. 

The power reflection coefficient from the first boundary, or from a semi- 
infinite plasma slab, is given by (Bachynski, Johnston, and Shkarofsky 1960) 

(1—8/k)*+(a/k)° 


(1+8/k)*+(a/k)” 
The energy which does penetrate the first boundary is, however, all absorbed 
in its first passage through the slab since ad is very large (~100—200) in the 
examples considered above. 


(c) The Effect of Slab Thickness 

In general, the absorptivity increases with slab thickness (Fig. 10) except 
in the stop-regions when the attenuation path length, ad, is large and all 
the energy which penetrates the first boundary is absorbed (Fig. 11). In the 
pass-bands undulations due to internal reflections are to be expected when 
the thickness is varied (Fig. 10). 

A further characteristic is that the absorptivity of a very thin slab is often 
greater than that of an infinite slab at very low frequencies (Fig. 12). This 
effect, which cannot be predicted from a geometrical optics ray tracing 
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Fic. 10. Absorptivity (A,) as a function of slab thickness for the ordinary wave for a 
frequency w = 1.2 wp. 
Fic. 11. Absorptivity (Aj) as a function of slab thickness for the ordinary wave for a 
frequency w = 1.0 wp». 
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Fic. 12. Absorptivity (A,) as a function of slab thickness for the ordinary wave for a 
frequency w = .1 wp». 


approach, is a result of the exact solution of this boundary value problem. 
Thus the reflectivity of a slab which is thinner than the penetration depth of 
the electromagnetic wave is very low in the stop-bands. Consequently, such 
a thin slab absorbs and transmits a large portion of the energy incident on 
its surface. The penetration depth, ds, is given by the reciprocal of the attenua- 
tion coefficient (1/a), and for frequencies between the plasma frequency (w,) 
and the electron collision frequency (v) is ds = C/w, for the ordinary wave or 
for an isotropic plasma. This corresponds to a value of d/A, = 0.16. 

For thin slabs (0.05 < d/\,) the characteristic maxima in absorptivity at 
the boundaries of the stop-bands disappear (Figs. 6 and 9). The absorptivity 
spectrum for the electron cyclotron wave becomes perfectly symmetrical with 
a peak at the cyclotron frequency, w,. The frequency interval between half- 
power points on the ener spectrum when v/a, = 107 is equal to 4.5y 
for d/A,y = 0.05 and 2.3v for d/A, = 0.005. The reflectivity of the plasma to 
this wave has a maximum of 0.02 at w, when d/A, = 0.005 and consequently 
boundary effects are of little importance. The ‘‘pressure-broadened” (due to 
collisions only) line width at half-power points in a radiation spectrum is 
expected to be twice the collision frequency (2v) when boundary effects are 
neglected (Margenau and Lewis 1959). In general, boundary reflections may 
be neglected from very thin slabs (0.005 < d/A,). 


(d) Undulations in the Spectrum 
The effect of sharp plasma boundaries leads to internal reflections in the 
slab and undulations in the absorptivity and radiation spectra. Peaks in 
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the undulations occur when 6d = nz (man integer) and dips when 6d = (n+4)zr 
as shown in Fig. 13 for the radiation spectrum of the ordinary wave or for 
an isotropic plasma. The curve 6d/z is shown for comparison. Curve 1 shows 
the radiation spectrum when d/A, = 2.5 and v/w, = 10-*. Curve 2 shows the 
effect of increasing v/w, to 10—' and decreasing d/), to 0.25. The characteristic 
peak around the plasma frequency is not so pronounced in the second case. 
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Fic. 13. Radiation spectrum (arbitrary units) of an isotropic slab: curve (1) v/wp = 107%, 
d/» = 2.5, Bd/x from 0 to 20; curve (2) v/wp = 107, d/Ap = .25, Bd/x from 0 to 2. 


In the regions where undulations are apparent 8/k is not close to unity 
and therefore the plasma is badly matched to free-space. The absorption path 
length, ad, is small (<1), otherwise the wave would be highly attenuated in 
its first passage through the plasma. In addition, the thickness, d, should be 
large enough that Bd becomes greater than 7/2 for undulations to occur. 


(e) The Effect of Electron Density and Magnetic Field 
The effect of changing electron density is to change the plasma frequency 
in accordance with the relation 
») > 
2 216 
(16) w=— and »»W=—— 
€om Wp 
The absorptivity curves (Figs. 2 to 9) are normalized with respect to plasma 
frequency, wp, or free-space wavelength corresponding to the plasma frequency, 


\p, and are shown in terms of the parameters v/wy, w/wy, and d/d, for a value 


of w/w, = 2. 
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In general, changing the magnetic field changes the position of the stop- 
bands and the regions in the spectrum of enhanced radiation. For example, 
if the electron cyclotron frequency (w,) is much larger than the plasma fre- 
quency (w,) then the stop-band for the electron cyclotron wave and the 
upper stop-band for the extraordinary wave become extremely narrow and 
their positions almost coincide with the electron cyclotron frequency. Conse- 
quently, only a single peak about w, would appear in the absorptivity spec- 
trum. In such a case the effect of the ordinary and ion cyclotron waves on the 
total spectrum would be small. 

Alternatively, making the plasma frequency much larger than the electron 
cyclotron frequency stretches out the stop-bands considerably. The stop-band 
for all the waves would extend effectively from zero to w,, with a narrow 
pass-band for the electron cyclotron wave below w,. Consequently, enhanced 
absorptivity should be observed mainly around w, with a small peak around 
w, (if » < w)). Consequently, unless the electron cyclotron frequency is of 
the same order of magnitude as the plasma frequency the frequencies of 
enhanced absorptivity and radiation will be around o, and ap». 

In most plasmas of interest the electron density at the slab boundary will 
not suffer the sharp discontinuity assumed in this analysis. Diffuse boundaries 
are expected to lower reflections and lead to higher absorptivities, especially 
in the stop-bands. Reflectivities lower than those expected from a slab analysis 
have been observed experimentally by Bekefi, Hirshfield, and Brown (1959) in 
discharges in waveguides and are attributed to gradual boundaries. 


IV. CONCLUSIONS 

In general the absorptivity spectrum of a plasma is very low in the stop- 
bands since little energy can penetrate the first boundary. In the pass-bands 
the absorptivity depends mainly on the absorption path length, ad. The 
spectrum, in general, has a maximum around the edges of the stop-bands 
since in this region the plasma is fairly well matched to free-space and the 
attenuation coefficient is large. The effect of a high electron collision frequency 
is to broaden the peaks. The absorptivity of a very thin slab may be high in 
the stop-bands due to a low reflectivity. The effect of internal reflections is 
to produce undulations in the spectrum. 

Since it is possible to relate the absorptivity spectrum of a body to its 
equilibrium radiation spectrum the solution presented here predicts the effect 
of such fundamental parameters as electron concentration, electron collision 
frequency, and magnetic field on the radiation spectrum of the plasma. 
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ON BOURRET’S HYPOTHESIS CONCERNING 
TURBULENT DIFFUSION! 


P. H. ROBERTS? 


ABSTRACT 


An integrodifferential equation proposed on heuristic grounds by Bourret 
for turbulent diffusion is compared with a similar equation which was derived 
by Roberts from the application of Kraichnan’s approximation to the exact 
equations of turbulent diffusion, and which is itself closely related to an equation 
postulated by Bourret as a result of further heuristic arguments. A simple 
generalization of Taylor’s model of turbulent diffusion is discussed, and is 
found to obey Roberts’ equation but, unless the spacial correlation between 
the velocity of the diffusing particles is zero, it does not obey the equation 
initially proposed by Bourret. 


1. INTRODUCTION 


According to the theory of molecular motion in a dilute gas of discrete 
particles, the probability density P(x,t) of particular particles marked in 
some way is governed by an ordinary diffusion equation of the type 


9 


(1) < P(x,t) = «iy ae P(t), 

provided we consider distances x large compared with the mean free path and 
times ¢ large compared with the mean collision time. Over these large distances 
and times, molecular diffusion may be thought of as a random walk process 
in which the effective step length is the mean free path. Likewise, if we consider 
distances large compared with /, the turbulent macroscale, and times large 
compared with //vo, where vp is the r.m.s. turbulent velocity, turbulent diffusion 
may be thought of as a random walk process in which the effective step length 
is / and the effective velocity is v». We may expect, therefore, that over these 
large times and distances, turbulent diffusion in a statistically homogeneous 
and stationary flow will be governed by (1), and that the constants «;; should 
depend only upon / and vp, i.e., by dimensional reasoning, be of the order of 
luo. The correct form for x;; was found by Taylor (1922) to be 


(2) Ki = j M,,(t’)dt’, 
7 

where 

(3) M(t’) = (o,(t)v,(t+t’)) 


is the Lagrangian velocity correlation, v(t) being the velocity of a marked 
fluid particle at time ¢ and v(t+?’) being the velocity of the same particle at 
an instant ¢’ later, the average (...) being over ensembles. 
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Over distances small compared with / and times short compared with //vo, 
we would not expect (1) to be obeyed since, over such distances and times, 
the effects of persistence of velocity correlations over finite space and time 
intervals should play an essential part. In place of (1) we might expect an 
integrodifferential equation in which an average is taken (with appropriate 
weighting) over the past history of the particles’ motions. Such, in fact, is an 
equation proposed by Bourret (1959, 1960): 

(4) . P(x,t) = = fate) Pasa’ 

gO!” Andee ea 
Bourret derived his result by generalizing a property of a model due to Taylor 
(1922) which includes the effects of velocity persistence. We shall presently 
(Section 2) give a very brief resume of the steps which led Bourret to (4), but 
before doing so we note a second integrodifferential equation which has been 
proposed by Roberts (1960, 1961) as a replacement of (1): 


0 a s j 
5 eee a er ’ 11). (w—y’ 4—#'\) Ply’ #! 
(5) > P (x,t) mae, J dt’ | dx’Q;,(x—x’, t—t’) P(x’,t’). 


(The integration over x’ is over all space.) Bourret (1961) has also derived an 
equation similar to (5) by heuristic arguments. In order that (2) be a conse- 
quence of (5), the functions Q;,(x’,t’) appearing in (5) must be related to 
M(t’) by 


(6) Mult) = fax Qu(x'e’. 


However, Roberts found that, on using the approximation scheme of Kraichnan 
(1959), 


(7) Q:,(x’,t’) = P(x’,t’) Uij(x't’), 
where 
(8) Ui(x’,t') = (u(x,t) u(x+x’, t+0’)) 


is the Eulerian velocity correlation between the motions of (different) fluid 
particles at the positions and times indicated. The approximation (7) is not 
consistent with (6), although the agreement may be expected to be good for 
large t’ (see Roberts (1961) for a fuller discussion). In what follows we do 
not suppose that Q;,(x’,t’) is given by (7) but we shall interpret Q;,,(x’,t’) in 
terms of Green’s functions for the equation governing the diffusive flux (see 
Section 3 below). 

In comparing (4) and (5) we notice at once that (4) is non-local in time but, 
unlike (5), is local in space. This implies that (4) accounts for persistence of 
velocity correlation along the path of a particle (i.e. a finite Lagrangian cor- 
relation time), but it does not seem to account for the fact that a cloud of 
marked particles set down at random in a turbulent flow must also display 
velocity correlation in space at any instant of time. It would appear that 
such a correlation effect is included in (5) through the spacial non-localness. 
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In Section 4, a simple ‘‘two-particle’’ model is studied in which an equation 
of the form (5) is found to be obeyed if there is spacial coherence between the 
two particles concerned, an equation of the form (4) being obeyed if and only 
if there is no spacial correlation between them. 


2. BOURRET’S HYPOTHESIS 

Bourret (1960) derived (4) from a study of a simple model of turbulent 
diffusion proposed by Taylor (1922). In Taylor’s model, it is supposed that 
a marked fluid element moves with velocity v along the x-axis in either the 
positive or negative direction. It is further supposed that, in traversing a 
distance dx, there is a probability of dx// that its direction of motion is rever- 
sed, and therefore a chance 1—(dx//) that it remains unchanged.* Thus the 
velocity of the element will tend to persist in one direction for a characteristic 
time of the order of //v. 

Let R(x,t)dx be the probabilicy that, at a time ¢ after the release of the 
element from the origin, it lies within the interval (x, «+dx) and moves to 
the right. Let L(x,t)dx be similarly defined for motion to the left. Then it is 


easily shown that 


OR OR v 
Y ae ps wa: < _ 
(9) ot Ox l laa 
; OL @b 9 
(10) at ax 7 (R-L). 
Defining a diffusive flux J(x,t) by 
(11) J = v(R-L) 
and recalling that P(x,t) is simply given by 
(12) P=R+L, 
we see that 
a SO 
(48) ots x’ 
f oJ 2v 2 oP 
us) ni re 
On integrating equation (14) over all x we find that 
15 | 242] f" seas = 0 
(15) at 1 tJ_. x,t)dx = 0, 
1.e. 
(16) { J(x,t)dx = Ae*?t!* 


*As Taylor (1922) originally defined it, the model involved a lattice of points on the x-axis, 
each separated from its nearest neighbors by a fixed distance d. It was supposed that the 
direction of motion of a fluid element could be changed only at a lattice point, the probability 
of such an event being g. Following Goldstein (1951) and Bourret (1960), we have made a 
transition to the continuous case by letting d — 0, g > 0, and d/q — 1. 
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where 4 is constant. If we suppose that the motion of the element is initially 
to the right, then A = v and the left-hand side of (15) is /(t)/v, where M(t) 
is the Lagrangian velocity correlation. Thus 

(17) M(t) = ve 7"! 

Bourret has noted this fact, and has also observed that, in consequence, the 
solution to equation (14) satisfying the initial condition J(x,0) = 0 may be 
written in the form 


at 
(18) Fins -2{ M(t—0)P(x,")dt’. 


It follows therefore that 


( AP(xt) _ a” { 
(19) at Ox « 


t 
dt'M (t—t')P(x,t’), 


0 


an equation whose three-dimensional form is 
dP (x,t) a” fe 
20 See emo It’ M ,,(t—t’)P(x,t’). 
nr at nan ee 
As it stands (19) has been proved true only when A/(t) has the particularly 


simple form (17). However, Bourret (1960) made the hypothesis that the 
equation is generally valid. He has shown that, if P(x,t) satisfies (20), then 


et 
(21) (x x3) = 2 | (t—t')M,,(t')dt’, 


0 
a result which, as Taylor (1922) has shown, is exactly valid in actual turbulent 
diffusion. Also for a very general class of functions .V;,(t’), the solutions of 
(20) do, in fact, approach a Gaussian form with the correct variance ast — ©, 
3. ONE-PARTICLE ANALYSIS 
Equations (9) and (10) govern a particular case of a more general model 


defined by 


(22 ) aR _oR as 
ZZ ee Qa, 
OL OL 
2 oe ae Oe = 
i) at ax 


where the “collision terms’’ @ and 6 are linear functionals of R and L. These 
functionals are most easily expressed by a “‘scattering matrix’’ whose elements 


are linear operators: 


9 a Qi; Q12 R 
(24) = , 
B 2; Ae i 
Clearly, no matter what the instantaneous values of R and L may be, 


(25) a+ p=0, 
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i.e., the sum of any column of the scattering matrix is zero. It follows that 
(24) may be written 


(26) (:) ~ € 5G) 


For a system which treats right and left symmetrically, g = f, and so 


7) (3) = (4 AG). 


i.e. 
(28) a= —6= —f(J), 


where f is a linear functional. Clearly the case f = J// defines the Taylor 
model. Equations (22) and (23) may now be written 


(29) AR eR 


t = F lf m 
; ot Ox I) 


OL 6L 
aa a 
at ox 


(30) = f(J). 


Equations (11), (12), and (13) follow as before, but (14) is replaced by 


a/] »0P 
0 - 


(31) — +29f(J) = —v —. 
ct Ox 


We may ask what condition must be imposed on f(J) in order that (19) may 
hold. The answer appears to be that f(J) must be a local functional in space; 


i.e. we must have 
at 
(32) txt) = | E(t—t’) J (x,t')dt’, 


where E(t) is some kernel. It is easily verified that, when (32) holds, the relevant 
solution of (81) may be written in the form (18). It is clear that the model 
defined by (29) and (30) can include the effect of velocity persistence in a 
more general sense than Taylor’s model. More generally still, however, the 
collision term f(J) might include a non-localness in space. This might be 
interpreted as the effect of spacial velocity correlation in the motion of neigh- 
boring fluid elements. (In fact, this is certainly the case for one simple ‘‘two- 
particle’? model we shall describe presently.) Thus we may expect that a 
relationship of the form 


er ee 
(33) f(x,t) = | dt’ { dx’ E(x—x',t—t')J(x’,t') 


is more generally valid than (32). The relevant solution of equations (13) and 


(31) may then be written 


aP (x,t) Se Oh ght hs gs 
(34) ae = | dt’) | dx'O(x—x't—’) P(x’), 


ct Cx 
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where Q(x,t) is the Green’s function for the flux equation (31); i.e. it is the 
solution of 


a 0 : 

(35) 80 26f(@) = 6(x)6(t) 
which vanishes for negative ¢. On multiplying (34) by x? and integrating over 
all x, we find that 


ef Iw 
(36) (x") = 2 | dt’ (t—t’) | dx’Q(x’,t’), 
7 oY uw, 


( 


from which (cf. 21) we may conclude that 


(37) | > Q(x,t)dx = M(t). 


co 
Equations (19) and (34) share the property of predicting a Gaussian form for 
P(x,t) with the correct variance as t — , for a wide variety of kernels E. 
It is important to notice, however, that equations (19) and (34) are essentially 
different in character: equation (84), like equation (5), determines dP/dt 
in terms of a weighted average of P not only over previous times but also 


over neighboring positions. 


4, TWO-PARTICLE ANALYSIS 

We will now discuss a two-particle model in order to clarify the role played 
by spacial correlations. We will suppose that at time ¢ the probability that two 
fluid elements initially situated at xt and x» should lie in the intervals (es; 
x; t+dx,) and (x2, X2+dx2) is B(x1,x2:t\2x1, x2)dx,dx»; when no confusion can 
arise this probability will be denoted by B(x1,%2,t)dxidx-. As before we will 
suppose that the elements move along the x-axis with velocity +v. Let 
W(x1,X2,t)dx,dx» be the probability that at a time ¢ after their release the 
elements lie in the intervals (x1,x;-+dx,1) and (x2,x2.+dx,.) and are both moving 
to the right. Let X(x1,%2,t)dx,dx. be similarly defined for the probability that 
both are moving to the left; Y(x1,x»,t)dx,dx. for element 1 moving to the 
right and element 2 moving to the left; Z(x1,x2,t)dx,dx. for element 1 moving 
to the left and element 2 moving to the right. We allow for the effect of spacial 
correlation by supposing that the probability that the motion of either element 
is reversed in traversing a distance dx depends upon the position of the other 
element and on whether this element is moving in the same or in the opposite 
direction. 

It is clear that 


(38) W+X+Y+Z = B, 


and that the probability densities W, X, Y, and Z must be governed by 


equations of the type 


OW , OW . OW, 
—o, Ox, . OX» - 


(39) 


Sc: 
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OX OX dX 





(40) tte, "te * 
aY, aY a¥Y 
(41) ot e+ i g, 
: AZ | 
(42) OZ _ 92, 92 = y, | 


* Co 

ot Ox, 0X2 
where the “collision terms”’ &, 7, ¢, and yw are linear functionals of W, X, Y, 
and Z. These functionals are most easily expressed by a “scattering matrix” 
whose elements are linear operators: 


é Qi: Qi2 G13 Qi4 W 
43 n cS G21 G22 G23 Gea x 
(43) 7.7 : 

é 431 G32 33 34 } 

Me Q41 G42 G43 Dag Z 


Clearly, no matter what the instantaneous values of W, X, Y, and Z may be, 


(44) Et+ant+et+u=0, 
i.e., the sum of any column of the scattering matrix is zero: 
4 
(45) > ai = 0, (j = 1-4). 
i=1 


Also, since* 
(46) W(x1,x%2) = W(x"), X (x3,%2) = X(%9,%1), Y(x1,%2) = Z(%2,%1), 
and 
(47) &(x1,%2) = &(X2,%1), n(X1,X2) = n(X2,%1), €(x1,X2) = p(X2,%1), 
we must have 
(48) ai3 = Gu, a3; = Qa, 23 = Axz, Q3, = A42, 23 = 44, 
34 = Agp. 
It is reasonable to suppose that the probability of a simultaneous reversal 
of the motion of the two elements is negligibly small, i.e. 
(49) Ga = Gn = @a = Ga => OC 
For a system which treats right and left symmetrically, 
(50) 133 = A923, Q21 = Q32. 


It follows from equations (45), (48), (49), (50) that 


g —2a 0 b b W 

(51) n A Q —2a b b x 
c a a —2b 0 ye : 
ue a a 0 —2b zs 


*More precisely, (46) should be written W(x1,x2:t.x4,x2) = W(x2,x1:12%, x4), ete. i 








eee 
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where a and 6 are linear operators which are equal if and only if the motions 
of the two elements are uncorrelated. Using (51), the differential equations 
(39) to (42) may be written 








0 0 0 
529 — + y— +y—+ 9, a c = 
(52) | 24 _ re a |n b(Y+Z) = 0, 
0 0 
~e 2 i ae ) Fakes Fi = 
(53) es "ae, r+ 20 |x b(Y+Z) = 0, 
(54) & tee —0S +2 2 |y-aWw4x) = 0, 
1 
0 0 0 
55 - — Y——— 2 — = 
(55) ie ft? ra 2» |z a(W+X) = 0. 


In the simplest case in which a and 6 are constants, the solution of equations 
(52) to (55) presents little formal difficulty. In this case the spacial velocity 
correlation between the motions of the elements is independent of their 
separation and may be found from the limiting forms of W, X, Y, and Z as 
t — o., It is easily shown from equations (52) to (55) that, for times large 
compared with a~!, 6-!, and |x1—x:s!/2, 


J 7 Waxdx, * | | Xdx dx» = b/2(a+b), 


j { ‘¥dxydxs 


It follows that the cia correlation between the velocities “; and ws of the 


(56) 


| | Zdx dx. = a/2(a+6b). 


two elements at simultaneous instants is 
(57) ( Uylle ) = v*(b—a)/(b+a). 


This vanishes if and only if a = 6. 

To determine \/(t), we suppose that the element | is initially moving to 
the right. Then, by equation (56), the probability that, at that time, element 
2 is moving to the right is 6/(a+6), and that it is moving to the left isa/(a+0d). 
Thus, we solve equations (52) to (55) under the initial conditions 


(58) jwe b5(x1—21)5(x2—x2) /(a+b), Y = a6(x,—x})8(x2—x>) /(a +b), 
= LX =0, Z=0. 


Since 
R(x1) = J [W(x1,%2) + ¥(x1,x%2) |dxo, 


(59) . 
L (x1) = x [XG \AZ (x1,X2 ») |dx2, 


it follows that 


(60) M(t) =o" a} [(W-X+ Y¥—Zdxidx». 








ROBERTS: TURBULENT DIFFUSION 1299 


On solving equations (52) to (55) under initial conditions (58), we find that 


(61) M(t) = v' [be +ae "| /(a+b), 


(which reduces to the form (17) if a = 8). 

To determine the joint probability function B(x1,x2,f) we must specify the 
initial values of W, X, Y, Z in accordance with the spacial velocity correlation 
(cf. 56), and must therefore solve equations (52) to (55) under the initial 
conditions 


W=X 
Y=Z 


b8(x1—x1)5(x2—2x2) /2(a+), 


) 
(88) ad(x1—x4)6(x2—x2) /2(a+d). 


From the solution derived in this way, we can form 


(63) P (x,t) = J B(x 1,X2,t)dxe 


and J(x;,t) from equations (11) and (59). We find that these quantities obey 
equations (13) and (34), and that the Laplace transform of Q(x,t) is given by 


. tg an TE ia ie 
(64) QO(x,p) = J QO(x,be "dt = p+4ab/(a+b) 
_2ab(p+2a+2b)(b—a)'e7'"" 


yip+4ab/(a+b)]|*(b+a)> ’ 


where 
(65) y? = (p+2a)(p+2b)(p+2a+2b)/[p+4ab/(a+b)]. 


It is easily verified from equations (61) and (64) that Q(x,t) satisfies (37). It 
is important to observe that Bourret’s equation (19) is satisfied if, and only 
if, the motion of the two elements is uncorrelated, for it is in this case and 
in this case only that Q(x,t) as given by equation (64) assumes the simple 
form 


(66) Q(x,t) = M(t)6(x). 
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UNE ETUDE DU PHOSPHORE ROUGE PAR 
SPECTROMETRIE DE MASSE! 


JeAN-DENIS CARETTE ET LARKIN KERWIN 


RESUME 

Dans la premitre partie de cette communication, nous décrivons une source 
pour échantillons solides, dont la température peut étre rapidement variée de 
100° C a 800° C et contrélée 4 +3° C prés. Une analyse du lithium avec cette 
source donne 12.42 pour Je rapport Li?/Li®, et indique l’absence de discrimination. 
Dans la deuxiéme partie, nous présentons les résultats d’une analyse du phosphore 
rouge. Nous avons déterminé les abondances des ions P*, Pz, P3, Pg et Ps en 
fonction de la température des électrons ionisants, leurs courbes d’ionisation, des 
mécanismes de formation, des valeurs énergétiques, l’identification des ions 
doublement chargés, des ions négatifs, des bandes d’Aston et les chaleurs de 
sublimation. Ces résultats indiquent que toutes les espéces subliment directe- 
ment du solide. 


I. UNE SOURCE POUR ECHANTILLONS SOLIDES 

1. Introduction 

L’analyse des échantillons solides en spectrométrie de masse remonte a 
Dempster (1918) qui bombarda la surface de l’échantillon avec un faisceau 
d’électrons. Comme exemple d'utilisation récente de cette méthode, on peut 
citer Plumlee et Smith (1950). Aussi, trés t6t, Aston (1921) analysa des 
échantillons solides en les placant dans l’anode d’une décharge a gaz. Dempster 
(1935) perfectionna la décharge, et Gorman, Hipple et Jones (1951) décrivirent 
une technique pour éliminer les difficultés d’instabilité qui lui sont inhérentes. 
Une troisitme méthode d’échantillonnage consiste 4 évaporer des ions d’un 
filament sous des conditions rigoureusement contrélées. Inghram et Hayden 
(1954) décrivirent plusieurs de ces sources. Une autre méthode, qui est souvent 
employée, utilise un four pour évaporer l’échantillon. Dans certains cas, la 
vapeur sublime ionisée; dans la plupart des cas, elle est ionis¢e par bombarde- 
ment électronique. L’article de Chupka (1959) donne un excellent exemple 
de cette méthode. 

Depuis un an, nous utilisons une source a vaporisation dont la température 
peut ¢tre rapidement variée entre les valeurs de 100° C et 800° C et qui est 
stabilisée 4 +3° C. La manipulation de la source est simple, et les échantillons 


peuvent ¢tre changés rapidement. 


2, Description 

La figure 1 indique schématiquement le fonctionnement de la source. Un 
creuset S est chauffé par un filament F qui est alimenté par le bloc B. La 
force électromotrice du thermocouple T mesure la température du creuset. Elle 
est comparée en C avec une force ¢lectromotrice de référence R proportionnelle 
a la température désirée. La différence entre les forces électromotrices est 
amplifiée et réagit sur B jusqu’a ce que le creuset soit a la température 
voulue. 


1Manuscrit regu le 2 mai 1961. 
Contribution du Centre de Recherches en Physique, Université Laval, Québec, Qué. 
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Fic. 1. Diagramme de la source. 


(a) Le creuset S 

Le creuset et son support sont indiqués dans la figure 2 dessinée a 1l’échelle 
indiquée. Un manchon de porcelaine sert de support, étant lui-méme fixé 
dans un mandrin de laiton qui peut ¢tre rapidement inséré dans le spectrométre 
de masse. Le creuset est usiné a partir de cylindres d’Al.O; fournis par l’Ameri- 
can Lava Corporation avec une pureté de 99%. Cette porcelaine n'est pas 
poreuse. On trouve une cavité a chaque extrémité. L’une, destinée a recevoir 
léchantillon, est fermée au moyen d’une vis d’acier inoxydable, elle-méme 
percée d’un orifice de 0.5mm servant a laisser sortir les vapeurs. L’autre 
cavité sert de réceptacle pour le thermocouple. Le contact thermique entre ce 
dernier et le creuset est obtenu au moyen de porcelaine appliquée sous forme 
liquide et ensuite séchée. Le tout a une masse de 15 grammes. Le manchon 
glisse sur le creuset et dans le mandrin, permettant un assemblage rapide. 

(b) Le filament F 

L’élément de chauffage est un filament de wolfram de diamétre 0.05 cm et 
de longueur 25 cm. Il a une résistance d’environ 0.2 ohm a froid. Lorsque 40 
watts environ sont dissipés dans ses 12 tours, il peut maintenir le creuset a 
une température de 800°C. Le contact thermique entre le filament et le 
creuset est assuré par de l’alundum et un agent agglutinant de marque 
Sauereisen (Thinner No. 14). 

(c) Le bloc d’alimentation B 

Le filament F est alimenté par un transformateur Hammond 1144R60, 
fournissant de 0 4 10 volts, 60 c.p.s. Le primaire de ce transformateur est 
branché sur un autotransformateur, modéle Variac 250 VA, lui-méme branché 
sur une ligne de 120 v non stabilisée. Le bloc fournit de 0 a 100 watts au 
creuset, la variation du Variac étant manuelle ou encore contrélée par le 


bloc A. 
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Fic. 2. Le creuset et son support. 


(d) Le thermocouple T 
Celui-ci est constitué d’une jonction chromel—alumel formée de fil numéro 18, 


i 
j 
; 
i 
{ 
d’un diamétre de 1 mm. II fut calibré en regard des températures de fusion de | 
l’étain, du bismuth, du plomb et du zinc, en plus du point d’ébullition de 
l’eau. La force électromotrice fournie était 39+0.5 uwv/°C. Avec ce thermo- 
couple, nous pouvions faire des mesures relatives de température de +1° C, 
et des mesures absolues de +15° C. 
(e) La f.e.m. de référence R | 
Une cellule 4 mercure Mallory No RM4R fournit 1.34 volts. Elle alimente 
un diviseur de tension composé de sept résistances de précision de 100 ohms 
chacune, d’un potentiométre de 100 ohms et d’une résistance de précision de 
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32.7 kilohms. Dans le total de 33.5 kilohms, il circulait alors 39.5 wa. La chute 
de tension a travers chaque résistance de précision était donc 3950 uv, équiva- 
lent a la f.e.m. générée par une différence de 100° C dans le thermocouple. Les 
résistances et le potentiométre de 100 ohms étaient alors calibrés directement 
de 0 a 800° C. 

(f) L’élément de comparaison C 

La f.e.m. générée par le thermocouple et la différence de potentiel fournie 
par le diviseur de tension sont placées en opposition, tel qu’indiqué a la 
figure 1. Ainsi, lorsque le diviseur de tension était ajusté a ‘‘520°C’’, si le 
thermocouple ne fournissait pas une f.e.m. correspondant a cette température, 
une différence de potentiel résultante se manifestait 4 l’entrée de |’amplifi- 
cateur A, 

(g) L’amplificateur A 

Le signal d’erreur fourni par C est converti par un vibrateur en voltage 
c-a., dont la phase dépend de la polarité du signal. Un amplificateur (e.g. 
Minneapolis-Honeywell No 40XZ-2200) permet a ce signal d’actionner un 
servomoteur dont le sens de rotation dépend de la phase du signal. Ce moteur, 
aprés une démultiplication convenable (1 4 80) contréle la puissance fournie 
a l’élément de chauffage F. 
3. Caractéristiques 

Une fois que la température désirée est indiquée sur les boutons de mani- 
pulation de R, le systéme porte le creuset a cette température dans un délai 
de 60 secondes environ. La température est ensuite maintenue 4 +3° prés. 
Dans la figure 3, nous voyons le comportement du systéme suivant une série 
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Fic. 3. Comportement thermique de la source. 
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de commandes. Cette figure indique la f.e.m. générée par le thermocouple 
en fonction du temps (l’échelle du temps est repliée 4 zéro 4 environ toutes 
les 25 minutes). La partie (a) indique la réponse a une série de commandes 
de températures croissantes de 100°, la partie (b) 4 une série d’augmentations 
plus faibles telles qu’indiquées et la partie (c) indique la stabilité sur une 
période étendue a 800°C. On remarquera: la réponse rapide (environ 1 
minute/100°), le plus petit intervalle commandé (5° C), l’absence d’oscillation 
et la stabilité (43° C). Au besoin, une étape d’amplification supplémentaire 
dans A donnerait une sensibilité susceptible de fournir une stabilité de +0.5° C. 


4. Discrimination 

Un but important d’une analyse au spectrométre de masse est la mesure 
des abondances relatives des masses présentes dans l’échantillon. Or, si la 
vapeur issue d’une source telle que nous venons de décrire doit reproduire 
fidélement la composition de |’échantillon solide, il faut éviter la discrimina- 
tion de la source envers les différentes masses. Les causes de discrimination sont 
variées. Les principales sont: l’évaporation, l’écoulement effusif des vapeurs, 
la section d’ionisation et la différence de pression entre les plaques accéléra- 
trices de la source. Ces effets ont été discutés par Stevenson (1949) et Cogges- 
hall (1944). Les deux premiéres causes, généralement les plus importantes, 
ont des effets opposés, proportionnels a la racine carrée de la masse. Ainsi, 
dans l’absence d’effets supplémentaires tels que absorption préférentielle par 
la surface de la source, effusion moléculaire, etc., un creuset ne devrait pas 
produire de discrimination si les amplitudes des deux effets sont égales. 

Afin de vérifier cette caractéristique, nous avons effectué une série de 
mesures sur l’abondance relative des isotopes du lithium. L’écart entre les 
masses des deux isotopes (6 et 7) de cet élément étant relativement grand et 
leur abondance relative étant bien connue, cette expérience constituait un 
bon essai du systéme. II s’agissait alors de déterminer: 

(a) si, 4 une température donnée, il y avait discrimination au cours de la 
vie d’un échantillon; 

(6) s'il y avait discrimination, quelle que soit l'histoire de l’échantillon, 
c’est-a-dire qu’il ait été chauffé a n’importe quelle température pendant des 
temps différents. 

A cette fin, nous avons exécuté quatre séries de mesures: 

(1) analyse de l’échantillon jusqu’a extinction; 

(2) analyse sur des paliers de températures croissantes; 

(3) analyse sur des paliers de températures décroissantes; 

(4) analyse sur des paliers de températures variant de fagon quelconque. 

L’échantillon était chauffé 4 200° C pendant 2 heures afin d’éliminer l’eau 
de déshydratation du lithium. Cette élimination était contrélée par l’analyse 
continue des masses. La source ne fournissait pas de fond appréciable autre 
que celui occasionné par l’eau qui disparaissait aprés l’intervalle cité. On 
augmentait alors la température du creuset pour l’analyse du lithium. 

Les résultats sont indiqués dans le Tableau I ot chaque chiffre est une 
moyenne d’environ 15 mesures. Les numéros en téte correspondent aux séries 


d’expériences précédemment mentionnées. 
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TABLEAU I 


Mesures sur le rapport d’abondances des isotopes du lithium 


Type d’expérience 


Temp. (1) (2) (3) (4) 
300° C 12.45 12.47 12.45 12.43 
350° C 12.44 12.43 12.43 12.46 
400° C 12.47 12.438 12.45 12.41 
450° C 12.30 12.36 12.34 12.36 
Moyennes 12.42 12.43 12.42 12.41 


La moyenne de ces quelque 200 lectures est de 12.42+1 pour le rapport 
des isotopes Li?/Li® L’écart maximum de la valeur moyenne d’une série 
(a, A 450° C) est 0.12, soit 1°%. La valeur moyenne de plusieurs autres travaux 
(Ordzhonikidze et Shyuttse 1956; Friedman 1955) est 12.48+1. Notre valeur 
est alors trés satisfaisante, et le probléme de la discrimination ne se posait 
pas pour notre source. . 

Il. ETUDE DU PHOSPHORE ROUGE 
1. Introduction 

Le phosphore n’a pas été l’objet de beaucoup d’études par les spectro- 
métristes de masse. Aston (1920) a établi que Il’élément était monoisotopique 
a 3% prés. Nous avons vérifié l’'absence d’autres isotopes (Kerwin 1954) et 
avons établi la formation des ions P*, Pf, Pf et Pf par bombardement 
électronique des vapeurs du phosphore rouge. Nous avons également rapporté 
l’existence de l’ion Pf, ce qui confirmait les conclusions de Pauling et Simonetta 
(1952) quant au mécanisme de formation du phosphore rouge a partir du 
phosphore blanc (le rouge serait composé de chaines de tétrahédres ‘blancs’, 
et serait alors susceptible d’exister en blocs de 4, 8, 12... atomes). Kane et 
Reynolds (1956) ont repris l'étude du phosphore rouge ainsi que de I’arsenic. 
Ils ont vérifié la présence des ions P+, Pf, P} et P{ dans le spectrométre, et 
ont mesuré les courbes d’ionisation pour ces espéces. Ils ont aussi mesuré 
les chaleurs de sublimation, ce qui indiquait que seule la molécule P, était 
vaporisée du solide, les ions fragmentaires étant le résultat du bombardement 
électronique et de l’ionisation du dimére produit par l’action du filament 
chaud sur le tétramére. Malgré une recherche assidue, Kane et Reynolds n’ont 
pas pu confirmer la découverte du P{ quoique lion analogue Asf{ ait été 
observé. 

Nous avons repris l'étude du phosphore, utilisant la source 4 température 
contrélée décrite dans la section précédente. Nous avons utilisé le phosphore 
rouge commercial, mais obtenu de différents fournisseurs et de grades variés. 
Tous les échantillons ont donné des résultats semblables. 


2. Les abondances relatives des ions 

Dans les figures 4, 5, 6 et 7, nous reproduisons les mesures effectuées sur 
les abondances relatives des différentes espéces d’ions obtenus en bombardant 
les vapeurs issues du phosphore rouge par un faisceau d’électrons d’énergie 
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variable. Dans les titres des figures, on trouve les conditions utilisées. Chaque 
courbe représente une moyenne de 20 séries de mesures. 

Ces courbes sont en accord qualitatif avec les quelques données que nous 
avons déja rapportées (Kerwin 1954) et avec Kane et Reynolds (1956). Nous 
pouvons tirer d’une étude de ces courbes les conclusions suivantes: 

(i) Les ions obtenus 4 une température inférieure 4 300°C sont issus de 
vapeur venant de l’échantillon 4 une certaine phase. Au-dessus de 300° C, 
cette phase dans |’échantillon solide changerait. Passé 350° C, les faisceaux 
d’ions sont absorbés dans notre spectrométre. 

Cette conclusion est appuyée par les mesures faites en fonction de tem- 
pératures croissantes et décroissantes, et par l’absence d’effet de dimensions 
de l’échantillon. Plusieurs mesures faites avec des échantillons de diverses 
masses ont éliminé la possibilité d’effets de rapport masse/surface. C’est aussi 
a la température 350° que paraissent les bandes d’Aston, mentionnées plus 
loin. 

(ii) Le taux de croissance des différents ions (figure 5) est différent pour 
chaque ion, d’une région de température a une autre. Ceci implique nécessaire- 
ment que les ions autres que Pf ne sont pas entiérement formés par dissociation 
de ce dernier, autrement les taux de croissance seraient ceux du Pf. Cette 
conclusion, qui contredit celle de Kane et Reynolds, est appuyée par deux 
autres faits: (a) les mesures sur les potentiels d’ionisation, décrites plus loin, 
et (0) l’analyse du rapport: 

y= PLP, 
en fonction de la température. Ce rapport devrait décroitre avec la tempéra- 
ture si les ions P*, P}, P} étaient formés du Pf, étant donné que la dissociation 
de ce dernier devient plus probable avec l’augmentation de température. Mais 
au contraire, les graphiques résultant de ces calculs indiquent plutét une 


croissance (figure 8). 











a - ay 


300 350 400 450 





TEMPERATURE °C 


Fic: 8. Graphique de y en fonction de la température. 
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En résumé, les espéces moléculaires autres que le Ps vaporisent du solide 
qui passe par au moins deux phases, dont la température critique se situe 
vers 300° C. 


3. Les courbes d’ionisation 
Dans les figures 9, 10, 11 et 12, nous reproduisons les mesures de I’intensité 
des courants des ions P+, P}, Pf et Pf en fonction de l’énergie des électrons 
ionisants. Ces courbes sont en accord général avec les mesures de Kane et 
Reynolds, a l’exception du détail observé au début de chaque courbe. Nous 
appellerons ces graphiques des ‘‘courbes d’ionisation”’. Quelques “‘bris’’, bien 
définis mais difficiles 4 reproduire a |’échelle sur un seul dessin, ne sont pas 
indiqués dans les figures. 
Courbe d’tonisation du Pt 
Nous avons mesuré des “‘bris”’ a 14, 16, 18, 20, 23.5 et 26.5 ev a pression 
basse, et aussi a 11 ev a pression moyenne. Une considération de ces bris 
nous permet de proposer le systéme cohérent de réactions suivant: 
11 ev: P +e-— P* + 2e, 
14 ev: P, +e — P* + P; + 2e, 
16 ev: P, +e — Pt + P + 2e, 
18 ev: Py t+e—P*t+P + P, + 2e, 
20 ev: Py +e — Pt + Pi + 2e, 
23.5 ev: P,; +e — Pt + 3(P) + 2e, 
26.5 ev: P4 +e — Pt + 3(P)* + 2e. 


Dans ce systéme de réactions, nous utilisons les données suivantes qui sont 
connues: 


I(P) = llev (Moore 1949), 
D(P2) = dev (Herzberg 1950), 
D(P,— 2P2) = 2.3 ev (Cottrell 1954; Goldfinger et Jeune- 
homme 1958). 
On obtient par contre la valeur suivante, qui est nouvelle: 


D(P,— P3+P) = 3ev. 


La substitution de ces valeurs dans les équations des réactions permet de 
vérifier leur concordance avec les bris observés et leur cohérence. I] est a 
noter que les bris 4 20 ev et a 26.5 ev peuvent étre expliqués par une variété 
d’états excités. Ce sont les réactions 4 18 ev et a 20 ev qui contribuent le 
plus a la formation du P+ (45% et 50%). Les réactions 4 23.5 et 26.5 ev 
contribuent 5%, les autres moins que 0.1%. 

Courbe d’ionisation du Pt 

Dans ce cas, l’étude des bris permet de suggérer les réactions suivantes: 


10 ev: P3+e— Py + 2e, 


12 ev: P2 +e—-P2 + 2e, 

13 ev: P, +e — P{* + 2e, 
14.3 ev: Py +e — Pz + Po + 2e, 
15.5 ev: P, +e — PT + P32 + 2e. 
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Fic. 9. Courbe d’ionisation du P*. 
Fic. 10. Courbe d’ionisation du Py. 
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Nous utilisons les valeurs citées au sujet de P*, et tirons les valeurs suivantes 
qui sont nouvelles: 


EPs) = ey, 
D(P3— Pt+P) = 4ev. 


Le bris 4 10 ev ne se produit qu’a pression élevée. La détermination de 
D(P}) vient de l’équivalence J(P2)+D(Pf) = D(P2)+J(P), et donne une 
valeur (4 ev) quiest | ev inférieure 4 D(P2), chose typique des ions diatomiques 
du groupe V: Nf, As{ et Sbf. La valeur de 12 ev pour J(P2) vient des mesures 
a faible pression, et trouve de l’appui dans la valeur trouvée pour D(Pt) 
ainsi que d’une courbe empirique de Mitra (e.g. Kiser 1960) pour le trio 
I(N2), Z(Asz) et Z(P2). Les bris A 10, 12 et 13 ev établissent le P2 comme 
produit direct de l’évaporation. 

Courbe d’ionisation du Pt 

Dans ce cas, l’étude des bris nous permet de suggérer les réactions suivantes: 


11.5 ev: P; +e — Pf + 2e, 
13.2 ev: P; +e— Pj* + 2e, 
14.5 ev: P,; +e —P; + P + 2e, 


17.5 ev: P, te — P{* + P + 2e. 


En utilisant les données déja citées dans ces réactions, nous trouvons les 


valeurs nouvelles: 


I(P3) = ERS ev, 


D(Ps— P3+P) = 3 ev. 


Cette derniére confirme le résultat obtenu avec les bris du P*. Les bris a 
11.5 et 13.2 ev confirment |’évaporation du P; comme tel du solide, quoique, 
selon les pentes relatives des bris, l’ion P} vient surtout de la dissociation 


du P4. 


Courbe d’ionisation du Pt 
Ici, nous avons trouvé plusieurs bris, et nous suggérons la réaction suivante: 


Yev: P, +e—P; + 2e 


en accord avec Kane et Reynolds. Les autres bris ne permettent pas une 
interprétation définitive, plusieurs états excités de l’ion produit se prétant 
aux valeurs observées, soit 10, 11.5, 12, 12.5, 13.7 et 14.2 ev. 

En résumé, l'étude des courbes des figures 9, 10, 11 et 12 nous permet de 
conclure que le phosphore rouge s’évapore comme dimére et trimére, ainsi 
que sous la forme du tétramére, et probablement comme monomére. Nous 
avons confirmé plusieurs valeurs énergétiques et nous avons établi de nouvelles 


valeurs dans quatre cas. 


4. La molécule Ps 
En chauffant un échantillon de phosphore rouge 4a 300°C, un faisceau 


d’ions apparait a la masse 248, et persiste pendant des heures, jusqu’a l’ex- 
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Fic. 11. Courbe d’ionisation du Py. 
Fic. 12. Courbe d’ionisation du P;. 
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tinction de l’échantillon. Au début, il est accompagné d’autres ions qui sont 
attribués a des impuretés de surface (oxydes); ainsi 


236: P4sO;, 284: PsOro, 
252: POs, 299: PsOz, 
268: P4QOsg, ot: PsOr10. 


Il n’y a pas de doute que I’ion a 248 est le P$. Beaucoup d’expériences ont 
été faites pour vérifier cette question, soit: 

(i) Videntité de la masse 248 a été établie par référence 4 Pf} a la masse 
124, en changeant le champ magnétique, en changeant le champ électrique 
et en interpolant les masses d’autres ions connus; 

(ii) l’abondance de l’ion 248 suivait fidélement (a l'exception de la remarque 
a suivre) le comportement de |’échantillon; 

(iii) on a établi qu’il n’existe pas de réaction pouvant donner une masse 
a 248 par une bande d’Aston (voir section 6); 

(iv) on a établi que le faisceau ne pouvait pas parvenir de composés du 
phosphore, e.g. P2W. ; 

Cependant, cet ion n’a pas été observé par Kane et Reynolds. On remarque 
premiérement que l’énergie des électrons utilisés par eux (45 ev) leur cofitait 
un facteur 2 dans la sensibilité. Mais il semble évident que c’est leur traitement 
thermique de |’échantillon avant l’analyse qui a décomposé les chaines Ps 
en tétraméres: ils évaporaient le quart de leurs échantillons avant de procéder 
a l’analyse. Deux expériences nous ont permis de confirmer la chose. Premiére- 
ment, nous avons évaporé une partie d’un échantillon avant une analyse. Nous 
ne pouvions alors distinguer le P} dans le spectre, en accord avec leur expéri- 
ence. Deuxiémement, nous avons étudié l’abondance relative du P* en fonction 
du temps, a différentes températures. 

A 320° C, l’'abondance du Pt était de l’ordre de 1 partie dans 500 du Pj, 
et persistait jusqu’a extinction. Mais a 400° C, l’ion était moins abondant et 
diminuait A rien aprés } de vie de |’échantillon (environ 2 heures). A des 
températures critiques légérement inférieures 4 320° C, on observait parfois 
une abondance de 1%, en accord avec nos lectures datant de quelques années. 

Dans leur ensemble, ces résultats et ceux de Kane et Reynolds confirment 
la théorie de Pauling et Simonetta qui postule que le phosphore rouge est 
composé de chaines de tétraméres, dont les fragments auraient des masses 
de 4, 8, 12... atomes. Un traitement thermique élevé détruirait les liens 
entre les tétraméres. Nous avons aussi cherché un ion a la masse 372, i.e. 
P+,. Nous avons plusieurs fois observé un signal a cette masse, mais son 
abondance était faible (1/5000 de Pf), et sa durée trop courte pour que |’on 
puisse établir son identité avec certitude. Cependant, ce comportement est 
qualitativement ce que l’on attendrait de P},. Le potentiel d’apparition du 
Pf est 21 ev, et le faisceau est un maximum pour une énergie des électrons 


ionisants de 80 ev. 


5. Les tons négatifs 
Dans le spectre de masse des ions négatifs, nous avons observé P-, Py et 
P;, dont les abondances relatives étaient respectivement 1.5:10:1.5. L’abon- 
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dance de P; était de 5 parties sur 10‘ par rapport 4 P}. Nous n’avons pas 
observé de P7, la limite supérieure de son abondance étant de | partie sur 
10°. Ces valeurs ont été mesurées a une pression de 5X10-*mm Hg. Ces 
trois ions négatifs ont été observés précédemment par Dukelskii et Zanberg 
(1952). 

La courbe d’ionisation de Py est donnée dans la figure 13. Les courbes 
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Fic. 13. Courbe d’ionisation du Py. 


pour P~ et P> sont similaires a celle de Py dans la région de 2 a 7 ev, mais 
les courants d’ions de ces espéces sont nuls ailleurs, contrairement a Py. 

L’'examen de ces courbes et une étude des mécanismes possibles de forma- 
tion d’ions négatifs a la lumiére du principe Franck—Condon nous aménent 
a la suggestion que dans la région d’énergie des électrons ionisants de 2 a 
7 ev, les trois ions négatifs sont formés par dissociation spontanée de l’ion 
(P;). Parmi les faits conduisant a cette conclusion, il y a l’abondance relative 
constante des trois ions en fonction de l’énergie des électrons, la forme réson- 
nante des courbes d’ionisation, la valeur des potentiels d’apparition (inférieure 
a 10 ev, et trés voisine pour les trois ions) et la vérification d’une valeur 
d’affinité électronique. 

En effet, ce dernier peut étre établi a partir de l’équation suivante: 


En = D(Pa— Poat+Pa)+Ee—Am 


ot E- est l’énergie d’apparition de l’ion négatif, D(Ps) est l’énergie de disso- 
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ciation du tétramére en divers produits, £, est l’énergie cinétique des produits 
(~ 0 au seuil) et Ay, est l’affinité électronique. A partir des courbes d’ionisation 
telles que celles de la figure 13, et des valeurs obtenues de la section 3, nous 
trouvons: 


A(P”) = 0+0.5 ev, 

A(P:) = 0.3+0.5 ev, 

A(P3) = 1+0.5 ev. 
La valeur pour P~ est d’accord avec celles de Branscomb (1957) et de 
Bates et Moiseiwitsch (1955). Les deux autres sont nouvelles. La figure 14 
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Fic. 14. Formation résonnante de lion Pr. 


indique schématiquement le processus suggéré pour le cas de Py. L’ion PZ 
est aussi formé au-dessus de 8 ev par un autre processus résonnant et au-dessus 
de 12 ev par un processus non résonnant. Cependant, les ions formés dans ces 
régions dépendent de fagon critique de la pression. Nous suggérons alors que 
les ions P> formés a 9 ev le soient ainsi: 


P, +e Pz* 
tandis que ceux formés au-dessus de 12 ev le sont par dissociation du P, 
en paires. 
6. Les bandes d’Aston 


Lorsqu’un spectrométre de masse est opéré a pression relativement élevée 
(~10-4 mm Hg) il y a des collisions entre les ions du faisceau et les molécules 
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de l’atmosphére spectrométrique. Ces collisions entrainent des dissociations 
et des échanges de charges, ce qui donne un spectre de masse contenant de 
larges bandes. Nous sommes couramment a étudier ces bandes d’Aston dans 
notre laboratoire (e.g. McGowan et Kerwin 1960). Dans le cas du phosphore, 
nous avons considéré 30 réactions de collisions. La présence des bandes a six 
positions de masses non entiéres affirmait l’occurrence d’un certain nombre. 
Les bandes observées sont (avec les collisions qui peuvent y contribuer): 


7 5 ts 
Masse 7.5: pt > p+ 


Py > Py" 


Masse 10.3: Pp; — pi* 


Pz — Pt 


Masse 15.5: pt — p+ 
2 > Pe” 
Ry =F" 


Pi* > P;* 


pit pt+ 
Masse 23.3: pi — pt? 
Masse 41.3: py > P3 
Masse 139.5: P7* > P3 


Enfin, 18 réactions n’ont pas été observées dans les spectres de masse que 
nous avons obtenus. A une pression de 10-4 mm Hg, les bandes d’Aston qu’ils 
donneraient ont une abondance inférieure 4 10~° par rapport a Pf. 


P; — Pt+ ou P53", 

Pi* = Pi ou PFT <our Pe" cogs Pe, 
PP," — P;* ou Pf on Fs, 

P; — Pt, 

Pe =e Pt cou Pt ow Py 

pS — Ppt, 

P3* —P*+ ou Pt on P;*, 

ptt —» ptt, 


7. Les tons doublement chargés 

Les ions doublement chargés sont toujours beaucoup moins abondants que 
ceux portant une charge simple. Dans le cas du phosphore, les ions doublement 
chargés P}+ et P{*+ ont le méme rapport m/e que les ions P* et Pf, et ces 
derniers dominent le spectre de fagon 4 nous empécher d’observer les premiers, 
sauf indirectement par la courbe d’ionisation et de fagon qualitative seule- 
ment. Dans le cas des ions P+*+*+ et P+ cependant, les rapports m/e non 
intégraux permettent une étude relativement sfire. Nous avons observé les 
deux avec des abondances respectives par rapport 4 P{ de 5/10* et 1/10* 
environ. 

Nous avons pu avoir aussi les courbes d’ionisation pour ces ions. Dans le 
cas de P**, nous observons des processus initiés par des électrons ionisants 
d’énergie 17.5, 38 et 55 ev. Nous proposons les réactions suivantes: 
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17.5ev: P* +e-—Pt* + 2e (contribution faible), 
388 ev: Ps te—>P2e +P +Pt* +3650, 
55ev: Py te Pt +2P+Pt* + 4e, 
ouP; te—PF +P + Pttt + 4e. 


Les énergies qui justifient ce choix de réactions sont tirées de la section 3. 
Nous n’avons pas observé de bris 4 30 ev, et alors la formation du P*+* a 
partir de l’atome neutre n’est pas établie. 


8. La chaleur de sublimation 
Lorsqu’il y a évaporation libre d’un solide, on peut démontrer (e.g. Kennard 
1938) que la chaleur de sublimation est donnée par: 


(1) a R In(I*/Jo) 
os 
rT 


ou. AH est la chaleur de sublimation en kcal/mole, 

R est 0.002 kcal/mole /degré absolu, 

I* est la pression de la vapeur sublimée (dans notre cas, proportionnelle au 
courant d’ions), 

I) est une pression de référence (dans notre cas, la partie la plus faible dans 
la partie linéaire de la courbe), 

T est la température de |’échantillon en degré absolu, 

To est la température de |’échantillon a la pression de référence. 

Dans la figure 15, nous voyons les valeurs de In (J+/Jo) pour les différents 
ions observés en fonction de la température. Les pentes de ces courbes nous 
donnent les chaleurs de vaporisation, qui sont indiquées. 

De ces courbes, nous tirons les conclusions suivantes: 

(i) La valeur de 50 kcal/mole pour le Pf} est en accord avec la valeur 
acceptée pour le phosphore (Kane et Reynolds). Etant donné que c’est le Pf 
qui domine le spectre de masse, cette valeur est aussi approchée de celle 
de l'ensemble. 

(ii) Les valeurs de 30, 42 et 37 kcal/mole pour le P+, P} et P} sont nette- 
ment différentes de la valeur pour P}, indiquant que ces derniers subliment 
directement du solide. Ceci confirme les conclusions de la section 3. 

(iii) Conformément aux résultats de la section 2, nous voyons que les quatre 
espéces changent de valeur de sublimation au-dessous de 300° C, ot la valeur 
est de l’ordre de 10 kcal/mole. 

(iv) Lorsque l’échantillon approchait la fin de sa durée, la chaleur de 
sublimation pour le Pf augmentait, souvent brusquement, a 70 kcal/mole. 
Ceci semble correspondre a l’élimination du Pf, surtout a de hautes tempéra- 
tures, i.e. 4 un changement marqué de phase. 

Nous avons effectué une autre série de mesures sur l’ensemble des vapeurs, 
i.e. la pression totale en fonction de la température de la source. Cette expé- 
rience vérifia que la sublimation d’aucune espéce n’est limitée par une autre, 
la valeur trouvée pour la chaleur de sublimation de l'ensemble étant 34 
keal/mole, inférieure 4 la moyenne pondérée de l'ensemble. Cette valeur 
semble indiquer une discrimination envers le P} de la part de la jauge. 
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Fic. 15. Chaleur de sublimation du phosphore rouge. 


Kane et Reynolds avaient trouvé la valeur 56 kcal/mole et pour Pf et 
pour Pt. Nous attribuons ceci peut-étre au fait qu’ils avaient employé des 
électrons d’énergie faible (11 ev). Ceci exige ordinairement un filament trés 
chaud, ce qui aurait augmenté la contribution au faisceau de P} venant de 
la dissociation du P34. 
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SUMMARY 
The first part of the paper describes a source for solid samples which may 
be quickly heated to a desired temperature of from 100° C to 800° C and main- 
tained to +3°C. An analysis of lithium gives a ratio of 12.42 for the 7/6 
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isotopes and establishes that the source does not mass-discriminate. The 
second part of the paper reports the results of an analysis of red phosphorus, 
including the relative abundances of P+, Pf, Pf}, P{, and P{ as a function of 
temperature and bombarding electron energy, the appearance potential curves, 
suggested formation reactions and energy values, the identification of doubly 
charged ions, negative ions, and Aston bands, and values for heats of vaporiza- 
tion. The results indicate that the various molecular species sublime directly 
from the solid. 
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THE SUPERCONDUCTING TRANSITION TEMPERATURE OF LEAD! 


J. P. FRANCK? AND DoucGtas L. MARTIN 


ABSTRACT 


The temperature of the superconducting transition has been measured for 
high purity lead samples from three different sources. The transition temperature 
for all samples agreed to within 0.001° and was determined with a gas thermom- 
eter to be 7.193° K+0.005°. The transition width for all samples was roughly 
0.004°. Some comments on the use of carbon radio resistors as thermometers in 
the temperature range 3° to 30° K are included. 


INTRODUCTION 

A recent communication from this laboratory (Pearson and Templeton 
1958) reported a value for the superconducting transition temperature of lead 
which was significantly lower than the generally accepted value. The present 
authors had available an accurate gas thermometer, constructed for use in 
connection with low-temperature calorimetry, and the lead transition tempera- 
ture has been redetermined to see whether Pearson and Templeton’s value 
could be confirmed. Our result is not in agreement with Pearson and Temple- 
ton’s figure but is in substantial agreement with earlier determinations of the 
transition temperature. 

At the same time we have shown that the transition temperature of high 
purity lead from three different sources is essentially the same. The width of 
the transition zone for all samples was about 0.004°, which is substantially 
less than reported in most earlier work. These results suggest that the super- 
conducting transition temperature of lead may be used as a secondary fixed 
point. 

Carbon composition radio resistors are commonly used as thermometers at 
low temperatures, and the results of a calibration of such a thermometer 
against the helium gas thermometer have been included as an indication of 
the deviations of a commonly used empirical interpolation formula from true 
temperature. 


EXPERIMENTAL 

General 

The measurements were performed using an adiabatic calorimeter assembly. 
A schematic view of the part of the apparatus at low temperatures is given 
in Fig. 1. The copper clamp calorimeter (described in more detail by Franck, 
Manchester, and Martin 1961) carried the resistance thermometer, a vapor 
pressure bulb, and a heater. The gas thermometer bulb was soldered to a 
cylindrical copper block which fitted into the clamp calorimeter thus providing 

1Manuscript received June 19, 1961. 

Contribution from the Division of Pure Physics, National Research Council, Ottawa, 
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Issued as N.R.C. No. 6466. 
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Fic. 1. Sketch of the apparatus. 


a very good thermal contact. The calorimeter and gas thermometer bulb 
were surrounded by an adiabatic shield. This made it possible to keep the 
calorimeter at any desired temperature to better than 0.001°. The assembly 
was cooled to 20° K using helium exchange gas. The exchange gas was then 
pumped off and further cooling was achieved by passing cold helium gas 
through the vapor pressure bulb which for this purpose was connected to two 
capillaries (for this technique see Martin 1955). Exchange gas was never 
readmitted during a run and the resistance thermometer was therefore in a 
high vacuum during each complete run. 

The cryostat consisted of three stages usually kept at the temperatures of 
liquid helium, liquid hydrogen, and liquid nitrogen. For measurements from 
20° K upwards, however, the lowest stage was kept at the temperature of 
liquid hydrogen. 


Gas Thermometer 

The bulb of the gas thermometer was made of copper, wall thickness 1.6 mm, 
volume ca. 20 cm*. The bulb was connected by a cupronickel capillary to a 
mercury in glass manometer of 19-mm bore and 1-mm wall thickness. The 
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bore of the capillary was 0.3 mm in the cooled part and 0.5 mm in the part at 
room temperature. The over-all volume of the capillary was about 0.41 cm, 
with 0.31 cm* at room temperature, 0.064 cm® in the liquid nitrogen bath, 
0.022 cm*® between the bulb and the adiabatic shield (i.e. virtually at the 
temperature of the bulb), and the remaining 0.014 cm* between the adiabatic 
shield and the nitrogen can. The manometer was used as a constant volume 
manometer by always setting the lower meniscus to very nearly the same 
mark. The dead volume in the manometer was on the average 0.61 cm*. The 
readings of both menisci were taken relative to the nearest millimeter mark 
on a standard meter bar in the same plane as the menisci, using a cathetometer. 
The drum head of the cathetometer was divided into 0.01 mm. For the lighting 
of the menisci an optical lighting unit with a cylinder lens which provided a 
beam of light parallel to the meniscus surface was used (Collins and Blaisdell 
1936). Thus reflections on the meniscus surface were avoided. Repeated 
readings of a fixed pressure usually agreed to within 0.02 mm. 


Vapor Pressure Measurements 

A small cavity (ca. 0.1 cm’) in the clamp calorimeter provided the bulb 
for the vapor pressure thermometer. It was connected by a capillary (0.8-mm 
bore in the lower part) to a mercury in glass manometer. Pressure readings 
were taken against a calibrated glass mirror scale to about +0.2 mm. The 
capillary ran through the vacuum space; it was, however, thermally anchored 
to the adiabatic shield and to all stages of the cryostat. Care was taken to 
avoid cold spots. In this connection it is important to take due account of the 
difference in vapor pressure at a given temperature of the partly converted 
hydrogen in the refrigerant bath and the hydrogen sample in the vapor 
pressure thermometer. 

Calibration temperatures were determined from the vapor pressures of 
hydrogen and helium. The vapor pressure of the hydrogen sample was com- 
pared near the normal boiling point with the vapor pressure of a sample of 
room temperature equilibrium hydrogen. Corrections due to ortho-para 
conversion were usually found to be negligible. It was possible to avoid any 
appreciable amount of conversion of the hydrogen by passing it sufficiently 
rapidly through the high pressure charcoal cleaner. 

Resistance Thermometer 

A 10-2 Allen-Bradley carbon resistor of $-watt rating was used. After the 
plastic coating was ground off, the resistor was covered with a layer of varnish 
and then varnished into a cavity on the clamp calorimeter. A few centimeters 
of each lead to the thermometer were also varnished to the calorimeter. 
Resistance was measured with the isolating potential comparator circuit of 
Dauphinee and Mooser (1955) to about 0.0005 Q, using a measuring current 
of 50 wa throughout. The uncertainty in resistance is equivalent to 2.5 mdeg 
at 30° kK, 1 mdeg at 19° K, and less for lower temperatures. 

Lead Samples 

Lead samples from three different sources were investigated. Sample 1 was 
Johnson—Matthey spectroscopically pure lead, sample 2 was Tadanac brand 
high purity lead, supplied by The Consolidating Mining and Smelting Company 
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of Canada, Trail, B.C. Samples 3 and 4 were zone-purified lead, made from 
Johnson—Matthey spectroscopically pure lead (this last material was kindly 
supplied by Dr. A. Rosenberg of the University of Toronto). Sample 4 was 
actually taken from the batch used by Pearson and Templeton (1958) in 
their determination of the transition point of lead. 

The Tadanac lead is claimed by the manufacturer to be 99.9999% pure. 
Melting studies on another batch of Tadanac lead support this figure (McLaren 
and Murdock 1960). A spectrographic analysis of the lead samples showed the 
purity of the Tadanac lead to be about 99.9999%, of the zone-purified lead to 
be about 99.999%, and of the Johnson—Matthey lead to be about 99.997%. 
These figures are supported by the value of the residual electrical resistance 
ratio (R4.2° x/R293° x), obtained by extrapolation of magnetoresistance mea- 
surements to zero field. For the Tadanac lead R4.2° x/Rag3° x isabout 4.6 X 10-5, 
for the zone-refined lead about 9.6 X10-°, and for the Johnson—Matthey lead 
about 2.5X10-* (Okumura, private communication). 

Samples 1 to 3 were in the form of a cylinder 20 mm long and 6 mm diameter. 
Sample 4 was made out of wire ca. 0.1 mm thick formed into a bundle 20 mm 
long and 3 mm diameter. All samples were soldered with indium-lead solder 
to the bottom of the gas thermometer bulb (see Fig. 1). 

Detection of the Superconducting Transition 

The transition was detected magnetically using a mutual inductance 
method. A primary and a secondary coil were wound close to one another on 
one end of each sample. Each coil had 10 turns. In series opposition with the 
secondary was an exactly similar coil wound on the opposite end of the samples. 
This arrangement was chosen to eliminate the pickup of stray fields in the 
detection circuit. The primary coil was operated at 800 c.p.s. and a field of 
about 0.05 oersted. An electronic amplifier was used for observing the signal 


in the secondary coil. 
CALCULATION OF TEMPERATURES FROM GAS THERMOMETER READINGS 


Temperatures 7 were calculated from 


Up u% 
(1) T(U1+B/V)* pd, T(1+B,/V,) = constant 
where p is the pressure, v) the volume of the gas thermometer bulb, B the 
second virial coefficient, and V the molar volume. The second term represents 
the correction for the dead volume, the summation extending over all parts 
of the dead volume. v,; is the part of the dead volume at temperature 7). 
The constant depends on the amount of gas in the system and is determined 
by calibrating at a known temperature. This calibration was obtained near 
the normal boiling point of either hydrogen or helium. The normal boiling 
point of room temperature equilibrium hydrogen was taken as 20.378° K as 
recommended by van Dijk and Durieux (1958). In order to interpolate to the 
actually measured vapor pressures we used the equation of Woolley, Scott, 
and Brickwedde (1948), which would lead to a normal boiling point of 
20.390° K, and subtracted 0.012° from the temperatures given by this equation. 
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Temperatures in the helium region were calculated from the 1958 temperature 


scale (Brickwedde et al. 1960). 
The molar volume V has to be calculated from the equation of state: 


(2) pV = RT (1+ B/V) 


where R is the gas constant. Solving this quadratic equation for 1/V and 
expanding into a power series leads to: 


(3) T B/V = (Bp/R){1—Bp/RT + 2(Bp/RT)?— ...}. 


Higher-order terms are negligible for the filling pressures used. Several sets 
of second virial coefficients for helium have been published. Kilpatrick, Keller, 
and Hammel (1955) give calculations of the virial coefficient which are in 
excellent agreement with experimental values below 4.2° K (Keller 1955). 
For higher temperatures, however, the calculated values of Kilpatrick et al. 
are increasingly too large (at least up to 20° K). Here the ‘‘adopted values”’ 
of Keesom (1942, p. 49) are in much better agreement with experimental 
determinations. Both sets of virial coefficients agree near 4.2° K. Therefore, 
at the present time, it seems best to use the data of Kilpatrick et al. below 
4.2° K and those of Keesom above this temperature.” It might be noted that 
corrections due to the third virial coefficient are quite negligible. 

The correction for the dead volume was calculated from the known dimen- 
sions of the capillary and manometer. The co-volume of the meniscus was 
taken from a formula by Kistemaker (1945). The temperature distribution 
along the part of the capillary between the adiabatic shield and the liquid 
nitrogen bath and between the liquid nitrogen bath and room temperature 
was assumed to be linear. This is sufficiently accurate as these parts of the 
dead volume were kept very small. Non-ideality corrections for the gas in the 
dead volume, which are represented by the term 7, B,/ V; in eq. (1), were found 
to be negligible. Further corrections were applied for: 

(i) The ambient temperature of the mercury column. 

(ii) Capillary depression. The table of Cawood and Patterson (1933) was 
used. 

(iii) Thermomolecular pressure difference. Values for this correction were 
taken from Keesom (1942, p. 119). The magnitude of this correction is rather 
small. For the lowest filling pressure used it amounted to 0.002° at 4° K. 


ACCURACY OF THE GAS THERMOMETRY 

The following possible causes of error were considered: 
(t) Absolute Value of the Calibration Temperature 

At the present time there is an uncertainty of about +0.01° in the normal 

*NOTE ADDED IN PROOF: A critical re-examination of all available virial coefficient data has 
shown that the Keesom ‘‘adopted values’’ are the best values below 4.2° K also (van Dijk, 
private communication). Such a change would not alter the results reported in the present 
paper. 
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boiling point of hydrogen and about +0.002° in the 1958 helium vapor pressure 
scale. Calibration either near 20° K or near 4° K therefore always leads to an 
uncertainty of 0.003° to 0.004° near 7° K. This uncertainty is not included in 
the error limit given below for the superconducting transition temperature 
of lead. 
(it) Determination of the Calibration Temperature 

An error of 0.2 mm Hg in the vapor pressure determination amounts to 
0.001° at the normal boiling point of hydrogen and 0.0003° at the normal 
boiling point of helium. Repeated determinations of the calibration tempera- 
tures agreed to within 0.001°. We believe therefore that errors due to this 
source are negligible. 
(111) Pressure Measurement 

The precision of this measurement is about 0.02 mm Hg. This represents 
0.0013° for the lowest filling pressure used (independent of temperature) and 
less at higher filling pressures. The accuracy might be slightly worse than this 
figure since the apparent position of the mercury menisci will be affected by 
changes in the lighting conditions and non-parallelism of the manometer 
tubes. As far as these errors are random they are minimized by taking readings 
at more than one filling pressure. 
(iv) Correction for Dead Volume 

A total error of 15% in this term, due to incorrect volume and temperature 
determinations, would lead to errors of 0.001° at 4° K, 0.002° at 7° K, and 
0.006° at 20° K. This error is independent of the filling pressure of the gas 
thermometer. We believe that these figures are upper limits for the error due 
to this cause. 
(v) Correction for Non-ideality of the Gas 

For a filling pressure of 646 mm Hg at 20° K this correction amounts to 
0.159° at 4° K and for the lowest filling pressure used, viz. 303 mm Hg at 
20° K, to 0.075° at 4° K. If the virial coefficients of Kilpatrick et al. were 
used for the calculation of this correction the corresponding figures are 0.168° 
and 0.079°. These figures give some idea of the possible uncertainty in this 
correction. 
(vi) Correction for Thermomolecular Pressure Difference 

In the present work this correction is so small that any errors are negligible. 

The major part of the possible error seems to arise from possible errors in 
corrections for non-ideality and for the dead volume and possibly in the 
pressure determination. The over-all accuracy may be judged from a compari- 
son of temperatures measured near 4° K with the gas thermometer (when cali- 
brated near 20° Kk) and the same temperatures derived from the 1958 helium 
vapor pressure scale (73s). Such a comparison is made in Table I| and it 
will be seen that the average difference is 0.002;°. In view of this comparison 
it should be quite safe to assign error limits of +0.005° to the determination 
of the superconducting transition temperature of lead. 
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TABLE I 
Comparison with 1958 helium vapor pressure scale 











bs. tk 
(gas thermometer (helium vapor 

Filling pressure calibrated against pressure : 
(mm Hg at 20° K) hydrogen vapor pressure) scale) Difference | 
646 4.1085 4.113, —0.0045 | 

4.109 4.113, —0.004, 

303 4.113, 4.113; +0.000; 

4.112, 4.113, —0.001; 

TABLE II 


Superconducting transition temperature of lead 





7S rk 
Filling pressure (calibration against (calibration against 
(mm Hg at 20° K) _ helium vapor pressure) hydrogen vapor pressure) 
646 7.1952 7.1875 
7.195; 7.187; 
303 7.192, 7.191, | 
7.192, 7.191, 











THE SUPERCONDUCTING TRANSITION TEMPERATURE OF LEAD 
Transition curves have been measured against the calibrated carbon 
resistance thermometer for all four lead samples and are shown in Fig. 2 in 
arbitrary units of mutual inductance vs. temperature. By taking points with 
both increasing and decreasing temperatures it was shown that no hysteresis 


Mutual Inductance (arbitrary units) 





7.180 7.185 7.190 TI95 7.200 7.205 


Temperature °K 


Fic. 2. Superconducting transitions in the four samples investigated. 
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occurred. The transition widths are about 0.004°. The transition temperature 
is taken as the mid-point of the curves as measured on the inductance axis. 
As the transition temperature for all samples agreed to within 0.001° a direct 
comparison with the gas thermometer was made only for sample 3. Two 
independent comparisons with the gas thermometer were made, using filling 
pressures of 303 and 646 mm Hg at 20° K. Readings of the gas thermometer 
were made within 0.001° of the transition temperature, the small difference 
being derived from the calibrated carbon resistance thermometer. For each 
filling pressure the transition temperature was determined independently by 
two observers, the results being given in Table II. Temperatures derived from 
the gas thermometer by calibrating near the normal boiling point of either 
hydrogen or helium are given separately in Table II. Taking an average for 
both calibration temperatures and both filling pressures employed leads to 
7.191,° K for the transition temperature. This has to be corrected for the 
effect of the earth’s magnetic field, which was measured just outside the 
cryostat during the experiment and found to be 0.30 oersted. Assuming that 
the magnetic field at the sample is the same, the correction, to a first approxi- 
mation, is about 0.001°. We therefore find that the superconducting transition 
temperature of lead in zero field is 7.193° K+0.005°. 


USE OF CARBON RESISTORS AS THERMOMETERS IN THE 
TEMPERATURE RANGE 3° TO 30° K 

It is well known that some types of carbon radio resistors are very suitable 
for use as thermometers in the above temperature range, but have the dis- 
advantage that the calibration often shifts on warming to room temperature. 
This may be overcome by calibrating at two or more temperatures every time 
that the thermometer is cooled, and using some suitable interpolation formula 
for obtaining intermediate temperatures. It is found that deviations of inter- 
polated temperatures, obtained in this way, from true temperature remain 
sensibly constant on successive coolings. Thus only a single calibration of a 
given carbon thermometer against a gas thermometer is required. 

Two widely used interpolation formulas are those of Clement and Quinnell 
(1952) and Clement (1955). The former is a three-constant formula and its 
use within this temperature range for a 10-Q, }-watt Allen-Bradley resistor 
has been studied by Morrison, Patterson, and Dugdale (1955), who give a 
graph of deviation from true temperature. The two constant formula of 
Clement (1955) shows, as might be expected, larger deviations from true 
temperature. However, only two calibration temperatures are required, for 
which the hydrogen and helium vapor pressure scales, in the sensitive regions 
near the normal boiling points, may be used. Figure 3 shows deviations from 
true temperature for the 10-2, 3-watt Allen-Bradley resistor used in the 
present work when interpolating with Clement’s formula: 


(log R/T,.)! =a+blogR 


where R is the resistance, 7, the interpolated temperature, and a and 6b are 
derived from the calibration temperatures. 
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Fic. 3. Deviation of the interpolated carbon thermometer temperature scale from true 


temperature (see text). 


DISCUSSION 


A survey of previous determinations of the superconducting transition 
temperature of unstrained bulk lead is given in Table III. With the exception 
of the value given by Pearson and Templeton there is satisfactory agreement 
between all the determinations within the sometimes rather large error limits. 
Pearson and Templeton did not publish their original experimental data but 
they have kindly made it available to us and there seems no doubt that the 
published figures considerably overestimated the accuracy of their work. 


TABLE III 
Superconducting transition temperature of lead 











Superconducting 

















*95 


transition Width of 
Authors Material temperature, °K transition 
Kamerlingh-Onnes Kahlbaum 99.99% 7.2 
and Tuyn (1922) 
de Haas, de Boer, and Kahlbaum 99.99% 7.193 0.02 
van den Berg (1934) 
Daunt (1939) Hilger 99.999% 7.22+0.03 
Ziegler, Brucksch, Extruded wires, 7.20+0.01 0.03 
and Hickman (1942) Baker and Co. 
Boorse, Cook, and National Lead Co. 7.22,+0.04* 0.05 
Zemansky (1950) 99.999% 
Pearson and Johnson—Matthey 7.175+0.005 0.001 
Templeton (1958) 99 .999% 
This work Johnson-Matthey 7.193+0.005 0.004 


&% confidence limits calculated 


estimates. 


99.999% 


from author's published data. Otherwise stated error limits are author's 
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The present results on different lead samples indicate that slight differences 
in impurity content do not alter the transition temperature appreciably. 
The transition temperatures of common lead and radiogenic lead, however, 
differ by about 0.028° (Hake et al. 1958), the superconducting transition 
temperature being inversely proportional to the square root of the average 
atomic weight. Fortunately the average atomic weight of common lead from 
different sources differs only by about 0.05%, despite larger variations in 
isotopic composition (see Russell and Farquhar 1960, p. 13). This difference 
would alter the transition temperature by less than 0.002°. 


CONCLUSION 


The superconducting transition temperature for high purity lead samples 
from three different sources is 7.193° K+0.005°. 


ACKNOWLEDGMENTS 


Dr. C. A. Swenson suggested to us the method we have used for making 
the nor:-:deal gas correction. We are also indebted to Mr. E. Green, Metrology 
Section, N.R.C., for the calibration of various scales, to Mr. R. M. Morris, 
Radio and Electrical Engineering Division, N.R.C., for the loan of the mag- 
netic field measuring equipment, to Mr. D. S. Russell, Applied Chemistry 
Division, N.R.C., for arranging the spectrographic analysis of the lead samples, 
to Dr. W. B. Pearson for the lead samples, to Dr. I. M. Templeton for the 
loan of his mutual inductance bridge, and to Mr. R. L. Snowdon for assistance 
with the experiments. 


REFERENCES 

Boorse, H. A., Cook, D. B., and ZEMANsky, M. W. 1950. Phys. Rev. 78, 635. 

BRICKWEDDE, F. G., vAN Dijk, H., DuRIEUX, M., CLEMENT, J. R., and LoGaN, J. K. 1960. 
J. Research NBS, 644A, 1. 

Cawoop, W. and Patterson, H. S. 1933. Trans. Faraday Soc. 29, 514. 

CLEMENT, J. R. 1955. Temperature, its measurement and control in science and industry, 
Vol. 2 (Reinhold, New York), p. 380. 

CLEMENT, J. R. and QUINNELL, E. H. 1952. Rev. Sci. Instr. 23, 213. 

Couns, S. C. and BLAIspELL, B. E. 1936. Rev. Sci. Instr. 7, 213. 

DaunT, J. G. 1939. Phil. Mag. 28, 24. 

DAUPHINEE, T. M. and Mooser, E. 1955. Rev. Sci. Instr. 26, 660. 

VAN Dijk, H. and DurtEux, M. 1958. Physica, 24, 1. 

FRANCK, J. P., MANCHESTER, F. D., and Martin, D. L. 1961. Proc. Roy. Soc. In press. 

DE Haas, W. J., DE Borer, J., and VAN DEN BERG, G. J. 1934. Commun. Kamerlingh 
Onnes Lab. Univ. Leiden, 2336. 

HAKE, R. R., MApotuer, D. E., and DEcKER, D.L. 1958. Phys. Rev. 112, 1522. 

KAMERLINGH ONNES, H. and Tuyn, W. 1922. Commun. Kamerlingh Onnes Lab. Univ. 
Leiden, 160d. 

Keesom, W. H. 1942. Helium (Elsevier, Amsterdam). 

KELLER, W. E. 1955. Phys. Rev. 97, 1. 

KILpaTRICK, J. E., KELLER, W. E., and HAMMEL, E. F. 1955. Phys. Rev. 97, 9. 

KIsTEMAKER, J. 1945. Commun. Kamerlingh Onnes Lab. Univ. Leiden, 268c. 

McLaren, E. H. and Murpock, E.G. 1960. Can. J. Phys. 38, 577. 

MartTIN, D. L. 1955. Phil. Mag. 46, 751. 

Morrison, J. A., PATTERSON, D., and DuGpALeE, J. S. 1955. Can. J. Chem. 33, 375. 

Pearson, W. B. and Templeton, I. M. 1958. Phys. Rev. 109, 1094. 

Russet, R. D. and FARQUHAR, R. M. 1960. Lead isotopes in geology (Interscience Publ., 
New York). 

Woo Ley, H. W., Scott, R. B., and BricKWEpDDE, F. G. 1948. J. Research NBS, 41, 379. 

ZIEGLER, W. T., BruckscH, W. F., and HickMAN, J. W. 1942. Phys. Rev. 62, 354. 








THE EFFECT OF AN EXPLOSION IN A 
COMPRESSIBLE FLUID UNDER GRAVITY! 


R. A. Ross 


ABSTRACT 

In this paper an investigation is made of the effect of an axially symmetric 
explosion at any depth in a semi-infinite, compressible, non-viscous fluid, acted 
upon by gravity. The explosion is represented by a line source of the form 
4(x)5(2—h)6(¢), where h is the depth of the source. An exact solution is given 
using the linearized theory. This solution is studied in detail by asymptotic 
methods, for the special case of a surface explosion. It is found that compres- 
sibility results in the gravity waves being propagated with a speed less than c, 
the speed of sound in the fluid. If x is the distance from the explosion and ¢ the 
time that has elapsed after the explosion, then for 43/3 < x/ct < 1 only “pre- 
cursor’’ waves are noticed at the point of observation. For x/ct < 44/3 large 
amplitude waves are present, similar to the waves predicted by the incom- 
pressible theory. 

INTRODUCTION 

An important problem concerning gravity waves on the surface of a fluid is 
the following. A semi-infinite, incompressible, non-viscous fluid acted on by 
gravity lies completely at rest. At a certain time a disturbance takes place 
at a point on the surface of the fluid; physically this would correspond to a 
stone being thrown into a still pond. The problem is to study the waves 
produced by the disturbance. The original analysis was done by Cauchy and 
Poisson, and the results of their work as well as later developments are dis- 
cussed extensively in the books of Lamb (1932) and Stoker (1957). 

A considerably more difficult problem is to study the waves when the fluid 
is assumed to be compressible. In the first part of this paper, with the use of 
the linearized theory, an exact solution is obtained for the effect of an explosion 
at any depth in the fluid. In the next part of the paper approximate results 
are obtained from this rather complicated exact solution, for the particular 
case of the explosion situated on the surface, using the method of steepest 


descent. 


FORMULATION OF THE PROBLEM 


Let z = ¢(x,t) be the free surface of the fluid for which z = 0 is the undis- 
turbed free surface. Let the x-axis be taken in the undisturbed free surface, 
and the z-axis pointing down into the fluid. Since the explosion will be in the 
form of a line source in the y-direction, the y-co-ordinate will be completely 
absent throughout. 

If v = the velocity of the fluid, g = the gravitational constant, ¢ = the 
velocity potential, v = —V@, c = the velocity of sound, and Q = the potential 
of other external forces acting on the fluid in addition to gravity, then it is 
shown in Longuet-Higgins (1950) that 


‘Manuscript received March 10, 1961. i 
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ab as a6_8 (1,3) (1,*) - dQ 
(1) ere te ar ar he 


If v is so small that v? is negligible (1) reduces to 


ad 22 a¢_—ss 


= oe, “V'o-g= Sa 
(2) ap ° Y o- 85, dt 


To the same order of approximation it is found that the boundary condition 


on the free surface becomes 


ao _  a¢ st 
(3) an = a on z=0 


and the displacement of the free surface is given by 


_ _1(9¢ 
@) a” (3 o. 


Our problem reduces to solving the differential equation (2), with a 
suitable choice of dQ/dt to represent an explosion. In addition, the solution 
must satisfy the boundary condition (3) and the initial conditions that ¢ and 
0¢/dt both vanish at ¢ = 0. It will be assumed that these conditions together 
with a boundedness condition at infinity give a unique solution to the problem. 
We shall choose 


(5) dQ/dt = c°6(x)6(2—h)6(t) 


to represent the explosion, where the right-hand side is the product of three 
Dirac delta functions. This corresponds to considering the effect on the fluid 
of a radial line impulse at (0,4) at time ¢ = 0. 


SOLUTION OF THE BOUNDARY VALUE PROBLEM 


If we define 


(6) $(x,2,p) = | -¢Méilaee dt, 


the Laplace transform of ¢, where p is a real variable sufficiently large to 
ensure convergence, and 


 —— sS 7 _— téz f — in(z—h) 

(7) o= on Ry dx E* ¢ dz, 

the double Fourier transform of ¢, we find from (2) and (5) that satisfies 
1 2b 8 a mz 

(8) (e944 banjo=57- 


Using the inverse of (7) and carrying out the n integration by residues, we find 


elt expligy—|s—hiA] 4p 


Ze 3 


(9) jin 


—oa 
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where 
(ese) 

eae 

We must set up a two-sheeted Riemann surface in the é-plane with branch 

points at 


2 2\3 
2. tage 8 
g= +i(248,) ‘ 


We take the branch cuts extending from these branch points out to infinity 
along the imaginary £-axis (Fig. 1). With these branch cuts we can show that 





Fic. 1. The two-sheeted Riemann surface with branch points at — = +7£, £1 = [(p?/c?)+ 
(g?/4c*)|4. The branch cuts are the imaginary £-axis from ig, to 1@ and from —ié to —i™, 
The point A corresponds to s = 0, and the point B tos = o, 

Re(A) > 0 everywhere on the upper sheet of the Riemann surface. Thus the 
integral in (9) is convergent. The general expression for ¢ satisfying the bound- 


edness condition at infinity and the initial conditions will be 


: S 1 z—h) | (°° exp[itx—|z—h|r 
= feel ep] fremlecle—a gy 
eg eaite re . 
+ = Ji ex (tx — Az)a(E)dE 


where a(&) is an arbitrary function of & The Laplace transform of (3) is 


9 


po= = on z = 0. 


(11) 85. 


Substituting (10) into (11) we determine a(é) and obtain 
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se *co aw , —_ ie— 
(12) $= exp| -2 aay Re J explir soi dé 


0 


1 (z—h) * exp[itx— (z+h)d] 





i Re [ “explity=G+h)Al 9 


1 
0 ‘otarrs” 
g 2 
To find ¢ we adopt the following procedure. We attempt te write the right- 
hand side of (12) as a Laplace integral, i.e. an integra) of the form 


J “ey (dt, 


where y(t) is a known function. Then we can say ¢ = y. This is the method 
used by Cagniard (1939) in an elastic wave problem. 

To carry out this method we deform the path of integration in the second 
and third integrals in (12) into the new path given by 


5 
x 22 2 2 . 
| B—(e+H) elles —x°—(z+h) | +is| £ (=+h)-+xp | 
:* - x’°+(s+h)? 
where s is a real parameter ranging from 0 to ~ (Fig. 1). The first integral is 
deformed into the path given by expression (13) with |z—h| replacing z+-h. 
In both cases the deformation can be carried out by Cauchy’s theorem with 
no contributions from poles and vanishing contribution from the circle at 





(13) 


infinity. 
Since s = 0 corresponds to 


f= | #,- (s+h) ale (+h)* | 


part of the new path of integration is along the negative imaginary axis. 
However, since the integrands are pure imaginary along the imaginary axis 
there is no contribution from this part of the path of integration. In terms of 
the parameter s the first two integrals of (12) become 


exp, —ps+= (s’—k})*” is 
(14) a exp| —p ( <4 Re — aan a 
: 25 2c" / 0 (s°—k})'”” 
1 [ is ] exp| — ps + (F—#)"* | ds 
—5, P| ~2 92? Re , ee ee 
where 


] 9 a 
k= : (x°+ (s+h)*)}, 


(x°+(s—h)’)}. 


Sl 


: 
f 
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We see that the first term in (14) vanishes if 0 < s < k,, and the second 
vanishes if 0 < s < ko. Thus the contribution at (14) to ¢ is 


‘4 (o--h) i bw 1s ot. ays 
1 exo} —  2¢? | cos| (t'— ha) |ue-a 








(15) 3s Hy)” 
1 exp| —g —) cos £ (tay | H(t—k2) 
Qa (t°—k2)'”” 
where 
H(t)=0, ¢t<0, 
= 4, t>0 


The third integral in (12) becomes 


@ 


F\(s) exp| -ps—@ (s*—ki as 
(16 a e oe Re Se te A oe i 
) og XP & 2c" F kt [p+Gi(s)](s° — ki)” 
Be ees _ gt 2p 24/9 
oa | - ep] a | FG) 000 pf ED i 
a ae i [p+Ga(s)l(s°— kt)” 


where 


l 
bol — 
es, 
8 
oa 
tol 
Ie 
| &, 
| + 
iS! 
Pa 
Io | 


F\,0(s) a 


1/2 - 
(17) Gi 2(s) = gl (w'+4) +(w-‘) : 


ixs — (2+h)(s*—ki)"" 
Ck ; 







v= 









Substituting 


ps 





er 

p+Gi2(s) 

in (16) and inverting the order of integration we obtain Laplace integrals and 
we can write down the contribution of (16) to ¢. The complete expression for 
¢ is, from (15) and (16), 


°%aw 
—pt —G1,2(s)(t—s) 
- | goers 
$s 












1 (z—h) 
i) Qn * “Se 
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—f exp] -_ Sh) Nanak, )Re J F,(s)(s°—R})~ exp] 8 (s*—ki) 


— od ds 


8 exp] =¢ ‘en naa, ) Re J F3(s)(s°—Rj)~  exp| & (s° *_Ri)t 


™ (-s)Ga(s) | ds 


where Fi(s), F2(s), Gi(s), Go(s) are given in (17). 

The first two terms of (18) are of the nature of a source and image. The 
first term is, in fact, the elementary solution of equation (2) (Van der Pol and 
Bremmer 1955). The second and third terms represent the effect of the surface 
of the fluid. The wave aspect of = solution is apparent since all terms are 
multiplied by H(t—k,) or H(t—k2). Thus nothing is noticed at a point of 
observation in the fluid until a pia interval of time has elapsed. 

The case of the source and point of observation situated on the surface of the 
fluid is obtained by letting h — 0 and z — 0 in (18). A difficulty arises in that 
the two integrals that result are divergent at the lower limit. However, if the 
operation of taking the real part is performed and simplifications made, 
convergent integrals are found. We avoid these difficulties by writing the 
lower limit as (x/c)+e, and ultimately letting « + 0. Letting h and z tend to 
zero, (18) reduces to 


Eni) Ref’ rio (st-%) ex] ( 
ss - “—" es 4! (z/e)+e Fils) ApoE: \ 
— (366s) | 
~£n(1-*) Re { Fuls)(s*-%) ex E Lar 
v c Sista | N e sins 2c Cc 


~ (o-s)Gu(s) | ds 


. 2\4 - 
Gi,2(s) < s “) +(s-*)], 


: i[(2_x) x__sx_)| 
F;,2(s) Sel. a a(s+2 #)|. 


If we let c — © in (19) we obtain the solution for an incompressible fluid. 


This gives 


ai et) 
(20) ay sn) - ds 


which by a change in variable reduces to the expression given by Stoker 
(Stoker 1957a). 
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ASYMPTOTIC EVALUATION 
In order to obtain a clearer understanding of (19) we shall find an asymptotic 
estimate using the method of steepest descent (Erdélyi 1956a). 
Substituting z = s/t (the variable z not to be confused with the distance 
co-ordinate), w? = gl?/2x, and a = x/ct (0 < a < 1), (19) becomes 


(21) o= e- H(1—a) Im I 


“= 





a) Im I» 


where 


a1 iw? g(z) 
a ee 


+< (2—a) (2 —a a 


i J ee way (z)dz 


ate @—a) (ar) 

g(2) = (2°—a”)(e+a—1)+(2—a)(z—1), 
gi(z) = (2°—a*)*(s+ta—1)—(s—a)(z—1), 
f(z) = (2-a) (2*—a*)*+ (2’—a*—az), 
fi(z) 


If we regard w* as large, the integrals 7; and J; in (21) are of the form to 
which the method of steepest descent can be applied to obtain the asymptotic 
representation. We deform the path of integration into the steepest path from 
a to 1 in the complex z-plane. We shall set up a two-sheeted Riemann surface 
with branch points at g = +a and branch cut the real axis joining s = +a. 
The path of integration for 7; is on the upper sheet of the Riemann surface. 
We see that 7; can be written as 


ial a ‘*?F(a)ds 


ate ae —a@ a 


Io = 


(z—a) (2?—a’)— (2° —a’°—az). 





where the path of integration is now on the lower sheet of the Riemann surface. 

The saddle points are the values of z for which g’(z) = 0, and there are three 
for all values of a. Call them 2:(a@), 22(a@), 23(a). If 0 < a < 4, 20(@) and 23(a) 
are on the lower sheet and z;(a@) on the upper sheet, and they are all real. 
As a increases through 4, z2(@) moves up to the upper sheet. When a = 43/3 
= (0.5291 the two upper sheet saddle points coincide. For @ in the range 
B/3<a < 1, z:1(a) and 2:(a@) are complex conjugates. 

In Figs. 2 and 3, the deformed paths of integration for J; and J; for a ranging 
from 0 to 3, and from 0.53 to 1 are shown. These deformed paths are the 
steepest curves, that is the curves Re g(z) = k, going from a to 1. In forming 
the asymptotic estimate for J; there will be three parts: (a) the integral from 
ato ~, (b) the integral over the saddle point, and (c) the nee from © tol. 
In ee ase of J. these three parts are present only for0 <a <3.Ifi <a<l 
there is no integral over the saddle point. 
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Fic. 2. The curves Re g(z) = k on the upper and lower sheets, for 0 < a <}. Curves on 


the upper sheet are solid, and on the lower sheet dotted. The branch cut extends from —a to a. 
The saddle points are at x:(a) and x2(@). The deformed path for J; joins in succession the 


points a, A, x:, B, 1. The deformed path for J2 joins a, C, D, x2, E, 1. 


We shall show that the two integrals from a to ~ cancel. For J; the steepest 
curve from a to © will be g(z) = 7s and for J, it will be g(z) = —7s, where 
s is a real parameter ranging from 0 to ~. Thus J; becomes 


J e*'F(s)ds 
and J, becomes 
j eS F*(s)ds 


where F(s) is a complex function of the real variable s and F*(s) is its complex 
conjugate. The contribution to ¢ will be 


-h. H(1—a) im f €°"'G(s)ds, 
TCA € 


where G(s) = Re F(s), and this is equal to zero. 
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A 


ee 
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Fic. 3. The curves Re g(z) = & on the upper and lower sheets, for 0.53 <a@ < 1. Curves on 
the upper sheet are solid, and on the lower sheet dotted. The saddle point is at 2}. The deformed 


path for J; joins in succession the points a, A, B, C, 21, D, 1. The deformed path for J2 joins 
a, E, F, 1. 


The integrals from © to 1 may be treated as follows. 

For J; the appropriate steepest path is g(z) = av/(1 — a?)—is, and for 
I, it is g(z) = —av/(1 -- a®?)—is, where s is a real parameter ranging from 
0 to ~, with s = 0 corresponding to z = 1. J; becomes 


(22) — eiwtavil oy oF (s)ds 
0 
where 
F\(s) = f(x).dx/ds 
1S) = re 


(x—a) (x’—a’)'’ 
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with x a function of s determined from g(x) = avV/(l — a*)+is, and I: 
becomes 


(23) gto fl e° ' Fo(s)ds 
0 


where 
. (x). dx /ds 
Fels) = 
2(s) (x—a@)(x°—a*)'”*’ 

with x a function of s determined from g(x) = —av/(l — a*)—is. The 
asymptotic evaluation of (22) and (23) is found by expanding F;(s) and F2(s) 
in a power series in s, and then integrating term by term. 

The imaginary part of the sum of the leading terms in the asymptotic 


expansions of (22) and (23) is 


(24) 


To evaluate the contributions from the saddle points, the various ranges of 
a must be handled separately. 
I. The Case of 0 ca <4 

In this case x;(a@) contributes to J; and x2(@) to J. By applying Laplace’s 
formula (Erdélyi 19560) we find the leading term in the asymptotic contribu- 
tion to ¢ to be 


(26) aninaasllaesin hy(a) sin (w"pi(a)+ } 2 _ ho(a) sin(w pola = 1) 
27Caw 2 1rcaw 


where 


(x—a) Vx" -a' V/g x) ; 
Pi,2(a) 2(xX1,2(a@)). 


In Figs. 4, 5, 6, 7, ti(@), he(a), pila), p2(@) are plotted as functions of a. 


x = X1,2(a), 


hy2(a) = 


II. The Case of a = 3 
We see that the application of Laplace’s formula to the evaluation of J, 
when a is close to 3 is invalid, since lim h2(a) = ©. Under these circumstances 
a—1/2 
it may be shown, by using the method of stationary phase (Erdélyi 1956c), 


that the contribution from x2(a@) to ¢ is 
ee 953 g 
(27) 39 ee 
ae V reaw 
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a 


FiG.4. p,(a) plotted asa function of afor0 <a <4} 
Fic. 5. p2(a) plotted as a function of a for 0 <a < 3. 
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3 *\ 2 3 “4 5 
a 
Fic. 6. /,(«) plotted as a function of a for 0 < a < 0.53. Ii(a) — — © asa > 44/3 = 0.53. 


III. The Case of } < a < 44/3 

In this case both x;(a@) and x2(a) contribute to J;, and there is no saddle 
point part of J». The contribution to ¢ will be the same as (26). 
IV. The Case of a = 44/3 

Since x;(a) and x2(a) coincide when a = 4!/3, Laplace’s formula again 
breaks down. However, using an extension of Laplace’s result (Erdélyi 19560) 


> 


we find for the saddle point contribution to ¢ in the vicinity of a = 4'/3 


48.62 = fe 
~ oe “an sin(w pi(a)). 


(28) : 
=#TCQA W 
V. The Case of 44/3 <a <1 
In this range of a, x:(@) and x(a) are complex conjugates. Only x;(q@) 
contributes to J;, while J, has no saddle point part. The saddle point contribu- 
tion is 


w? g(a) 


(29) hy(a) sin(w pi (a) +0(a))e 


g 
V2 reaw 
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9 “| “2 “3 “4 ‘5 
a 


FiG. 7. he(a) plotted as a function of a for 0 <a < }. Inla) ~ ~ asa >}. 
hy(a), pila), Q(a), and g(a) are shown graphically in Figs. 8, 9, 10, and 11. 
The expression (29) is, of course, small compared with the contribution from 
a= i. 


Thus we have, in summary, 


Ia —-—5 4 9 9 ———5 9 
(30) ¢= ey Vi —a’(2a*°— ta“ —a+4) ‘cos(warv/1 —a’)+ S(w", a) 
TCaW 
where 
S(w', «) = & hy(a) sin(w pi(a)+i2) 
2 reaw 


— —2 — fy(a) sin(w*p2(a)—}4) (0<a < 4!/3) 


2 mcaw 
except near a = 3 and a = 43/3. 
lia:=%, 
Sj - ~<a be a date —he: 
V2 caw 32 V reaw" 
Ifa = 43/3, 





Shel . 4%; 
———8 35 sin(u'p:(a)). 


“7.2 
S(w',a) = 3 
“7CQM W 








ROSS: EXPLOSION IN COMPRESSIBLE FLUID 


S(w*, a) = —&— hy(a) sin(w' pi(a)+0(a))e "(44/8 <a < 1). 


V. 27Caw 


G *7T a ~_ 3 


Fic. 8. h(a) plotted as a function of @ for 0.53 <a <1. 
Fic. 9. pile) plotted as a function of aw for} <a <1. 
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“6 “7 8 3 Oo 
a 


Fic. 10. 6(a) plotted as a function of a for} <a <1. O(a) > —3r/4as a — 3. 
Fic. 11. The damping factor g(a) plotted as a function of a for 0.58 <a <1. g(a) 0 
as a > 0.53. 










ROSS: EXPLOSION IN COMPRESSIBLE FLUID 





THE ELEVATION OF THE FLUID SURFACE 
From (4) and (30) we have approximately (omitting all except the leading 
term) 


a1) ¢ = 4/24 in(@px(a) cos(a'ps(a)+4x)— 4/22 hala) pala) 


cos (w"po(a) —in 










and this expression shows that the surface elevation is large and increasing 
with time. 







COMPARISON WITH THE INCOMPRESSIBLE CASE 
The expression for the velocity potential for the incompressible case given 
in (20) can be written as 








t BP witli: 
(32) e== Im f gerne. 
TX 0 
Applying the method of steepest descent to (32) gives 
Bah g ; “ =) 
9 a 2 
and ~~ ——— 2) 
c ae +Te$ * cos( 4 4)? 
the first term being the sum of the contributions from z = 0 and z = 1, and 






the second term arising from the saddle point at z = }. 


We would expect to get this same result by letting ¢c—> © (i.e. a—> 0) in 
(30). We find from (30) 


fa A ee n( wo 2) 
lim ¢ = Li g/Lsin aT 4 : 


That this is not the same as (33) is the result of a lack of uniformity in the 
asymptotic representation (30), with respect to the variable a. Using expression 
(26), which is derived from Laplace’s formula, the saddle point part of J, 
apparently tends to zero, if we let a — 0. In actual fact the saddle point part 
of I, tends to —1/zt, a result which now checks with (33). Thus (26) is invalid 
near a = 0. 

It is seen that the compressibility of the fluid produces the following results: 

(a) If a time ¢ has elapsed after the explosion, the disturbance on the fluid 
surface has reached no further than distance x = ct from the explosion. 

(b) Ifx/ce < t < x/0.53c the waves on the surface are of small but increasing 
amplitude as ¢ increases. These waves are of a nature similar to the precursor 
waves arising in the propagation of an electromagnetic disturbance through a 
dispersive medium (Stratton 1941). 

(c) When t = x/0.53c there is a sudden rise in amplitude as the main 
disturbance arrives. 

(d) If x/0.538¢ <t < © there are two distinct sets of large amplitude 
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surface waves. The set arising from x;(@) has an amplitude indefinitely increas- 
ing with ¢, while the other set becomes insignificant as ¢ tends to infinity. 

(e) As t tends to infinity the surface waves become the same as the waves 
predicted by the incompressible theory. 


The author is grateful to the Canadian Mathematical Congress for a summer 
research grant. 
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APPLICATION OF FOURIER TRANSFORMS TO 
ANTENNA PATTERN SYNTHESIS! 


ALFRED T. VILLENEUVE? 


ABSTRACT 


The two-dimensional Fourier transform of the field quantities is employed to 
synthesize source distributions over plane rectangular apertures when a field 
distribution is prescribed at an arbitrary distance from the source. Several 
examples of application of the technique are included. The approximations 
involved are discussed. 


I. INTRODUCTION 


When antennas of large aperture are used in radio transmission systems, 
the receiving antenna may not be located in the Fraunhofer region of the 
transmitter, or vice versa. Therefore, in certain antenna applications it may 
be necessary to obtain a prescribed variation of field intensity over a region 
which is not in the antenna far-zone. Under these conditions the known 
methods of far-field synthesis may not be applied directly and new techniques 
must be found. 

If the antenna in question is a plane rectangular aperture and the field 
variation is prescribed over a plane parallel to this aperture, then the two- 
dimensional Fourier transform may be employed to synthesize the required 
aperture distribution. Therefore, let us first determine the relationships which 
exist among the Fourier transforms of the various field quantities. 


II. FOURIER TRANSFORM APPROACH 
An arbitrary harmonically time-varying electromagnetic field may be 
derived from two parallel vector functions whose direction in space is fixed. 
Let these functions be denoted by 


(1) A =t,"; F =i 


where y” and y¥% are solutions of the scalar wave equation. In rectangular 
co-ordinates general solutions may be obtained by linear superposition of 
plane wave solutions. Such a superposition may be expressed for radiation 


problems as 


7 I a. p ju jr 
(2a) y= Qn): j f(uv)e7"e’"e"? du dv, 
ov Je 


(2b) y= Gey | j g(uvje re’ e"”* du dv, 


1Manuscript received March 16, 1961. 
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Hughes Aircraft Co., Culver City, California. The research reported in this paper has been 
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contract No. 500-7119 and was performed at Syracuse University. 
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where 


9 


(2c) w= BP — 4 — 2 
and 
(2d) B= wV/ (ue). 


The electric and magnetic fields may then be obtained from these by the 
following relations: 


(3a) E= + uxvxd-vxF, 
Jwe 


(3b) A = vxd4+—- vxvxe. 
jou 


Performance of the above operations gives the following expressions for the 
rectangular components of # and H: 


2 2 as 
(4a) £E,(x,9,2) = =a = f(u,v)e"e’e"”* du dv, 


(4b) E,(x,y,2) = a fff-s =f (u,v) — jw g(u |e Wz pi oie? Jy dy, 


2 (x.y.2) = jux joy ,jwz . 
(4c) E,(x,y,2) a JIL-E = J(u) 4430 (ua) |e ee” du dv, 


1 
(2xr)* jou 


(4d) H,(x,y,2) = Om vere %e"”* du dv, 


4 


(4e) H,(x,y,2) = oa JJ | -#* g(u,v) iu fue) ee? ei”? du dv, 


(4f) H.(x,y,z) = ay J | = g(u,v) jo flu | eee? du dv. 


If we choose our paths of integration along the real axis from — © to + 
we have the Fourier integral inversion theorem (Sneddon 1951). Thus, if we 
denote the Fourier transform of E;(x,y,z) by & ;(u,v,z) we have the transform 
pair 


(5a) C,(u,v,2) = f J E,(x,y,2)e ““e?™ dx dy, 


(5b) E,(x,y,2) = ars J E,(u,v,2)e™"e'™ du dv 


with a similar pair for #;,and H;. From the uniqueness of the Fourier transform 
pair we then have 


w 49" 


(6a) C6 ,(u,v,2) = — Fue SM ver” . 
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(6c) € ,(u,v,2) = -| #2 7tu —ju e(us) on, 


2452 
WwW v jwz 
3 uve *, 
jap 8°) 


(6e) H,(u,v,2) = -|2 g (u,v) jw fue) | en, 


(6d) KH, (u,v,2) 


(6f) HK (uv,2z) = _[= g(u,v) +ivflus) | on. 


Two of the field quantities may be chosen independently and the rest are 
automatically determined. In what follows we shall make FE, and E, the 
independent choices. Let E, and EF, be specified over a plane z = 2;. From 
the foregoing equations we have for the transforms of the field quantities for 


any value of z: 


(7a) E (uy,z) = &,(uv,ziye"*™, 


(76) E,(u,v,2) = €,(u,v,2i)e"?*, 


ai a u = v ai 
(7c) & ,(u,v,2) _ fs O',(u,v,21) + E (us) | eiv’? me, 


j 1 2 2 2, @ ss 
H,(u,v,2) ouw & (u,v,21) + (B°—u")&, (u,v,2) le ™, 


1 


HF hana) = = [(8°—v*)&,(u,v,21) +u0e, (u,v,21) Je, 


WP ‘ 1 2 : _ 
HK ,(u,v,3) = ia [v&,(u,v,21) —u&,(u,v,21) Je" oem 


If the sources of the field are located at z = 0, then for values of both z and 
2, > 0, the above expressions must represent waves travelling or attenuating 
in the positive z direction. We must therefore choose the branch of w which 
gives this characteristic, i.e. for real uw and v 

(8a) jw = — jv (8 -- v — v*) a +o < #, 
(8d) ju=—-V(w+v — B) u® -+-o* > fF. 


For both z and 2; < 0, w in the field equations is replaced by —w. 


III. RELATION TO SOURCES AT z =0 


If the sources of the field are electric and magnetic current sheets in the 
x—y plane the fields are simply related to them. The boundary conditions 


at 2 = O are 


(9a) J,(x,y) = a X[A(x,y,0+) —A(x,¥,0-)], 
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(90) M,(x,y) a —u, X(E(x,y,0+) — E(x,y,0-)], 


where /, and M, are electric and magnetic surface current densities respect- 
ively. The transform quantities also satisfy the same boundary conditions. 
In component form these are 


(10a) JZ (tv) = H,(u,v,0+)-—H#’,(u,v,0-), 
(10d) A (Uy) = —H,(u,v,0*)+ H,(u,v,0-), 
(10c) My (uv) = —&,(u,v,0+)+&,(u,v,0-), 
(10d) Me(uv) = &,(u,v,0t+) —&,(u,v,0-). 


In terms of &, and &, the first two boundary conditions become 


(a) Suy(ue) = = {uolé (u,0,0*) +62(u9,07)] 


+ (8°— u’) [e, (u 0,0") + &,(u,v,0 )] i 
(11b)  Ser(u,v) { (8°—v°)[E,(u,v,0*)+6,(u,v,0" )] 


+uv[o ,(u,v,0t)+&,(u,v,0-)]}. 


If at z; > 0 and at z. < 0 the field components £, and E, are prescribed the 
transformed boundary conditions then become 


(12a) Wy (uv) = —&,(u,z1)e +6, (u,v,22)e", 


jw22 


(12b) MW ,(u,v) = ©, (up,z)e™" —€ , (u,v,22)e", 


(12c) JS w(U,v) = {uv[O,(u,v,2:)e 2” +6 ,(u,v,22)e""] 


+ (8°—u*)[E, (up, ze "+6, (u,v,22)e"" ]}, 
(12d) Jar (u,v) = {(8°—v°)(E,(uv,2z:)e "+6, (u,v,22)e"*] 


+uv[O,(u,v,2:)e 7 +6 ,(u,v,22)e""]}. 


From equations 10 and 11 we may examine several cases of interest: 

(a) If a,XE is discontinuous across z = 0 both electric and magnetic 
source currents are required in general. 

(b) If a@.XE is an even function of z there will be only electric current 
sources required. 

(c) If a@,XE is an odd function of z then only magnetic current sources 
are required. 
Case (a) is the general case. If the required sources are of finite extent this 
case is mathematically equivalent to the problem of diffraction by an aperture 
in a perfectly absorbing screen. Also, if the sources are of finite extent case (b) 
is mathematically equivalent to radiation from an aperture in an infinite 
plane of infinite impedance or from a sheet of electric currents in free space. 
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Under the same space limitations of the sources case (c) is mathematically 
equivalent to radiation from an aperture in an infinite perfectly conducting 
plane or from a sheet of magnetic current* radiating in free space. 


IV. SYNTHESIS OF SOURCE CURRENTS 

From equations 12 we have the transform of the source distributions. The 
evaluation of the inverse transforms, however, is complicated by the appear- 
ance of w in the expressions. Since this has branch points the resulting integrals 
cannot be easily evaluated. If the sources are confined to a rectangular aperture, 
however, it is possible to evaluate them as complex Fourier series. 

The procedure for synthesizing the source distributions is as follows. In 
this discussion let S; and .“7; represent a source function and its transform, 
respectively, and let A; represent the right-hand sides of equations 12. The 
situation is illustrated in Fig. 1. If the sources vanish everywhere outside 


Fic. 1. Orientation of rectangular aperture. 


the rectangular aperture of dimensions L, and Ly they may be expanded in 
double complex Fourier series as follows: 


Bi —. — ,J(2mx/Lz)z j(2ne/Ly)y |x| < L,/2, 
(13) Silesy) = 2 DS Cone é ly| < L,/2, 


Si(x,yv) = 0 all other x,y. 

The Fourier transform of S;(x,y) is then given by 

a = a «SI L,/2)—mr] sin[(vL,/2)—nx] 
14) L(uv) = LoL , inl(uls/2)—ma] sin[(oL,/2)—n) 
( 4) (a e oe 2, [(uLy ‘2)—mr] [(vL,/2)—nx] 
From equations 12 then we have 

- . sin[(uL,/2)—m-a] sin[(vl,/2)—nx] 7, | 

(15) Lely 2D Con [(uL,/2)—mr] [(vL,/2)—nx] whee) 


*Physically, a magnetic current element may be replaced by an equivalent electric current 
loop. 
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If the values u = 2mxr/L,, v = 2nm/L, are substituted into equation 15 we 
get 





(16) LeLyCun = {2m 7 ae) 
and the Fourier coefficients C,,, are completely determined. Thus, the sources 
are synthesized. 

Equation 16 points out the existence of certain restrictions on the Y;, and 
therefore on the prescribed field, for fixed values of L,; and Ly, if one is to avoid 
supergaining with its attendant high Q and high loss conditions (Harrington 
1958, 1959, 1960). Conversely, for given Y, the minimum dimensions of the 
aperture are fixed, for a non-supergain aperture distribution. This results 
from the fact that C,,, corresponding to values of m and n lying outside the 
ellipse 


9 2 


m nN 
(L,/a)?* (Ly/ay? 


result in primarily non-propagating fields which contribute principally to 
energy storage. Therefore, if L, or L, is too small, the points given by 


(17) 1 


27 2n7 


’ Un = 
i. 


may be too widely separated and the series on the left of equation 15 may not 
converge to .¥; between these points without the inclusion of these primarily 
non-propagating modes. To insure the convergence it becomes necessary that 
some Cy, which give primarily evanescent fields be added to the aperture 
distribution or that the aperture dimensions be increased to bring the wu, and 
v, Closer together, thereby giving greater control through increasing the number 








be = 






of propagating fields. 

For a given size of aperture this gives a method of checking the assumed 
field distributions to determine whether they can indeed be obtained from the 
aperture with a low Q distribution. If the series of equation 15 does not con- 
verge to the .A, between the sample points then, for a non-supergain distribu- 








tion, a larger aperture is required or the prescribed fields should be modified. 






Unfortunately there seems to be no simple a priori way to determine whether 
a prescribed field is obtainable from apertures of fixed dimensions. The type 
of approximation to be expected when the aperture is not large enough to 
produce the prescribed field is discussed in the Appendix. 








V. APPLICATION OF SYNTHESIS TECHNIQUE 






The synthesis technique has been applied to the case (c) previously 





mentioned, for several prescribed fields. In this case the sources are given by 


(18) M, = —a,X2E(x,y,0*). 







The prescribed fields have no y component of E. They were prescribed over 
a plane a distance 10D from a square aperture D on a side, with D = 100 X. 
The prescribed fields are as follows: 
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Case No J 
z/D«= /0 y/0 =O 
— Compureo Fieio 


—— Normaiizeo PrescriaeD FitLd 


Ql 02 03 04 as . 4 0.8 ag 
x/D —— 


Fic. 2a. Comparison of prescribed and computed fields. 


: 


Ss 


w 
y 
q 
9 
g 
% 


§ 


04 05 
 —— 


Fic. 2b. Synthesized aperture distribution for case I. 
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Case No I 
2/D =/0 y/D=0 
— (Compureo Fiéeto 


—— Normazeo Peesceiszo Fittp 





Dé Ex | ee 




















ny 
§ 


y 
¥ 
q 
Q 
z 
% 








AMPLITUDE 




















Os Ob 
x/D —— 


Fic. 3b. Synthesized aperture distribution for case II. 
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L2 
Case No I 

ig 2/D=1/0 ~—-y/D=Q25 

— Comeureo Fieto 
—— Noemaizen PeesceiBen feo 

a 
uy a8 
3. 
Q 














oO 0./ 02 a3 QO4 Qs 06 Q7 ag ag 40 
2/0 —— 


Fic. 4. Comparison of prescribed and computed fields. 
(1) Gaussian distribution: The normalized prescribed curve is 


(19) D’E,(x,y,10D) = eet =O 


(II) Superposition of Gaussian distributions: The normalized prescribed 


curve is 


¢ 2m (uy = (4.01)° —(4.01)2(z/D—.195)? —(4.01)2(z/D+.195)2 
(20) D*E,(x,y,10D) = 1.08 “~~~ [e +e | 


—(4.01)2(y/D—.195)2 —(4.01)2(y/D+.195)? 
X |e P +e ’ ]. 


The resulting fields computed from sources synthesized by the above technique 
are compared with the prescribed fields in Figs. 2a, 3a, and 4. Figures 2a and 
3a are evaluated for y = 0, Fig. 4 for y = .25D. The corresponding sources 
are shown in Figs. 26 and 3b. They are also evaluated for y = 0. The variations 
in the computed patterns indicate that the aperture must be somewhat larger 
in order to reproduce the prescribed patterns without introducing some 


degree of supergaining. 
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APPENDIX 


If, when only non-supergain aperture distributions are considered, the 
required sources of a prescribed field variation extend beyond the available 
physical aperture, then the field that results from approximate sources deter- 
mined by the transform synthesis method will only approximate the specified 
field. One may wish to know what type of approximation is obtained by this 
method. It has been shown (Harrington e¢ al. 1958) that in general it will not 
be an approximation in the least-squares sense. However, for non-supergain 
aperture distributions, the least-squares approximation is obtained if the 
approximate sources coincide with the exact sources over the physical aperture. 
Since in general these are unknown, this cannot be directly accomplished. It is 
possible, however, to approach this latter technique by mathematically 
extending the aperture until it is sufficiently large to encompass essentially 
all the necessary sources. The coefficients of the Fourier series over the extended 
aperture may then be determined. These give the correct sources which may 
then be chopped off to fit into the physical aperture. If this resulting distribu- 


tion is non-supergain the pattern will approximate the prescribed one in the 
least-squares sense. 









THE APPLICATION OF A FLUXGATE MAGNETOMETER TO AN 
AUTOMATIC ELECTRONIC DEGAUSSING SYSTEM! 






R. L. GRAHAM AND J. S. GEIGER 






ABSTRACT 


The three-component fluxgate magnetometer developed by Serson has been 
adapted to provide continuous correction signals to the degaussing system of 
the Chalk River iron-free 8-ray spectrometer. Improved electronic circuitry has 
been developed for the magnetometer which minimizes the zero error and 
reduces to <10~4 gauss long-term drift caused by component aging. The degaus- 
sing coil arrangement used to generate uniform magnetic field components 
opposite to those of the earth is indicated and the current regulation system is 
described briefly. The method in which the magnetometer correction signal is 
introduced into the current regulators is shown and an example is given of the 
performance of this degaussing system during a magnetic storm. 















I. INTRODUCTION 


The Chalk River 1-meter radius, iron-free, double-focusing 2/2 8-ray 
spectrometer (resolution 20.01%) (Graham, Ewan, and Geiger 1960) was 
designed on the assumption that the earth’s and other stray magnetic fields 
could be cancelled completely in the region of the focused electron orbits. 
On this assumption, the only magnetic field which influences the electron 
trajectories is that created by the coil system of the spectrometer (Lee-Whiting 
and Taylor 1956, 1957); hence the magnetic field strength and the momentum 
of the focused electrons are both strictly proportional to the current flowing 
in the spectrometer coils. This spectrometer current is regulated to, and can 
be measured with, a precision of a few parts in 10°. The spectrometer is located 
in an iron-free building whose location was chosen to avoid any local distortion 
of the shape of the earth’s magnetic field. Two sets of large degaussing coils 
are used to create magnetic field components (vertical and north-south) 
which are equal and opposite to those of the earth’s magnetic field in the 
region of the electron orbits of the spectrometer. In initial operation the 
currents in these degaussing coils were held constant, the values being set 
empirically to give zero field at the spectrometer when it was turned off, this 
zero field condition being established with a fluxgate magnetometer detector. 
Any subsequent variation in the earth’s magnetic field produced an uncertainty 
in the spectrometer calibration. When electrons of energy ~80 kev (momen- 
tum ~1000 gauss cm) are being focused, the spectrometer magnetic field 
strength on the optic circle (radius 100 cm) is ~10 gauss. Thus, if the earth’s 
field changed by, say, 10-* gauss the ‘‘apparent’’ momentum of the focused 
electrons differed from the ‘‘real’” momentum by 1 part in 10°. Diurnal varia- 
tions of the earth's magnetic field are typically 510-4 gauss and disturbances 
due to magnetic storms following solar flares are sometimes as large as 10-? 
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gauss. In order to eliminate the calibration uncertainty resulting from these 
fluctuations we have recently installed an automatic compensation system 
which maintains the cancellation to <10~* gauss on a long-term basis. 

Our original intention was to replace the constant current regulators by a 
current regulation system which received direct feedback from a three-com- 
ponent fluxgate magnetometer placed at the center of the spectrometer coil 
system. The coil system of the spectrometer had been specially designed to 
have “‘null’’ field at its center of symmetry (Lee-Whiting, ref. #21 in Graham 
et al. 1960). An unexpected difficulty was encountered, however, which made 
this direct method of control unacceptable. The “null” field condition at the 
center of the spectrometer can be upset by several factors, one of which is the 
generation of eddy currents in the cooling channels and coil support structure 
that results from any change of the spectrometer coil current. These eddy 
currents give rise to transient magnetic fields at the center of the spectrometer 
coil system which are picked up by the magnetometer and give a false indica- 
tion of a degaussing error. Another difficulty is that the cancellation of the 
spectrometer fields in the central null region is quite sensitive to ambient 
room temperature whereas the magnetic field in the electron orbit region is 
independent of room temperature (Graham et al. 1960). 

An alternative approach to this automatic control problem was therefore 
adopted which assumes as its basis that the variations in the components 
of the earth’s magnetic field at the spectrometer are proportional to those 
observed at a location some 400 feet from the spectrometer building. The 
constant current regulators used during the initial operating period have been 
retained with the minor addition of voltage adding networks which permit the 
standard cell reference voltages of these regulator units to be varied in pro- 
portion to the changes in the earth’s freld observed at the remote magnetometer 
location. The remotely located magnetometers are used to measure the 
magnetic field strength in this arrangement and must therefore have long-term 
gain stability in addition to a good stability of their zero setting: an over-all 
long-term stability of <10~-* gauss is required. 

We found that the long-term stability of the fluxgate magnetometer 
developed by Serson (1957) was not adequate for this application. Tests showed 
that the stability of the magnetometer was related to the magnitude of its 
zero error, i.e., to the field strength indicated by the device when no field was 
present. In this paper we discuss some of the factors which contribute to the 
zero reading error of this instrument and present a modified version of the 
control circuitry for this instrument which experience shows to have a long- 
term stability of <10-‘ gauss. The use of this magnetometer for the automatic 
compensation of variations in the earth’s magnetic field strength is briefly 
described and an example is given of the performance of the complete 
degaussing system during a magnetic storm. 


Il. MAGNETOMETER AND ITS CONTROL CIRCUITRY 
The reader is referred to the earlier papers of Serson and Hannaford (1956) 


and Serson (1957) for a detailed description of the design and operating 
principle of the fluxgate magnetometer which we have adapted to the problem 
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outlined above. In brief, this fluxgate consists of two identical mu-metal 
strips on which are wound the primary coils (see Fig. 1, upper right). The 
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Fic. 1. This figure illustrates schematically the construction of the fluxgate magnetometer 
as developed by Serson and the associated electronic circuit components. The detail of the 
primary winding drive connections to the fluxgate head, of the frequency doubler circuit, 
and of the power amplifier circuits are given in Fig. 2. The circuits of the a-c. amplifier, the 
phase sensitive detector, and the d-c. amplifier are given in Fig. 3. 








sinusoidal exciting current in the primary coils is such that the mu-metal is 
‘saturated’ for part of each half cycle. If no steady magnetic field is present 
the saturation field pattern is symmetrical about zero and has only odd 
harmonic components. The presence of a steady field component in the 
direction of the strips introduces asymmetry causing even harmonic compo- 
nents to appear. The phases of the latter depend upon the sign of the steady 
field. The secondary winding is tuned to the second harmonic frequency. This 
signal is amplified and compared with a second harmonic reference potential 
in a phase sensitive detector, Fig. 1. The resulting d-c. signal is further ampli- 
fied and fed back as a direct current to the secondary solenoid which cancels 
out the steady field component. The fluxgate unit itself is used as a null 
device in which the steady magnetic field along the mu-metal strips is kept 
close to zero. 

In our early tests we constructed a number of fluxgate units in this laboratory 
according to Serson’s design (1956). Subsequently we purchased several 
commercially made units from Canadian Applied Research Limited (Toronto). 
The performance of the first commercial units we acquired, although superior 
to units of our own construction, did not match the performance of units 
constructed by this firm at an earlier date for the IGY program. We are 
indebted to Mr. R. L. Adams and Mr. E. B. Moss, of this firm, for re-examin- 
ing, in detail, their manufacturing techniques and for subsequently supplying 
us with fluxgate units of a high standard. 
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We tested the magnetometer circuitry given by Serson (1957) and a number 
of variations of it in our attempts to achieve the freedom from drift and the 
low zero error required for the present application. The basic functions of 
these circuits are shown in block diagram form in Fig. 1. The circuit character- 
istics which critically influence the zero stability of the magnetometer are the 
following. The drive current fed to the primary winding of the magnetometer 
must be entirely free of any second harmonic component. The presence of 
even a small second harmonic content will give a false indication of a steady 
field component even though the “push—pull’’ primary coils in the fluxgate 
unit are designed to minimize inductive coupling between primary and 
secondary coils. Since the secondary circuit usually detects a small amount 
of the fundamental frequency component, the second harmonic reference 
signal driving the phase sensitive detector must likewise be free of fundamental 
component or the phase sensitive detector will register such a component as 
a false signal as well. In the original three-channel magnetometer the three 
phase-sensitive detectors are supplied by the same second harmonic reference 
supply which has a rather high output impedance. We found that there is 
appreciable ‘cross talk’ between the channels because of this, i.e., a signal in 
one channel can cause a small ‘false’ signal to appear in the other two channels. 
In the original design there was no provision for zero adjustment of the d-c. 
amplifier. Hence a zero output signal did not necessarily correspond exactly 
to a zero field component at the fluxgate unit, the offset being a function of 
the d-c. amplifier tube characteristics and loop gain of the system. 

The circuits now used with the Serson fluxgate unit are shown in Figs. 2 
and 3. They are based on Serson’s (1957) circuits but have been modified to 
overcome the limitations discussed above. A commercial power oscillator of 
high quality provides the fundamental voltage signal of 660 c/s. The manu- 
facturer, Elin/Electronics International Co. of Burbank, California, guarantees 
a pure sinusoidal waveform with <0.2% harmonic distortion, a frequency 
stability of <0.1%, and a long-term amplitude stability of <0.1%. The 
revised frequency doubler circuit of Fig. 2 produces a purer 1320 c/s wave- 
form and has a lower output impedance than did the original circuitry. It 
includes provision for eliminating any fundamental (660 c/s) component, i.e., 
by adjusting P2. A higher gain a-c. amplifier is used to amplify the secondary 
coil signal (Fig. 3). The phase sensitive detector has been revised and works 
at a higher voltage level on a current adding, rather than a voltage adding, 
principle. This minimizes cross talk between the three channels via their 
common source of second harmonic reference voltage. Adjustment of P;, with 
S» closed, permits one to set the d-c. bias level for V; so that there is no inherent 
zero error, as indicated by M1, arising from d-c. unbalance at the phase 
sensitive detector. C; acts as a ripple filter to minimize direct feedback of any 
small possible 1320 c/s component in the d-c. output signal which could give 


rise to a zero error. 





Fic. 2. Excitation current control and frequency doubler circuits for single- and three- 
component magnetometers. The 600 c/s supply is an ELIN power oscillator type DK-115 
(Electronics International Co., Burbank, California) which supplies both the single- and 
three-component magnetometers. 
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In the course of our experimenting with fluxgates, it became apparent that 
the mu-metal can ‘remember’ temporary exposure to a steady field of 210 
gauss even though it is repeatedly being ‘saturated’ during each excitation 
half cycle. The nature of this difficulty was not appreciated at first. We 
explored the use of other magnetic materials as fluxgate elements, e.g., super- 
malloy, in unsuccessful attempts to achieve more reproducible behavior. 
Sokolowski (1959) has also noted that the performance of the magnetometer 
employed by the Uppsala 8-ray group to adjust their spectrometer degaussing 
system deteriorated after exposure to magnetic fields of a few gauss. We have 
added a simple feature to Fig. 3 which is reminiscent of the technique used by 
watchmakers to demagnetize watches. By closing S; and varying the setting 
of P; slowly from full voltage to zero, the mu-metal ‘memory’ is wiped out by 
the large (100 gauss 0) a-c. magnetic field generated by the secondary 
winding. In this way memory ‘zero errors’ of up to 5X10~* gauss or more can 
be reduced to <3X10-* gauss. The method used to establish the ‘zero error’ 
experimentally is that described by Serson and Hannaford (1956). If the 
fluxgate is rotated through 180°, the output signal should have the same 
magnitude but opposite sign. If the magnitudes are unequal, the ‘zero error’ 
is taken to be one-half the difference. 

A further difficulty which we encountered in using the fluxgate to measure 
precisely the vertical component of the earth’s field results from the use of the 
secondary coil of the fluxgate unit as the solenoid that generates the steady 
bias field. Our experience shows that the performance deteriorates noticeably 
for field components of 20.1 gauss. The indicated field then becomes dependent 
upon the magnitude of the excitation current. With the fluxgate unit oriented 
vertically in the earth’s field (B, = 0.57 gauss) the observed shift is ~150 
gamma (1.5X10-* gauss) for a 10% (~10 ma) change in excitation. This 
effect is apparently due to the non-uniformity of the bias field generated by 
the secondary solenoid and is more pronounced for some units than for others. 
The secondary winding length is only slightly greater than that of the 4 in. 
long mu-metal strips. Non-uniformity of the bias field at the ends of the 
solenoid is believed to be responsible for this effect. Since the vertical compo- 
nent of the earth’s field is of greatest concern to us we have chosen to circum- 
vent this difficulty by placing the vertical fluxgate unit inside a longer solenoid 
(2 in. I.D. X7 in. long) which generates a more uniform vertical bias field. This 
solenoid is wound on a pyrex glass tube in order to minimize variations due 
to thermal expansion. The current for it is derived from a 105-volt source 
which is stabilized to <1 part in 10° rather than the OG3 voltage reference 
tube used by Serson. The use of the larger solenoid greatly improves the 
performance at 0.57 gauss (vertical component at Chalk River) but does not 
quite match that observed in the absence of a steady field component, e.g., 
when oriented in the east-west direction. 

The performance of the C.A.R.L. commercial Serson fluxgate units in 
conjunction with the circuits of Figs. 2 and 3 now meets our requirements. The 
zero error, i.e., indicated minus actual field, in a field-free region, is 5S X10-5 
gauss after demagnetization. The indicated field reading is essentially inde- 
pendent of excitation current over a wide range, e.g., 70 to 120 ma with a 
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typical slope of <2X10- gauss per 10-ma excitation current. The exciting 
current is normally maintained at or near 90 ma which gives the optimum 
second harmonic signal from the tuned secondary winding and hence the 
highest loop gain (2100) of the complete system. The indicated field reading 
is insensitive to primary current phase which can be reversed by S;. Long-term 
drifts due to tube aging are also negligible; the most serious effect of this type 
can be expected to result from drift in the d-c. characteristics of Vs. While 
provision has been made for a manual correction of such a variation (Ps), the 
long-term zero drift has been found to be <10~ gauss per month without 
readjustment. 


III. COMPENSATING COIL ARRANGEMENT AND CURRENT REGULATORS 


The array of coils used to cancel the vertical and north-south components 
of the earth’s field is shown schematically in Fig. 4. The coil parameters 
listed in the caption are due to Lee-Whiting (1957). This array is designed to 


SPECTROMETER 
AXIS 


5 
FEET 


Fic. 4. Schematic diagram of coil array used to cancel the earth’s magnetic field in the 
region of the 6-ray spectrometer. The coil parameters are as follows: 
Diameter Separation Number of Current 
Coil No. (ft) (ft) turns (amp) 
z 14.54 9.46 22 | 
2 3 27.50 9.2: 66 | ~2. 85" 
5, 8 27.50 5.85 52 | 
6,7 27.50 5.6 23 | eae? 
The currents listed are for (a) a vertical field component of ~0.572 gauss and (6) a north- 
south component of ~0.13 gauss. 
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generate a compensating magnetic field which is uniform to <10-4 gauss 
within a 9 ft diameter sphere surrounding the center of symmetry. The 
reader is referred to our earlier publication (Graham et al. 1960) for a more 
detailed discussion of the magnetic properties. In addition to the two-com- 
ponent coil system shown in Fig. 4, a pair of square coils (26.4 ft26.4 ft 
separated by 14.39 feet) has been installed for the purpose of cancelling an 
east-west component. As yet these are not needed. 

Figure 5 gives a simplified schematic diagram of the current regulation 
system and coil arrangement used to cancel the north-south component of 
the earth’s magnetic field. The current regulators themselves are of con- 
ventional design. The loop gain is sufficient to maintain the degaussing coil 
current constant to within a few parts in 10° on a long-term basis. The standard 
cell and the resistance elements shown in Fig. 5 are located in a constant 
temperature oil bath (28.8+0.1° C) in order to minimize drift. Potentiometer 
P. (Fig. 5) is adjusted empirically to produce a zero field component at the 
spectrometer. This is measured using a single-channel magnetometer oriented 
in the appropriate direction. 

A three-component magnetometer measures the vertical, north-south, and 
east-west components of the earth’s magnetic field and provides the correction 
signals that alter the currents in the vertical and north-south control systems. 
The three fluxgate units are located in a small remote constant temperature 
(+2° C) enclosure. Its location, in the woods some 400 feet from the spectrom- 
eter building, was chosen to minimize stray magnetic signals from other 
laboratories and vehicular traffic as well as the magnetic field generated by the 
8-ray spectrometer. The vertical unit does, however, detect a small field of 
~5 X 10~ gauss from the spectrometer when it is focusing energetic (~3 Mev) 
B-rays. The magnetometer units are adjusted to give zero output signal by 
appropriate adjustment of the baseline controls, Py, in Fig. 3. Subsequent 
changes in the earth’s field result in non-zero (positive or negative) signals of 
1 volt for a change in field of 1 X10-* gauss. The potentiometer, P7, of Fig. 5, 
is adjusted empirically so that a potential of V volts applied between S and G 
results in a change of VX10-* gauss at the center of the spectrometer. The 
resistance ratio shown in Fig. 5 (1.5 M to 11 K+0.5 K) is appropriate for the 
north-south component of ~0.13 gauss at the spectrometer. A larger ratio 
is used for the vertical component of ~0.57 gauss. 


IV. PERFORMANCE 


The performance of the automatically compensated degaussing system 
during the magnetic storm of March 5 and 6, 1961, is illustrated in Fig. 6. 
The upper part is a portion of the recording voltmeter chart for the vertical 
channel of the three-component magnetometer. The lower part is the corre- 
sponding chart record of the single-channel magnetometer whose fluxgate 
unit was oriented vertically in the null field region at the center of symmetry 
of the spectrometer. Note that the lower chart sensitivity is five times that of 
the upper. The maximum fluctuations in the lower lie within the desired limits 
of +10- gauss while the maximum variations in the vertical component 
caused by the magnetic storm are about 50 times this limit. 
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The fluctuations in the lower chart are due only in part to instabilities of 
either or both magnetometers. The periodic large spikes were caused when the 
spectrometer current was altered stepwise in an automatically programmed 
cycle. There is also a small periodic variation caused by the temperature 
cycling of the air-conditioning system which results in cyclic changes in the 
spectrometer coil spacings. This upsets the cancellation of fields from the 
spectrometer coils at the center of symmetry. As has been discussed in detail 
in our earlier paper (Graham et al. 1960) these dimensional changes do not, 
however, influence the spectrometer calibration, i.e., product of field strength 
at the optic circle multiplied by the optic circle radius. The central magne- 
tometer is thus, indirectly, a fairly sensitive monitor for fluctuations in the 
spectrometer room temperature as well as for those in the degaussing system 
and the spectrometer current regulating system. 

Experience with the automatic degaussing system over the past five months 
has been rewarding. It has not been necessary to reject any data because of 
magnetic storms as was the case previously. It is now found that the position 
(spectrometer current setting) for a given conversion line from a particular 
source has a typical day-to-day reproducibility of <1 part in 10°, i.e., that 
expected from counting statistics and other causes of uncertainty. 
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THE TEMPERATURE DEPENDENCE OF THE SHAPE OF 
PARAMAGNETIC RESONANCE LINES. PART IP 







MatcotmM McMILLAN? AND W. OPECHOWSKI 







ABSTRACT 


At sufficiently low temperatures the temperature dependence of the shape of 
paramagnetic resonance lines of crystals is mainly due to the temperature 
dependence of the occupation numbers of the energy levels of the system of 
interacting paramagnetic ions. The theory of this effect, recently discussed by 
the authors, is extended to the case of anisotropic exchange interaction and 
arbitrary shape of the paramagnetic crystal. 









In a recent paper (McMillan and Opechowski 1960, hereafter referred to as 
Part I*) we developed a theory of the shape of paramagnetic resonance lines. 
We assumed the effect of the variation of temperature on the resonance 
absorption to be the result of the temperature dependence of the average 
occupation numbers of the various energy levels of the absorbing system, 
thereby neglecting the effects connected with lattice vibrations, the “‘spin— 
lattice relaxation’, and thermal expansion, which assumptions should be 
valid at sufficiently low temperatures. 

In Part I we were concerned mainly with understanding and explaining the 
relationships between the theoretical expressions for the moments of an 
absorption line given by Van Vleck (1948), Pryce and Stevens (1950), and 
Kambe and Usui (1952). In addition we derived quite general but applicable 
expressions for the first and second moments of a line, and as an application 
we specialized them, for effective spins } and 1, to the case when the spin-spin 
interaction is the sum of magnetic dipole interaction and isotropic exchange 
interaction.f To give some simple numerical results we then considered the 
case of a spherically molded sample of nickel fluosilicate crystal. 

In this paper we shall extend the work of Part I to include an anisotropic 
exchange interaction, and as will be seen, the resulting expressions for the 
moments will, in general, be somewhat different from those given before. Again 























‘Manuscript received June 27, 1961. 
Contribution from the Department of Physics, University of British Columbia, Vancouver 8, 
B.C. 
2Present address: Department of Physics, Massachusetts Institute of Technology, 
Cambridge, Mass. 

*We have noticed the following errors in Part I: 

(i) The last line of equation (5.10) should read 


+ 2 (61 — 523) (be2 — 573) Aj By. 


(ii) The numerical values of the ordinates in Fig. 2 should be multiplied by 10. 

tWhen writing Part I we were not aware of the existence of a paper by Kopvillem (1959) 
in which he treats the case of effective spin } starting out from the general formulae of Pryce 
and Stevens. Kopvillem does not discuss the first moment at all. On the other hand, he gives 
an explicit formula for the second moment for the case of ethylsulphates and double nitrates 
of rare earths assuming isotropic exchange interaction. He also calculates the numerical values 
of the second moment for some of these substances taking T = 2.5° K and the resonance 
wavelength A = 3.25 cm, and compares these values with the scanty experimental data. 
We should like to thank Dr. Kopvillem for drawing our attention to his paper. 
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for purposes of illustration, the case of nickel fluosilicate will be considered 
but with the restriction on sample shape removed. It will be seen that using 
the results of this paper the signs and magnitudes of the exchange constants 
for this salt can be determined. 

We shall use the same notation and terminology as in Part I so only the 
necessary minimum of explanation of symbols and ideas will be given here. 

We are concerned with calculating the first and second moments of a line 
shape function which we defined in Part I to be proportional to the imaginary 
part of the high-frequency susceptibility. In particular, we are concerned with 
applications of the general equations (4.1) and (4.2) of Part I (which we called 
the ‘‘first approximation formulae”) to the case of a spin system wherein the 
spin-spin interaction is the sum of magnetic dipole interaction plus anisotropic 
exchange interaction. Thus, we take the truncated spin-spin interaction 
operator to be 


(1) Sf 1 7. HW o) 
ij 
(2) OY 2 AChS. Bist Bich. ptCiS.as 


where 
(3) CY = al +4216; Brvis— 1), 


(4) Ci = af —gi6'r; Byij— 1). 

The bar over the operator # denotes that part of WM which commutes 
with the unperturbed energy, the eigenstates of which are assumed also to be 
eigenstates of the operator corresponding to the z-component of spin. For 
nickel fluosilicate the z-direction will be taken to be in the direction of the 
optic axis of the crystal. In equations (3) and (4) a‘! and a’ are real numbers 
denoting the exchange constants perpendicular and parallel to the z-axis. 
There are only two such constants for the general anisotropic interaction 
because of the requirements of hermiticity and invariance under interchange 
of two spins. For isotropic exchange, 


(5) 
and the truncated spin-spin interaction operator above then reduces to 
equation (5.2) in Part I. 

Following the same procedure as was outlined in Section 5 of Part I but 
now using the interaction operator given by equation (1) above, we find that 
for effective spin S = 3 the first and second moments of the absorption line are 


(6) h (Av), as tanh( ) >, (c¥—ct), 
oe tj 


ats é eee sij ij\2 
h’(Av"), 4N = (Cy—Cy’) 
] . 


+ tanh’ = \ricy- cyYy(ce#-Cci yt , 


ijk 
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where JN is the number of spins in the system and where a; and az are the 
energy values of an unperturbed spin, that is, they are the energy values when 
the spin-spin interactions are “switched off’’. We label the levels so that 
a; < @2. 

It will be recalled from Part I that we use the notation (Av’) to denote the 
Ith moment of the line shape function about the unperturbed frequency which 
can be taken to be the position of the peak of the absorption line for high 
temperatures (for example, room temperature). It should be noted that 
(Av) is not, in general, equal to the /th central moment of the line shape function. 
We denote central moments by (Av’), and in particular we have the relationship 


(8) (Av*), = (Av?) — ({Ar))?. 


For effective spin S = 3 the second central moment can be written as 


(9) h(Av?)e, = i 1-21) Lci-ci* 


For effective spin S = 1 the first, second, and second central moments of the 
‘‘ground state line’”’ are, respectively, 


es ij / 


(11) h*(Av*), = 
(NZ) 2 { (Ci)? (6 +b" +b") —2C PC jb" + (C7)? (6 +6") } 
+(NZ’)* So’ (cC#(b" —b")?—2CPC} (6 — b%) (67 —b") 
ijk 
+C7C7F (6" —b")’}, 
h? Av"), = =| (ci4)*[1—2Z77(6" —b%)] 
J 
—2CC#[Z-b" — Z-* (6 —b™) (0 —b™)] 
te (C¥)*[Z-* (6 +6") —Z7(6" —5**)*}} ‘ 
and the first, second, and second central moments of the ‘‘excited state line”’ 


are, respectively, 
(13) h(Av); = (vz) (o—6") Dd Ci— (6"—") ricy! 
ij Wy ¢ 


(14) A*Ay*), = 
(NZ)~* SO" (CY)*O" +b" +b") —2C7C fb" + (CY) O" +0") | 


+ (NZEY* LECCE (Ob)? 2CLCF(O" —b") (67 —B") 
tj 


+C¥C¥ (6% —5")*}, 





1372 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


(15) WX Av*)ey = QO'E (CL) ZO" —5)"} 


J 
—2CCY[Z- "6" — Z-*(6" —b") (6% —b")] 
+4. (cyyiz a +b") = za —b**)?] } : 


By “ground state line’’ we mean the absorption line corresponding to transi- 
tions from the ground state to the first excited state of an unperturbed spin, 
and by “excited state line’? we mean the absorption line corresponding to 
transitions from the first excited state to the second excited state of an 
unperturbed spin. 

In the latter six equations we use the following notation for Boltzmann 
factors which we introduced in Part I: 


b* om en tl/kT 


We also use @), d2, @3 to mean the energy values of an unperturbed spin and 
we take a; < ad. < a3. Further, Z = 6% +6" +6%3. 

From this point on, we shall be concerned mainly with equations (10) and 
(13). In order to make these equations more amenable to analysis of experi- 
mental results we rewrite them in a different form. We define quantities 
x, y, and 6 by 


(16) x=N" DCH, 
ij 

(17) a > ee, 
ij 


(18) @ = arc tan (y/x). 
Then if we define F,, and F,, by 
bac 
FF... = ie ie for the ground state line, 
_h{Av) _ 


— 5 for the excited state line, 
V/ (x -+-9 ) 


Fex = 


we can write equations (10) and (13) in the form 
(19) Fy = Z7(b% —b®™) sin 6—Z—'(b% —b%) cos 8, 
(20) Fx = Z-(b™ —b%) sin @A—Z—!(b™ —b%) cos 8. 


It is convenient also to rewrite x and y in terms of more familiar quantities. 
Using equations (3) and (4) and making the usual assumption that the 
exchange constants are different from zero only between nearest neighbors 
we have 


(21) x=N* DI’ Cy = LCi = zay—24 6 Me, 


Jj 
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(22) y=N* Dd 'cCY = DY’ CH = 2a, 4+-3256°M1, 
ij 


J 






where z is the number of nearest neighbors and where (see Van Vleck 1937) 


Dts Bvis— 1) = N;®; 
j 





N, is the number of spins per unit volume and @ is the demagnetizing factor, 
a quantity which depends on the shape of the sample. In particular @ is zero 


for a spherically molded sample. 
We now consider as a point of illustration the application of equations (19) 


and (20) to nickel fluosilicate. The unperturbed energy values are 







a; = —g\BH—6, a, = 0, a3; = g\BH—4, 





where H is the applied magnetic field and 6 is a positive constant which, for the 
temperatures we are interested in, can be taken to be temperature independent. 
Since we have introduced the convention that a; < a2 < a3, we assume that 
gyBH > Os 

At this stage we must make one point clear. In Part I and in this paper we 
have assumed that when the resonance experiments are performed the applied 
magnetic field H is kept constant throughout the experiment and that the 
frequency of the oscillating magnetic field is varied. When this is the case the 
unperturbed energy values do not change since they depend on H but not on 
the frequency of the oscillating magnetic field, and therefore the only quantity 
occurring in equations (19) and (20) that changes during the experiment is 
the temperature 7. For the sake of convenience, however, the experiments 
are performed by fixing the frequency of the oscillating field and slowly varying 
the ‘‘constant” magnetic field H. In this case the unperturbed energy values 
change, this leading to difficulties in using (19) and (20). In order to use these 
equations, however, we shall assume that it is a good approximation to treat 
the unperturbed energy values as being constant. This assumption is certainly 
questionable when the lines are broad or when there is a large shift in the 
position of the lines with the temperature. 

Let us suppose that the frequency of the oscillating magnetic field is fixed 
and equal to »*. Then for the ground state line we shall take 























ad2—aQ, = hv* 






that is 







gy8BH+6 = hv* 


so the quantity H can be eliminated from equation (19). Thus we can write 
equation (19) as 







(23) Fer = Zz'(1—B"”) sin 9B—Z71(1—b"”"-”) cos 8, 





where 





Z1 = 140" +70", 
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For the excited state line we shall take 


a3—a. = hi*, 


that is, 


gi8H—6 = hy* 


so that quantity H can be eliminated from equation (20), which then can be 
written 


(24) Fen = Z73(b"" +? — 9?" *) sin 9— Z31(1—b°*" ™) cos 8, 


where 


Z. = 1 pre +25 4 p2cnr* +8) 


Thus equations (23) and (24) are the theoretical expressions for the first 
moments of the ground state line and the excited state line for nickel fluo- 
silicate when the frequency of the oscillating magnetic field is fixed and equal 
to v*. 

In Figs. 1 and 2 we have plotted these functions for various values of the 
parameter 6 for the case when v* = 24 kMc/sec. To illustrate the change in 


100 
Ry uo 
80 
Y= 24 kMc/sec 


0 


Y=24 kMc/sec 


6=135 


SS 
ie 


@=90° 


! 
TeK)~ 3 4 


-20 945° 
- 40 
-60 
-80 
-100~7 -100 


Fic. 1. The reduced first moment F,, of the ground state line as a function of the tempera- 
ture for various values of the parameter 6 for the case where v* = 24 kMc/sec. 
Fic. 2. The reduced first moment F., of the excited state line as a function of the tempera- 


ture for various values of the parameter 6 for the case where v* = 24 kMc/sec. 


the functions when another frequency is used we have plotted them for 
v* = 48 kMc/sec in Figs. 3 and 4. In all cases we have taken 6 to be tempera- 
ture independent and equal to the value of 0.12 cm which the results of 
Penrose and Stevens (1950) indicate for temperatures below 20° K. We have 
only considered values of 6 up to mw since the /’s are odd functions about 
6= nf. 

Now if experiments were performed using a frequency of 24 kMc/sec, say, 
then the angle @ and the normalizing factor (x?+y*)? could be determined by 
fitting the experimental results with the curves in Figs. 1 and 2. Hence the 
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y = 48 kMc/sec 





1——e- 
TK) 













; rns ~.20 @-45° 
leer. 

- 40 —6=30° we - 40 4 

he / 
-60 / 020° -60 / @=0° 

/ / 

3 od 

~.80 -.80 / 4 
“ie -100-—~ 









Fic. 3. The reduced first moment F,, of the ground state line as a function of the tempera- 
ture for various values of the parameter @ for the case where v* = 48 kMc/sec. 

Fic. 4. The reduced first moment F,x of the excited state line as a function of the tempera- 
ture for various values of the parameter @ for the case where v* = 48 kMc/sec. 











quantities x and y given in equations (16) and (17) could be determined, and 
thus using equations (21) and (22) the exchange constants a, and a, could 





be determined. 
As a further check on the values of the exchange constants found, the 
expressions for the second central moments of the two absorption lines could 






be used. 

Before concluding this paper we would like to emphasize a few points about 
our theory which must be realized before it can be used to analyze experi- 
mental data. Firstly, our theory applies only to the case when there is not too 
large an overlap in the absorption lines, this point being contained in the 
assumption leading to equations (2.10) and (2.11) of Part IL. Secondly, we 
have calculated expressions for the first and second moments of an absorption 
line and these quantities should not be confused with the “‘mean shift’ and 
the “half width” of a line. The mean shift of a line is the shift in the position 
of the peak, so only if the absorption line remains symmetrical about its 
maximum can the mean shift of a line from its position at room temperature 
be taken equal to the first moment. On the other hand there seems to exist no 
simple relationship between the half width and the second moment. In fact, 
if the line shape were described by a cutoff Lorentz curve it is quite con- 
ceivable to have the half width and second moment change in opposite 

















directions with a change in temperature. 
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NEW SPECTRUM OF CF+! 


R. D. VERMA 







ABSTRACT 

A spectrum (A'II-X'!+) of CF*+ has been obtained for the first time in the 
region 2025-1900 A. The rotational and the vibrational analyses are given. The 
principal constants of the two states as evaluated are given below in cm™, except 
for r. which is given in angstroms. 









> We WeXe B, Qe le 
1yp+ 0 1402.18 11.807 1.5172 0.01877 1.2290 
NT 51157 .08 1266.77 13.73 1.4239 0.01945 1.2686 






A brief discussion of the electron configurations is given. 






INTRODUCTION 


Recently (Verma and Mulliken 1961; Verma 1961) in work on the spectra 
of CCl and CCI* it was noticed that the ground state ro values of the CCI and 
CCI+ molecules are almost equal. This suggested that the vz electron in the 
electron configuration of CCl is nearly non-bonding. On the other hand our 
experience with the spectra of NO and NO*, which should be quite analogous 
to CCl and CClI* in their electron configurations, shows that the vz electron is 
strongly anti-bonding. Before any theoretical conclusion is drawn regarding 
this change of bonding character of the vm electron, it is desirable to have 
analogous experimental results on CF, SiF, SiCl, NS, and PO which fall in 
the same group of molecules as CCl and NO. The successful identification 
of a CCl* spectrum led to a search for the analogous spectrum of CF* which 
ended with the positive results that form the subject matter of the present 















paper. 







EXPERIMENTAL 






The spectrum was excited in a flowing atmosphere of helium containing a 
trace of CF,, at a few millimeters pressure, by a high-frequency discharge 
(2450 Mc/sec). This source emitted the new spectrum very strongly along 
with the B—X system (Andrews and Barrow 1951) of the CF molecule. On 
increasing the relative pressure of Cl’, in the mixture of He and CF,, it was 
noticed that the new spectrum at about 1960 A was only very weakly excited, 
whereas the B—X system of CF molecule remained unaffected. A photograph 
of the spectrum in the region of 1900-2100 A was taken with a 21-ft vacuum 
spectrograph in the third order of a grating with 30,000 lines/inch. The over- 
lapping second-order spectrum was absorbed by a chlorine filter (5-cm cell 
filled with chlorine at 1 atmosphere). Eastman Kodak SWR and Ilford Q, 
photographic plates were used. The exposure times ranged from 20 minutes 
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TABLE I 








J R 





(a) The (0,0) band, Ry = 51109.58 cm 


51091 .68 | 17 51060.23 —51012.33 
94.27 | 18 56.88 06.10 
96.80  51088.40 | 19 53.31 50999.75 
99.04 87.73 51079.49 || 20 51108.6 49.48 93.24 
101.12 86.96 75.76 21 07.5: 45.56 86.47 
03.04 86.07 1: 1} 2 06.: 41.48 79.59 
04.75 84.92 67.96 23 04. 37.25 72.47 
06.24 83. 66 63.47 || 24 03. 32.74 65.11 
07.53 82.19 59.62 || 25 O1.1: 28.05 57.73 
08.65 80.49 55.09 || 26 23.14 50.05 
09.58 78.62 50.32 || 27 18.20 42.25 

76.54 45.56 || 28 12.67 34.08 
74.31 40.51 29 07.40 25.93 
71.87 35.21 || 30 01.76 
69.21 29.72 || 31 50995 . 92 
66.37 24.07 || 32 90.02 
63.47 18.20 || 33 83.78 


Ooo RWNK OC WHO ON eK © 


(6) The (1,1) band, R, = 50971.32 cm 


50955. 21 | 15 
57 81 95 | 16 
59.56 34 ae 
61.49 56 | 18 
63.37 5.65 50932 19 
65.14 5.52 28. 20 
66.38 44.21 24 21 
67.86 2.48 22 
68.95 02 5.87 || 23 
69.91 39.12 27 (|| 24 
70.85 37.11 5.28 25 

53 || 26 
71.32 3: 96.17 || 27 
76 || 28 


] 
2 
3 
4 
5 
6 
7 
8 


) The (2,2) band, R, = 50826.61 cm™ 
50801 13 50788 
01.3: 14 85 
798 15 83 
7 ie 79. 
95.36 17 16 
93 18 73 
90.92 19 69 


The (1,0) band, R, = 52346.66 cm™! 


— 


Nowe Goss! 


Ot Or OO te 


S 


) The (0,1) band, R, = 49737.89 cm7 


15 49692 | * 49649 
49709 16 8Y 8 Ad 
09 49700 29 17 49737 89 87.3: 39. 2: 
O8 697 3: 18 84 33 
07 93.7 19 81.8: 28 
07 8Y 89 20 78 22.4 
05 86.1! 2) 75 16: 


TS8S. OR ON eS 














VERMA: SPECTRUM OF CFt 


TABLE I (Concluded) 











Q P J R Q P 

. 30.22 04.77 81.86 22 72.36 10.31 
9 31.60 03.45 78.00 || 23 69.05 03.98 
10 33.02 01.95 73.50 || 24 65.29 597.85 
il 34.17 700.29 69.05 || 25 61.57 90.74 
12 35.22 698.54 64.62 || 26 57.41 83.97 
13 36.05 96.58 59.85 27 53.59 

14 36.82 94.49 54.93 | 28 49.65 


(f) The (1,2) band, R, = 49620.71 cm= 


2 49601.85 15 49575.89 49534.12 
3 03.98  49592.87 | 16  49620.71 73.48 28.88 
4 06.29 92.47 | 17 70.87 23 56 
5 08.34 91.58 18 68.20 18.06 
6 10.31 90.74 49573.91 | 19 65.43 12.48 
7 12.02 89.68 70.23 | 20 62.23 06.67 
8 13.48 88.45 66.05 || 21 59.18 500.74 
9 14.90 87.12 62.23 | 22 55.89 494.77 
10 16.30 85.63 57.71 || 23 52.45 88.56 
rr 17.43 83.96 53.32 | 24 48.78 2.19 
12 18.34 82.17 48.78 25 45.05 75.72 
13 19.15 80.23 43.98 | 26 41.07 
| 


19. 78.06 39.17 27 37.03 








(g) The (2,3) band, R, = 49498.82 cm7 


0 49475 .68 14 49457 .02 49418. 56 
1 77.81 15 54.74 13.53 
2 80.67 49472 .55 16 49498 . 82 52.25 08.38 
3 83.10 72.08 17 49.66 02.99 
4 85.20 71.47 18 46.90 397.55 
5 87.11 70.71 49457 .06 19 44.07 91.91 
6 88.98 69.82 53.27 20 41.09 

7 90.67 68.74 49.31 || 21 37.79 

8 92.23 67.48 45.47 22 34.52 

9 93.68 66.12 41.34 23 30.85 

10 94.76 64.58 37.13 24 27.26 

11 95.81 62.95 32.90 25 23.45 

12 96.74 61.10 28.16 

13 92.55 59.14 23.45 


i 
to 3 hours-for the Av = 0 and —1 sequences while an exposure of 6 hours for 
the (1,0) band was not sufficient for the development of the whole band. The 
second-order hollow-cathode iron spectrum was superimposed for the standard 
spectrum. The hollow-cathode iron wavelengths were taken from the recent 
(Crosswhite 1958) interferometric measurements. Three different plates were 
measured on a photoelectric comparator (Tomkins and Fred 1951) and an 
ordinary comparator (Swiss, SIP). The wavelengths were evaluated by using 
a linear interpolation formula and the corrections obtained from the second- 
order iron wavelengths were adjusted so as to give correctly the third-order 
Si lines which were excited in the CF sources. This must have taken care of 
the coincidence error if it were there at all. The accuracy of the measurements 
was estimated as better than +2 ma (equivalent to +0.05 cm~'!). The vacuum 
wavenumbers of the lines of all the bands are shown in Table | along with their 


analysis. 
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ROTATIONAL ANALYSIS 


The spectrum consists of two well-developed sequences which could be 
identified easily from their intensity as Av = 0 and —1 containing (0,0), (1,1), 
(2,2) and (0,1), (1,2), (2,3) bands respectively. The (0,0) and (0,1) bands 
were quite strong. The (1,0) band was very weak and only a few Q lines 
could be recorded in the longest exposure given. The R branches of the bands 
show pronounced head formation. As is evident from Fig. 1, the spectrum 
was well resolved, and no difficulty was encountered in identifying the various 
branches. Each band was found to have an intense Q branch along with R 
and P branches. In the (0,0) band, the R(Q) line was identified without any 
ambiguity. One might doubt that the R(O) line shown in Fig. 1 might as well 
be R(29). However, R(29) is expected to be very weak because the corre- 
sponding line in the P series is not recorded. Further, in a number of lightly 
exposed plates in which P lines of high J value were not developed, the R(0) 
line could be seen. Hence it was concluded that the spectrum is due to a 
‘TI'S transition. This is also consistent with the fact that in the 'lI-'Z tran- 
sition the R branch should have higher intensity than the P branch for a low 
J value. The emitter of this spectrum can safely be stated to be CFt, taking 
into account the structure of the bands, the experimental conditions, and the 
magnitudes of the constants. 

The assignment of J numbering in the analysis of the bands could be made 
without any difficulty and it was checked with the well-known tests described 
by Herzberg (1959). The determination of the constants was performed by 
standard methods (Herzberg 1950). However, since the R branches were not 
resolved for high J values, the combination differences A,/ could not be formed 
for all the J values. Therefore, an evaluation of the B, and D, values could 
be effected from the combination differences 4,F obtained from the Q and 
P lines. 

One might, of course, expect that A-splitting in the "II state might cause a 
slight error in such a determination of B, values. However, it was possible to 
examine the A-splitting up to J = 14 in some of the bands by comparing the 
differences R(J) —Q(J) and O(J +1) — P(J +1). The results showed no evidence 
for A-splitting up to J = 14. Furthermore, in the (0,0) band, the returning R 


g 
branch seems to overlap the J lines of lower J value. If this is true, which 
obviously it is, one can verify that there is no A-splitting in the II state within 
the experimental error, by extending the above checkup to the highest J 
level observed. It is thus certain that the A-splitting in the 4II state is less 
than the experimental error. Similar results are shown by the author in the 
case of CCI*. The cause of the absence of this splitting was discussed at length 
in the case of CCl* and the same considerations can be applied to the case 


$y was evaluated by averaging the sets of 4,7’ values obtained from the 
(0,0) and (0,1) bands, the lines of which were measured more accurately than 
others due to their higher intensities. The remaining B, values (B,’, B,’, and 


b,'’, v = 0-3) were then determined by using only Q lines. If one plots 





VY 0002 4v9t 19 Jo wiaisss purg TO 


(I‘I)d (I'l) (O°O)0 (0*O)Y 


] ALVId 
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Q(J) —CJ(J+1) against J/(J+1), where C is an approximate value of (B’— B”’), 
the slope of the straight line obtained gives the necessary correction to obtain 
the exact value of (B’—B’’) and the intercept with the ordinate gives the 
vo of the band. The difference (B’ — B’’) can thus be determined very accurately. 
If one of the B’s is known accurately by another method, the other can be 
determined in this way in the cases where the measurements of the R and 
P lines are not good or their series are not complete. In the present case, all 
the remaining B, values were determined in the above manner relative to 
B,’ from the various bands. The B, values so evaluated for all the bands are 
collected in Table II and are represented by the following formulae for the 
two states: 


Il 


(1) B,’ = 1.4239—0.01945(0+4), 


(2) B,” = 1,5172—0.01877(v+3). 


TABLE II 
The rotational constants 
B, for ‘II and '* states of CFt 





v B,’ 


B ur 
v 





0 1.4142 1.5078 
l 1.3949 1.4888 
2 1.3753 1.4701 
3 1.4515 


All the principal rotational constants calculated are shown in Table III. 


TABLE III 


Principal constants of ‘Il and !Z* states of CF* 


B,’ = 1.4239 cm B,”’ = 1.5172 cm™ 
a,’ = 0.01945 cm7 ae” = 0.01877 cm" 
D,’ = 6.1X107* cm D,.”” = 6.7 X10-* cm™ 
re’ = 1.2686 A re’ = 1.22900 A 
we = 1266.77 cm™ we’ = 1402.18 cm™ 
wx, = 13.73 cm“ we Xe’ = 11.807 cm™ 
ve = 51157.08 cm™ weve’ = +0.0517 cm™ 


VIBRATIONAL ANALYSIS 


The band origins of all the bands determined from the Q branches in the 
manner described above are arranged in a Deslandres scheme shown in 
Table IV. The following formula represents their values very well: 


(3) » = 51157.08+[1266.77(v’ +4) —13.73(v—3)"] 


— [1402.18 (0 +4) — 11.807(v” +4) +0.0517(0” +4)3]. 
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TABLE IV 


The Deslandres scheme for the band origins 


at! 
w 


0 1 
51088.89 1378.7 49710.16 
1239.31 1239.31 
52328 .20 1378.73 50949 . 47 1355.58 49593 .89 
1211.85 
50805. 74 


1332.74 49473 .00 


ELECTRON CONFIGURATION 
The present spectrum is quite analogous to that of the CCl* molecule 
described recently. Therefore all the remarks about the electron configurations 
and the electronic states given in the case of CCI* are applicable to CF* also 
and hence it is very probable that the lower state is the ground state of the 
molecule. Then the electronic transition can be written as 


wri xovm, A Il > wr xo”, X St. 


It is interesting to note that r, decreases on ionization and increases on 
excitation of the vz electron appreciably in the case of CF (see Table V), as 


TABLE V4 


CF or_CF* CCl.or:CCl* 


Electron configuration re (A) ro (A) 


Wes ea 
wr xo* vr 


CF 653 CCI 
pee pacbina CF* 634 CCI* 


I 
I 
wr xo? 2290 CF* 685 CCI* 
wri xo ur . 2686 CFt* 634 CCI** 


“Note: The B-X system of CF (Andrews and Barrow 1951) has been 


photegraphed at high resolution by the author. The rotational structure of 


the spectrum indicates that B is probably a 22 state and perhaps has the 


electron configuration wz‘ xo vzx?. 


it does in NO, whereas in the case of CCI] the changes in 7» are small and 
of opposite sign. This suggests that the um electron is definitely anti-bonding 
in the case of CF and NO, but becomes slightly bonding in CCI. One also 
finds from the comparison given in Table V that r, or ry of the wz! xo? ur 
electron configuration of CF or CCl and that of the wx xo um electron con- 
figuration of CF* or CCI* are almost equal, indicating that the xo electron 
in both CF and CC] is nearly non-bonding, while it is fairly bonding in the 
case of NO. 
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THE REVERSIBLE TEMPERATURE COEFFICIENT OF 
REMANENCE OF A Pt-Co PERMANENT MAGNET* 


W. L. ZINGERY AND T. M. WIRT 


The reversible temperature coefficient of remanence of a Pt-Co permanent 
magnet has been measured. Experimental data on this outstanding per- 
manent magnet alloy is meager because of its relatively higher cost. Martin 
(1957) has studied the metallurgy and optimum heat treatment for this 
material and the Johnson, Matthey and Company (1959) have made measure- 
ments on the use of this material at 100° C. 

Because Pt—Co has the highest coercive force known and would seem to 
offer unusual magnetic stability, a 50:50 atomic per cent multipole magnet 
assembly of Pt-Co was measured with great precision (Zingery 1961) for 
long-term magnetic stability (Zingery et al., to be published) and for the 
reversible temperature coefficient of remanence. The magnets were found to 
have superior magnetic stability and to have a reversible temperature co- 
efficient of —0.042% /°C. This is about twice the value obtained for a similar 
Alnico V assembly. The magnets were operating at a permeance coefficient 
of 20, and were temperature-cycled five times from 100° C to —60° C before 
the reversible measurements were made. The Johnson, Matthey and Company 


1800+ =A INCREASING TEMPERATURE 
© DECREASING TEMPERATURE 


THERMISTOR RESISTANCE 


a 
1300'\ 6 


9.86 9.87 9.88 9.89 
RESONANT FREQUENCY (Mc) 


Reversible magnetic field change with a change in magnet temperature, Resonant 
¥ 6 > 


y is proportional to magnetic field. Thermistor was a type 3$2A1 and was oven-aged 


s work was supported by the U.S. Air Force. 


sy 496) 











NOTES 1385 












reports that their Pt-Co (Platinax II) magnetization suffers a decrease of 
“about 3.5% when first heated to 100° C’’. If one assumes an initial tempera- 
ture of 20°C and a linear coefficient of —0.042%/°C, an increase in tem- 
perature to 100° C would yield a 3.36% decrease in magnetization. This leaves 
very little irreversible change. Figure 1 illustrates the reversible nature of our 
measurements. 
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SPECIFIC HEAT ANOMALIES ASSOCIATED WITH THE 
RESISTANCE-MINIMUM PHENOMENA IN 
MAGNESIUM-MANGANESE ALLOYS* 










DovucGLas L. MARTIN 






Small amounts of transition metal impurities in solution in some pure 
metals cause an anomalous variation of electrical resistance with temperature 
giant’’t thermoelectric 






oe 


at low temperatures. Usually associated with this are 
powers and large specific heat anomalies (De Nobel and Du Chatenier 1959; 
Zimmerman and Hoare 1960; Martin and Franck 1960; Crane and Zimmerman 
1961; Franck, Manchester, and Martin 1961). It has been known for some 
time that impure magnesium shows a resistance minimum at low tempera- 
tures (Meissner and Voigt 1930) but previous attempts to detect a corre- 
sponding specific heat anomaly led to a negative result (Logan, Clement, and 
Jeffers 1957). 

Recently Hedgecock, Muir, and Wallingford (1960) and Gaudet, Hedgecock, 
Lamarche, and Wallingford (1960) have published a series of measurements 
on the electrical resistance of dilute alloys of manganese in magnesium. They 
showed that for very dilute alloys, only a resistance minimum was observed 














in measurements down to about 2° K. For slightly more concentrated alloys 






(greater than about 0.1 atomic ©% manganese) a maximum was also observed. 
; ( x 






Since previous attempts to detect a specific heat anomaly in this alloy 





system might have failed because measurements were made at too high 






temperatures, a preliminary investigation has been made in the lowest tem- 






perature range at our disposal (O.4° to 1.5° K). 






*Issued as N.R.C, No, 6451. 
TThis terminology is due to Dr. M. Bailyn. 
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The results of measurements on two alloys are given in this note. The 
calorimeter is described elsewhere (Martin 196la). The first alloy was No. 
729 of Hedgcock et al. (1960) and contained nominally 0.018 atomic % of 
manganese. The alloy sample and strips rolled from the top and bottom were 
annealed for 16 hours at 450° C in a helium atmosphere. The ratios of electrical 
resistance at 4.2°K to that at 273° K for the test strips were 0.025; and 
0.025,, indicating rather good homogeneity. From these figures the composi- 
tion is estimated at 0.025 atomic % Mn.* This sample weighed 44.741 g and 
the atomic weight was taken as 24.328. The second alloy was prepared by 
melting a magnesium—manganese alloy containing 1.5 wt. % Mn (from the 
Dow Chemical Co.) with Johnson-Matthey spectroscopically pure magne- 
sium, in a graphite crucible under a helium atmosphere, in order to produce 
an alloy with composition of a nominal 0.15 atomic ©Z manganese. The alloy 
was cooled slowly after melting and then annealed in a helium atmosphere 
for 11 days at 630° C in an attempt to homogenize the ingot. A number of 
test strips were then rolled from various places on the ingot and the strips 
and ingot were annealed at 450°C for 24 hours and water-quenched. The 
ratio of electrical resistance at 4.2° K to that at 293° K for the various test 
strips varied from 0.116 to 0.235 indicating that some parts of the ingot 
contained twice as much manganese in solution as did others. Thus the 
sample was not very homogeneous. This sample weighed 19.600 ¢ and the 
average atomic weight of the alloy was taken as 24.366. 

The results of the specific heat measurements on the alloys (two independent 
sets of measurements on each) are shown in Fig. | together with a line repre- 


senting the specific heat of pure magnesium (Martin 19610). It will be noted 


that for the more concentrated alloy the maximum of the specific heat anomaly 
is probably in the same region as the maximum in the electrical resistance 
(Gaudet et al. 1960). For the more dilute alloy the maximum of the specific 
heat anomaly appears to lie below 0.4° K. It is therefore interesting to specu- 
late whether a maximum in the electrical resistance might be found in this 
temperature range also. Such evidence as is available (Hein and Falge 1957) 
would indicate that this is not the case. As far as can be decided from the 
limited temperature range of the measurements, both the anomalies are of 
similar shape to those found for the copper-iron system (Martin and Franck 
1960; Franck et al. 1961) which apparently shows only a resistance minimum 
and no maximum. Certainly the maximum of the specific heat anomaly in 
the copper-iron system occurs in a region where there is no maximum in the 
electrical resistance. It is clear that some work remains to be done before a 
satisfactory correlation between the resistive and thermal anomalies in these 
dilute alloys can be obtained. 

The temperature range of the present measurement is too small to permit 
an estimate of the entropy contents of the specific heat anomalies to be 
made. Jt does seem clear, however, that the spin of the manganese ions in 
solution in magnesium must be rather large, perhaps the same as for man- 
ganese in the monovalent metals copper and silver (De Nobel and Du Chatenier 
1959; Zimmerman and Hoare 1960; Collings and Hedgcock 1961), 

During the specific heat measurements on the more concentrated alloy it 


*! am indebted to Mr. Ek. Wallingford, of Ottawa University, for this estiniate 
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SPECIFIC HEAT OF MAGNESIUM AND MAGNESIUM—MANGANESE ALLOYS 
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Fic. 1. The specitic heats of pure magnesium and two dilute alloys of manganese in mag- 
nesium., 


was noted that a time of the order of 1 second was required for thermal 
equilibrium to become established within the sample after a heating period. 
This appeared to be a characteristic of the sample and not due to poor thermal 
contact between sample and calorimeter. The effect was most marked at the 
lowest temperatures and might correspond to the exchange of thermal energy 
between the free electrons and the spin system of the manganese ions. 


| am indebted to Dr, F. T. Hedgecock of Ottawa University for letting me 
have the alloy samples and to Messrs. J. W. Fisher and R. L. Snowdon for 


preparing the calorimeter samples. 
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THE ELECTRONIC SPECIFIC HEAT OF LITHIUM-MAGNESIUM ALLOYS* 


Dovuc.Las L. MarTIN 


Pure lithium metal has a body-centered cubic structure at room tempera- 
ture, and up to about 70 atomic % magnesium may be dissolved in lithium 
without altering this structure (Hansen 1958). On cooling to liquid helium 
temperatures, alloys containing less than about 20 atomic % magnesium 
undergo a partial transformation of the martensitic, or diffusionless, type to 
a faulted hexagonal lattice (Barrett and Trautz 1948; Barrett and Clifton 
1950). It seems probable that at 1° K most (~90%) of both the pure lithium 
and the alloys which transform is in the low-temperature phase (see Barrett 
and Trautz 1948; Martin 1960a, 6; Hull and Rosenberg 1960). Measurement 
of the electronic specific heat is therefore made on a two-phase system and 
since the electronic properties of the two phases may differ it is clear that 
detailed information cannot be expected from the present experiments. How- 
ever, some general features do emerge and the results are of some interest 
in connection with other recent work on these alloys (Dugdale and Gugan 
1961; MacDonald, Pearson, and Templeton 1961). 

The measurements were made in a *He cryostat using a technique described 
elsewhere (Martin 1961a). The alloys were made by heating the constituents 
(Lithium Corporation of America low sodium grade of lithium (purity 
~99.99%) and Johnson-Matthey spectroscopically pure magnesium (purity 
~99.99+%)) in a sealed iron container to 800° C for 1 hour. The 1 atomic % 
magnesium in lithium sample weighed 8.908+0.014 g (after buoyancy cor- 
rection) and the weight of oxide and other contaminations was estimated at 
less than 0.011 g. The average atomic weight of this alloy was taken as 
7.114. The ratio of electrical resistance at 4.2° K to that at 293° K was measured 
for pieces taken from each end of the sample and found to be 0.1208 and 
0.1192 respectively. Thus the sample appears to be rather homogeneous and 
the residual resistance ratio confirms that the composition is very close to 
1 atomic % magnesium (Dugdale, private communication). The 10 atomic % 
alloy weighed 11.42+0.03 g, the weight of oxide and other contaminations 
was estimated at less than 0.05 g, and the average atomic weight of the alloy 
was taken as 8.678. The residual electrical resistance ratios of pieces taken 
from each end of the sample were 0.558 and 0.545, indicating reasonable 
homogeneity and that all the magnesium was probably in solution. 

The results of the specific heat measurements are shown in Fig. 1. Two 
separate sets of measurements were made on each sample, the carbon thermo- 
meter being recalibrated each time it was cooled to low temperatures. The 
results have been analyzed by the method of least squares, assuming that 
they may be represented by the equation 

C= yp leek T 
Sb ¢ 
*Issued as N.R.C. No. 6450. 
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SPECIFIC HEAT OF LITHIUM-MAGNESIUM ALLOYS 
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Fic. 1. The specific heats of two lithium-magnesium alloys. 


where C is the total specific heat, y the electronic specific heat coefficient, 
R the gas constant, and @ the Debye temperature. The values found for y 
and @ are shown in Table I, together with the 95% confidence limits. 


TABLE | 


Results for lithium—magnesium alloys 


y (ucal/°K? g-atom) 6 (°K) 


Lithium — 1 atomic % magnesium 


20 @ aon +49 
Ist set 386.8+4 33073 
2nd set 386.845 3497 fn 
. ace 719 g9q +38 
All points 386.7+3.2 339 56 
Lithium — 10 atomic % magnesium 
os . roe t+ O7 
Ist set 363. 7+6 327 — 
—wo 
«AO 
2nd set 368. 8+6 379 b 220 
-bo 
; aan 0.0 6 249 +60 
All points 366.34+3.8 349 bo 


Note: Error limits are 95° confidence limits. 


In Table I] the results are compared with recent values of y and @ for pure 
lithium (Martin 196la@) and magnesium (Martin 19616). The Sommerfeld 
free electron values for y for the pure metals are also shown. 

The main feature of the results is the smooth variation of both y and @ 


between the values for the pure metals. If it is assumed that the lithium 
and alloys have almost completely transformed to the hexagonal phase, then 
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TABLE II 


Summary of results for pure lithium and magnesium and lithium— 
magnesium alloys 


y (ucal/°K? g-atom) 6 (°K) 
Pure lithium ie 
(Martin 1961) 38943 a6 to 
Lithium —- 1 atomic % +38 
magnesium 386.7+3.2 339 156 
Lithium — 10 atomic % +60 
magnesium 366.343.8 349 
Pure magnesium +110 
(Martin 19616) 292.7+42.5 440 | 56 


Note: Error limits are 95% confidence limits. 


Sommerfeld free electron values for 4+ 
Pure lithium 180 weal/°K? g-atom 
Pure magnesium 237 ucal/°K? g-atom 


all the figures in Table I] may be approximately interpreted as referring to 
hexagonal lattices (since magnesium has a hexagonal lattice). The lattice 
parameters of the hexagonal phase of pure lithium, the alloys, and pure 
magnesium are very similar (see Pearson 1958). Then the variation of y on 
going from pure lithium to pure magnesium must be due mainly to the 
variation in the average number of electrons per atom. The variation might 
be attributed to an alteration in the density of states at the Fermi level and 
to alterations in the effects of the electron-phonon and electron-electron inter- 
actions on the electronic specific heat. It should be noted that the observed 
electronic specific heat decreases with increase in electron density while the 
Sommerfeld free electron value increases with increase of electron density. 
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THE Xe'®s YIELD IN THERMAL FISSION OF Pu*°* 
J. G. Bayty,! M. F. Duret,' N. B. PouLsen,':? anp R. H. ToMLtinson® 


Experiments to measure the absolute yields of the products of Pu*** fission 
have been described by Fickel and Tomlinson (1959a). In their paper the 
xenon yield had been determined from two samples, a PuO. sample and a 
Pu-Al alloy sample, the quoted result being assigned an error of about 15% 
because of the uncertainty in estimating how much Xe! was removed from 
the 135 mass chain by neutron absorption during the irradiation. We have 
used the data on the Pu-Al sample and an analogue computer to make a 
more accurate estimate of the correction due to neutron absorption in Xe!. 
The Pu-Al sample was selected since the self-shielding factor was small and 
could be estimated reliably (Fritze, McMullen, and Thode 1959). 

The details of the reactor power history were obtained from the NRX 
Operation Log. These were corrected to agree with the NRX Power Log, a 
reliable record of the relative daily energy output. The history, in digital 
form punched on paper tape, was fed into an analogue computer which 
simulated the neutron flux and the iodine, xenon, and cesium differential 
equations. The neutron flux was assumed proportional to the reactor power 
level, the normalizing factor being determined from the total energy output 
of the reactor and the irradiation of the sample which had been determined 
from cobalt monitors by Fritze et a/. (1959). 

The equations describing the amount of iodine, xenon, and cesium in the 


sample as functions of time are: 


dl'/dt = —d, 1+ y 24, 
dXe™/dt = d,1' + yx. —Ax Xe! — oxgXe", 
dCs"5/dt = Ax Xe", 
dXe"6/dt = ox¢Xe"™, 


where the symbol Xe'* does not include that formed directly in fission. 

The xenon cross section, iodine yield, and xenon yield were taken from 
Fickel and Tomlinson (19596), who measured the xenon capture cross section 
in a neutron spectrum similar to that of the Pu-Al sample irradiation. The 
numerical values of constants appearing in the equations were as follows: 


ox = 3.2 10° barns, 

y, = 0.1084 per hour, 

Ax = 0.0752 per hour, 

*Issued as A.E.C.L. No. 1308. ; 
1Atomic Energy of Canada Limited, Chalk River, Ontario, 


*Present address: Texas Agricultural and Mechanical College, College Station, Texas, U.S.A. 
3McMaster University, Hamilton, Ontario. 
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yx/¥1 = 0.0536, 

total irradiation = 4.83X10!8 n/cm?, 
total time = 564.3 hours, 

total reactor energy = 15,750 Mw hours, 
ox@ = 9.82X10-4 P per hour, 

P = reactor power level, Mw. 


Figure 1 shows the reactor history and the relative concentrations of 
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Fic. 1. Reactor history and relative nuclide concentrations. 
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Xe!, Cs!5, [135 and Xe'*®, The results show that 28.9 atoms of Xe! are 
destroyed to form Xe'** for every 78.2 atoms which decay to Cs" so that the 
total Xe! produced was 37% greater than the measured Cs! yield. This 
result is insensitive to assumption concerning the direct yield of Xe'® in 
fission, being reduced by 0.3% if yx is reduced to zero. 

Combining these results with other results determined by Fickel and 
Tomlinson we can obtain the yield of Xe! in fission of Pu”. The absolute 
yield of Cs!*8 was determined to be 6.90% and the ratio y135/¥133 = 0.786+.011 
was determined mass spectrometrically. Thus these results lead to a Xe! 
yield of 0.786 X6.90% X 1.37 = 7.43% for this sample. 

Estimating possible contributing errors to be 3% in the Xe! cross section, 
3% in the irradiation and assuming errors of 0.3% because of the uncertainty 
in the direct yield of Xe! and a 1% computing precision, we conclude that 
our estimate of the amount of Xe! destroyed by neutron capture has an 
accuracy of 4.9%. This corresponds to an error in the Xe! yield of 2.7% 
relative to the Cs! yield. 


FIckEL, H. R. and Tomiinson, R. H. 1959a. Can. J. Phys. 37, 926. 
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national Conference, Geneva, September, 1958, 15, p. 436. 
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NOTE ON A MECHANISM FOR THE VERTICAL TRANSPORT 
i OF IONIZATION IN THE IONOSPHERE 


W. |. AxForD 


The object of this note is to point out a possible mechanism for the vertical 
transport of ionization in the ionosphere, which involves in turn the process 
of formation of dense ionized layers described by Dungey (1956) and the 
gravity waves of Hines (1960). The main feature of the mechanism is that it 
causes ionization to be carried down from a reservoir such as the E region, 
in a series of layers which may be quite dense and well defined even when the 
density of ionization in their immediate neighborhood is not large. This is 
essentially a convective process, and is therefore considerably more efficient 
in transporting ionization than turbulence, which is diffusive and merely 
broadens the existing layers. 
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In its simplest form, the Dungey process can be most easily explained by a 
consideration of Fig. 1. As a rough approximation, the electrons and ions can 





Fic. 1. The simplest form of Dungey’s process. Ionization is constrained to move on the 
lines of force of the magnetic field (dashed lines). The neutral wind profile shown therefore 
forces the ionization along the field lines in the direction of the small arrows, causing it to pile 
up around the dotted line. 


be regarded as being constrained to move on the lines of force of the magnetic 
field, like beads on a string; thus the component of the neutral wind parallel 
to the field drives the ionization along the field lines, while the perpendicular 
component simply distorts the field slightly. If the wind profile is of the form 
shown in the figure, it is obvious that the ionization must pile up at the point 
of zero velocity until its partial pressure gradient balances the force exerted 
by the neutral particles. This description is applicable to / region levels; in 
general the process is more complicated due to effects such as recombination 
and the increased mobility of ions transverse to the magnetic field in the 
lower regions of the ionosphere.* 

As a simple extension of Dungey’s process we see that vertical transport of 
ionization can occur if a node in the wind profile (such as that shown in Fig. 1) 
moves vertically, since the ionization is made to follow the node by the force 
exerted by the neutral gas. The vertical motion of such a node is a natural 
consequence of the presence of gravity waves in the upper atmosphere, as 
described by Hines; the phase velocity of the dominant waves is directed almost 
vertically, and in the downward direction if (as is expected) the energy is 
propagating upward. Hence if a wave moves downward from a region of 
ionized gas it will carry ionization with it as a series of layers situated at every 


*\ more detailed analysis is given by Whitehead (1961); the problem is also under study 
by L. R. O. Storey and his colleagues, and by the author. 
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second node. This may be called the ‘corkscrew’ effect, since ionization is 
drawn down from the ionosphere just as a corkscrew draws a cork from a 
bottle. Note that the separation of the layers must be an integral number of 
wavelengths, since they can only occur (by the simple Dungey process) at 
nodes where the neutral wind profile and the magnetic field are related in the 
manner shown in Fig. 1. 

Downward moving reflecting layers are often observed between the E and 
F regions, but it is not clear whether a significant amount of ionization is 
actually being transported, or whether the effect is purely wavelike with no 
net transport of ionization. Sharply defined sporadic-E layers of a similar 
nature are also commonly observed; here the vertical transport of ionization 
is likely to be more important, since a node moving down from the E region 
encounters a rapid decrease in the density of ionization produced in situ. 
Some ionization may be carried down to the D region in this way, but at these 
levels the efficiency of the process is greatly impaired. 

The tendency of the ionization to follow a moving node in the neutral wind 
profile is reduced by collisions with neutral particles, and for a given wind 
velocity there is a limit to the vertical velocity which can be imparted to the 
ionization; thus if a node changes its position too quickly the ionization will 
not be able to follow. Hence for any given phase speed of the wave, there is a 
critical height above which ionization will follow and below which it will not 
follow, a node in the wave pattern. For typical gravity waves the corkscrew 
effect appears to become inefficient in regions below 110 km or so. It is observed 
that sporadic-E layers are most dense on the average at this level and this, 
we suggest, is due to an accumulation of ionization at the point where it is no 
longer completely entrained by the waves. 

As a result of reflections at certain heights, standing waves can occur in 
addition to the travelling waves. For a pure standing wave, the nodes (and 
therefore the layers of ionization) occur at fixed heights which are determined 
by the structure of the atmosphere. If a downward travelling wave is combined 
with the standing wave, the resulting nodes move downward in a rather 
intermittent fashion, their motion being slowest in the vicinity of the nodes 
formed by the standing wave pattern above. This can be seen from Fig. 2 
where the variation with time (¢) of the height (#) of a particular node is 
shown for superimposed standing and travelling waves of equal wavelength 
(#1) and period (7°), and for a range of values of the amplitude ratio, R. Clearly, 
if we were to observe the instantaneous position of nodes in such a pattern 
(having a finite, non-zero value of R), it would appear that there were 
‘preferred heights’ corresponding to the position of the nodes of the standing 
wave component; this becomes more pronounced with increasing values of R. 
The effect is enhanced by the variations of the neutral wind shear at the nodal 
points; the shear is greatest when the nodes are moving slowly and so the 
ionized layers are most dense at the ‘preferred heights’. 

There is evidence that sporadic layers of ionization occur at ‘preferred 
heights’ in the £ region, with separations of about 10 km (Seddon 1960). 
This can be interpreted as a statistical effect caused by the superposition of 





CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


0 7 2T 


Fic. 2, Trajectories in the (4,t) plane of a node in a wave pattern consisting of super- 
imposed standing and downward travelling waves, for various values of the amplitude ratio R. 
When R = 0, the pattern is that of a pure travelling wave and the node moves steadily 
downward. When R = ©, the pattern is that of a pure standing wave, and the node remains 
at a constant height. 


standing and travelling Hines waves, as described above. Seasonal variations 
in the distribution of ‘preferred heights’ would then reflect nothing more 
than the variations in the structure of the atmosphere. 

Gregory (1961) has observed that D region scattering layers also seem to 
occur at preferred heights. The standing wave explanation is still in order for 
the D region, but the ‘corkscrew’ effect becomes very much weaker with 
decreasing altitude. Below 110 km the ionization is unlikely to follow com- 
pletely a particular node in the wave pattern, but a series of nodes moving 
downward must nevertheless entrain ionization to some extent. The mechanism 
may contribute a significant amount of ionization to the D region, since the 
total quantity that must be accounted for is relatively small. 
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